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0 INTRODUCTION 


Most of the past effort in the field of earth resources data processing 
has been research oriented* Earth resources imagery has been provided by NASA 
to a number of researchers who have processed the data in various ways in 
order to determine what, if any, useful information could be extracted from 
the given images* These experiments have demonstrated that useful information 
can indeed be extracted from aircraft and satellite multispectral scanner 
imagery of the earth’s surface* Economic studies have indicated potential cost 
effective systems based on these techniques. Consequently, it is anticipated 
that during the 1980-1990 decade earth. resources satellites will be designed 
and flown for specific purposes, i.e*, to monitor severe weather systems, to 
monitor water pollution, to survey and monitor world food production, etc* In. 
these applications it may be more cost effective to process the data on-board 
the satellite and transmit the data products directly to the users rather than 
transmit the raw data to a ground processing station for generating the data 
products and then distributing the data products to the users via another 
satellite system. 

The purpose of this study was to investigate the feasibility of an on- 
board earth resoxn?ces data processor launched during the 1980-1990 time frame. 
Since about five years are required to design, build, check out, and launch 
such a system, a 1980 system would be based on 1975 technology, and a 1990 
system would be based on 1985 technology. 

In order to determine the feasibility of on-board processing we must 
first define the on-board processor. This requires that we define both the 
technology available for use in the design and the computational requirements 
required of the processor* The computational requirements depend on the algo- 
rithms that the processor must implement which in turn, depend on the data 
products that must be extracted from the data to satis:^ the users. Consequently, 
in order to determine the feasibility of on-board data processors we must 
start with a study of projected user applications to define the data format 
(data throughput rate, number of spectral bands, etc*) and the information ex- 
traction algorithms the processor must implement. Based on these constraints 
and the constraints imposed by the available technology we can design some 
on-board processors and evaluate their feasibility. The study plan is sum- 
marized in Figure 0(1). 



p processor designs 



Figure 0 (.1) Study Plan 

A brief description of the content of each of the succeeding sections of 
this report follows. 

Section 1. In this section we present the results of a survey of earth 
resources user applications and their data requirements, earth resources multi- 
spectral scanner sensor technology, preprocessing algorithms for correcting 
the sensor outputs and for data bulk reduction, and a candidate data format to 
be used in subsequent sections. 

Section 2. This section contains the computational requirements required 
to implement the data analysis algorithms, a review of some computer archi- 
tectures and organizations, a design of some computer architectures capable 
of handling the algorithm computational requirements, and a discussion of the 
on-board processor environmental effects. 

Section 3. The ability of the on-board processors designed in Section 2 
to implement the algorithms described in Section 1 in real time for the re- 
quired throughput data rates depends on the components that will be available 
at the time of system design. The lead-time required for design, procure- 
ment, fabrication, checkout, and launch is about 5 years, so that 1980-1990 
launches will utilize 1975-1985 technology. Consequently, we require ac- 
curate component and system technology forecasts for the next 10 years. 

Section 4. This section identifies the pertinent performance parameters, 
isolates the independent and necessary parameters, and relates these parameters 
to the system requirements for each of the user requirements discussed in the 


10 







preceding sections. This will* allow us to determine the feasibility of on- 
board processing for each user type in the 19180-1990, time frame and to per- 
form a tradeoff anal^^sis to determine .the sensitivity of our results to each 
of the important system parameters . ‘ 

Section 5* This section contains an overview of the entire study re- 
ported in detail in the preceding sections. Significant results and conclu-' 
sions are discussed, and recommendations for future actions by NASA are pre- 
sented. 


11 



Table Jl(l) 


Liat of S3anbols 


Abbre- 

viation 

ATS 

ALU 

cm 

CPU 

CCD 

CMOS 

ERTS 

EIFOV 

EAROM 

GSFC 

HML 

HTLU 

HTC 

IR 

I/O 

IC 

I^L 

IMP 

JTL 

km 

LARS 


LSI 

max 

ML 

MTBF 

MSI 

MHz 

MOS 


Abbre-' 

Definition viation 


Advanced Technology 
Satellite 

Arithmetic Logic Unit 
Centimeter 

Central Processing Unit 
Charge Coupled Device 
Complimentary Metal- 
Oxide Semiconductor 
Earth Resources Techno- 
logy Satellite 
Effective Instantaneous 
Field of View 
Electrically Alterable 
Read Only Memory 
Goddard Space Flight 
Center 

Hardware Maximum Like- 
lihood 

Hardware Table Look-Up 
Hybrid Technology Com- 
puter 
Infrared 
Input /Output 
Integrated Circuit 
Integrated Injection 
Logic 

Interplanetary Monitor- 
ing Platform 
Josephson Tunneling 
Logic 
Kilometer 

Laboratory for Applica- 
tions of Remote 
Sensing 

Large Scale Integration 
Maximum 

Maximum Likelihood 
Meantime Between Failure 
Medium Scale Integration 
Megahertz 

Metal-Oxide Semiconduc- 
tor 


MOSFET 


m 

pm 

yp 

pPML 

yPTLU 

Ps 

mv 

mw 

rain 

MIMD 

MISD 

MSS 
ns ‘ 
NASA 


nm 

n-MOS 

OPAU 

OEM 

pixel 

pj 

PE 

PROM 

RAM 

ROM 

sec 

SPAU 

SOS 

SIMD 

SISD 

SSI 

TLU 

TTL 

SUMC 

V 

W 


Definition 


Metal-Oxide Semiconduc- 
tor Field Effect 
Transistor 
Meter 

Micrometer 
Microprocessor 
Microprocessor Maximum 
Likelihood 

Microprocessor Table Look-Up— 

Microsecond 

Millivolt 

Milliwatt 

Minimum 

Multiple Instruction Stream-Multiple 
Data Stream 

Multiple Instruction Stream- 
Single data Stream 
•Multispectral Scanner 
Nanosecond 

National Aeronautics and Space 
Administrat ion 
Nautical Mile 

n-Chahnel Metal Oxide Semiconductor 
Optimal Pipelined Arithmetic Unit 
Original Equipment Manufacturer 
Picture Element 
Picojoule 

Processing Element 
Programmable Read Only Memory 
Random Access Memory 
Read Only Memory 
Seconds 

Serial Pipelined Arithmetic Unit 
Silicon On Sapphire 
Single Instruction Stream-Multiple 
Data Stream 

Single Instruction Stream-Single 
Data Stream 

Small Scale Integration 
Table Look-Up 

Transistor Transistor Logic 
Ultra Reliable Modular Computer 
Volt 
Watt 
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1 EARTH RESOURCES ALGORITHMS AND DATA SETS 


In this section we present the results of a survey of earth resources 
user applications and their data requirements, earth resources multispec-* 
tral scanner sensor technology, preprocessing algorithms for correcting 
the sensor outputs and for data bulk reduction, and a candidate data for- 
mat to be used in si±>sequent sections* 

Section 1*1 contains the results of the user requirements survey and 
their projected data needs in the 1980-1990 decade* The survey is based 
on existing literature and on personal interviews with earth resources 
experimenters* A survey of existing algorithms for carrying out the user 
requirements was also conducted* The maximum likelihood, perceptron, 
table look-up and clustering algorithms were examined in detail. 

Section 1*2 deals with present-day and projected state-of-the-art 
technology relative to electro mechanical and solid-state scanners and 
their characteristics . 

Section 1,3 contains a discussion of preprocessing algorithms for 
radiometric, gain, and offset corrections* Preprocessing algorithms for 
reducing the data bulk passed to 'the on-board processor using data com- 
pression and redundancy removal techniques are surveyed and analyzed. 

In Section 1.4 a candidate data format is developed* This is used in 
later parts of the study as a baseline format for designing on-board 
computer architectures. 

1*1 USER REQUIREMENTS 

1*1.1 Applications Survey 

The purpose of this section is to provide a brief survey of applica- 
tion areas which are most likely to be affected by remote sensing and 
automatic data interpretation techniques* 

Two types of image-related information will be considered in the fol- 
lowing discussion: (1) spectral information, and (2) spatial information* 

Spectral information is that resulting from the intensity response of a 
scene in the spectral bands of a multispectral scanner. Spatial informa- 
tion is the relationship between features in a scene* For example, the 
automatic classification of image elements in a lake can be easily carried 
out using spectral information, while the shape of the lake is best 
determined using spatial Information* 
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The discussion in this section is of a general naxure. Specific data 
requirements for each of the application areas discussed are covered in 
Section 1.1.2. 

1.1. 1.1 Agriculture (A) 

Agricultural applications are receiving more attention by more inves- 
tigators than almost any other application of Earth observations from 
space. The most important application in this area is that related to 
agricultural food production. The utilization of remotely-sensed data for 
crop acreage and • production prediction will continue to play a central 
role in the design of sensors and processing hardware because of the social 
implications of this particular application. 

Data processing functions related to food production include such 
varied tasks as: 

Al. Agricultural census 

A2. Plant species identification 

A3. Plant stress (due to insects, drought, or moisture) - 
A4. Soil conservation practices 
A5. Crop yield estimates 

The most challenging single technical problem associated with multi- 
spectral data analysis in crop investigations appears to be in the area of 
developing adequate automated machine techniques (i.e. , the utilizaxion of 
sample data to design automatic recognition devices). This problem will 
be discussed in Section 1.4. 

1.1.^C2 Coastal-Zone Studies (C) 

The monitoring of the physical as well as the biological environment 
in the coastal-zone regions is of great importance in preserving the 
quality of life in these regions. Data processing procedures related to 
this application must be able to support the following functions: 

Cl. Mapping of shorelines 

C2. Mapping of shoals 

C3 . Wetlands inventory 

C4 . Bath 3 nnetry determination 

C5. Bottom topography studies 

C6, Mean high/low water line determination 

C7 . Pollution detection 
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The wide range of possible features of interest, in terms of their 
spatial and spectral properties, requires that the minimum-sized objects 
of interest, with their 'attendant scene contrasts, be stipulated to permit 
an effective specification of required resolution,' sampling rate, and 
quantization levels. Resolution and sample rate is largely determined by 
the minimum feature size and quantization level is determined by the multi- 
spectral analysis techniques used for identification and classification. 

1.1. 1.3 Forestry (F) 

This application is in many respects similar to the agricultural pro- 
blem. It is generally expected that a remote sensing system based pri- 
marily on multispectral analysis which will perform most agricultural sur- 
vey problems will also be capable of solving many of the forestry survey 
problems. Some of the most important applications of remote sensing in 
forestry are; 

FI . Forest-nonforest delineation 
F2. Forest typing 
F3. Detection of forest fires 
F4. Plant stress detection 

These applications span a wide spectrum of sensor requirements. For 
' example , forest-nonforest delineation can be accomplished v/ith fairly 
coarse resolution, while individual tree counts and classification would 
require resolution in the order of one meter or less. It is doubtful that 
satellites in the foreseeable future will possess the data requirement 
capabilities to solve the latter problem, 

1.1. 1.4 Geography (G) 

Some of the major applications of remote sensing to geography are: 

G1 . Land-use change 

G2. Earth resotirces location 

G3. Delineation of urban/rural areas 

G4. Detailed urban structure 

G5. Traditional map preparation 

Of these, the ones related to the inventory and classification of man's 
activities are receiving the most attention- Because most of the above 
application areas involvfe not only spatial and spectral signature investi- 
gations but also generic pattern recognition, the data processing 
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requirements tend to be very complex. Also, for political reasons, there 
will likely be a need for archival copies of portions of the imagery 
rather than summary data expressed in numerical form. These requirements 
limit the applicability of on-board data processing to these areas. 

1.1. 1.5 Geology (L) 

Some of the most important applications of remote sensing to geology 

are: 

LI. Structural Geology (faults, folds, lineaments) 

L2. Geomorphology (landform classification) 

L3. Lithologic mapping 
L4 . Geologic hazards 
L5 . Landslides 
L6 . Volcano studies 

The information required for most geological investigations will include 
both spatial and spectral information. However, in some cases only 
spatial information will be of use. Examples of this would be cases where 
the spectral information existing in the mineralogy of the viewed area is 
obscured by vegetation or snow. In other cases, spatial Information may 
be secondary as, for example, in the case of relatively featureless 
terrain . 

The potential exists for partially automated analysis of the complex 
interrelated spectral and spatial information which is o.f significance in 
geological surveys. For example, certain lineaments and various other 
patterns could, in principle, be detected automatically. Some limited 
spatial frequency analysis has been attempted for geological investiga- 
tions. Although this type of analysis has the potential for automation, 
it only represents a limited sector of all image interpretation which is 
of interest in geological studies. Automating the geological survey pro- 
blem will undoubtedly be much more challenging than automating other 
applications of remote sensing such as crop surveys because of the fre- 
quent importance of simultaneous processing of both spatial and spectral 
information. 

1.1.1. 6 Hydrology (H) 

Because of the importahce of water-resources data acquisition to many 
users and government agencies, present and anticipated Earth resources 
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missions will find significant application of remote sensing to this area. 
Some of the high priority hydrologic applications are: 

HI. Delineation of land-water boundaries 

H2. Delineation of hydrologically-related terrain hectares 
H3. Hydrod 3 mamics , including floods, reservoirs, and estuaries 
H4. Water quality evaluation 
H5. Snov? cover and run-off evaluation 

Black-and-white infrared photography and near-infrared scanner imagery 
have proven to be extremely effective in detecting water surrounded by 
land. In infrared and near-infrared imagery, water appears black because 
of its absorption of solar energy in these wavelengths. The high contrast 
with the surroundings makes it possible to easily detect bodies of water 
such as lakes, rivers, streams, and reservoirs. Coupled with the repeti- 
tive coverage offered by a spacecraft, such a remote sensing technique 
could be utilized to monitor changes in the bovindaries of surface water. 

In the infrared and the near-infrared bands, healthy green vegetation 
reflects strongly and can be readily distinguished from water. Sensors 
operating in these ranges can therefore be useful in detecting water/vege- 
tation interfaces. One particularly notable application of this feature is 
the delineation of wetlands areas. Wetlands are difficult or nearly impos- 
sible to map by conventional ground survey methods or by regular black- 
and-white or color photography. 

It is evident that considerable variation in data processing require- 
ments exist for the various hydrology applications . listed above . Some 
■can be expected to be accomplished by gross evaluation of the presence or 
absence of water. Such examples may relate to evaluation of regional 
water resources and various other large-scale water inventory applications. 
In these examples no firm requirement exists for precise, geometric loca- 
tion data or for high spatial- frequency data content. Hence, these gross 
evaluation applications appear to be prime candidates for on-board data 
processing. 

Other hydrology applications, however, may require somewhat precise 
location data as well as spatial or textural data. Among these applications 
are flood hazard evaluation and precision mapping of submerged land forms . 
These applications will generally require the use of spatial and spectral 
data. 
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1.1.1. 7 Meteorology (M) 

Two of the principal areas where remote sensing is expected to yield 
results of metereological importance are: 

Ml. Cloud cover siirvey 

M2. Prediction and assessment of natiiral disasters 

The information required for cloud cover survey is the three-dimen- 
sional structure of the cloud cover of the globe. At first glance, the 
scene radiance present in the features of interest ’ (clouds) seems to lend 
itself well to a fairly narrow djmamic range or nonlinear quantization, as 
most cloud cover images commonly published are relatively bright with 

respect to the underlying land or ocean surface. When the total range of 

4 

lighting conditions (from a few hundred to nearly 10 foot-lamberts ) and- . 
the apparent demand by some users for good radiometric resolution over the 
full range are considered, this prospect may diminish in importance. The 
relatively low spatial resolution required does offer good potential for 
optimization of resolution, quantization level, and sample rate over the 
very wide scan field. However, as the spatial resolution requirements are 
specified at the point of swath contiguity, existing sensors produce ex- 
cessive spatial resolution at the nadir point which might profitably be 
traded for radiometric accuracy. Some of the relatively simple applica- 
tions are appropriate for near-term automated data processing. One of 
these wouli be that of percent cloud cover evaluation. This can be accom- 
plished by multispectral analysis and area integration. 

Disasters resulting from hurricanes , tornadoes , flash floods , etc . , 
could be predicted more accurately and damage assessed more quickly by em- 
ploying remote-sensing technology. A new insight into the behavior of 
severe weather has been gained from the analysis of ATS satellite time- 
lapse photography, radar, and motion picture photography of tornado cloud 
tops taken from aircraft flying at about 15 km. A telltale '.'cloud turret" 
rises from the cirrus cloud shield, punctures the tropopause, collapses, 
and falls back into the smooth top of the cirrus shield. Then, another 
turret pops up. This tornado associated phenomena is shortlived, lasting 
8 to 20 minutes. 

These new scientific findings point the way to use remote -sensing 
technology for the observation of severe weather. An area scanner capa- 
bility to view the whole earth every 2P minutes needs to be added to a 
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geostationary satellite. This area scanner would be pointed toward an 
active weather area (approximately 800 x 1000 km) and view that specific 
area every 2 or 3 minutes. Targets of 200 to 400 meters should be dis- 
tinguishable. Such data would allow meteorologists to observe the turret 
cloud phenomena and refine his severe weather warnings both in time and 
location. 

However, pin-pointing the occurrence of severe weather more- accurately 
still cannot .prevent the disaster. Once it has occurred, other remote- 
sensing technology can be used for damage assessment. Side-looking air- 
craft radar (SLAR) is an example. The geostationary area scanner might be 
used for this purpose. With another set of optics having higher resolution, 
smaller areas might be viewed ‘for damage assessment after-. the clouds 
have cleared. In an area measuring approximately 80 x 100 km, those 
features approximating 80 to 100 meters in size could be imaged as a first 
gross assessment of the destruction. Damage assessment data can be used 
for many purposes: emergency routing of traffic (especially rescue . 

vehicles), evacuation steps, estimates of emergency housing needs, and 
location of emergency food stations. The application of sensor tech- 
nology to the disaster problem is a matter of timeliness in acquiring the 
observational data in usable form, 

1.1. 1.8 Global Oceanography (0) 

The needs of the community of ocean researchers for s 3 moptic environ- 
mental data are naturally divided into two categories: (1) the coastal 

environment and, (2) the global oceans. The former was discussed in 
Section 1.1. 1.2. Some of the most important applications of remote 
sensing to global oceanography are: 

01. Study of biological processes 

02. Sea-ice sxirveillance 

03. ’ Study of current patterns 

Assessment of features of biological significance in the more produc- 
tive waters of the global oceans (approximately 10 percent of the total 
global area) require resolution on the order of 1 km and remotely acquired 

signatures of a quality that will permit the determination of chlorophyll 

3 

to within a factor of two for concentrations of 0.2 mg/rn or greater. 
Observations of chlorophyll in the open oceans for concentrations as low 
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as 0,02 mg/m are important for global ecosystem analysis- Pollution detec- 
tion and environmental impact are particularly important in monitoring the 
natural and cultural stress induced in the coastal region. 

Space techniques have demonstrated an immediate sensor application for 
observation of the polar environment. Resolution as large as 10 km may be 
used to delineate major boundaries and ice movements. However , the most 
pressing problem is to define > monitor ^ and forecast the amount and loca- 
tion, of open water in polar regions, in particular, the Arctic. Resolution 
on the order of 30-100 meters is needed, with a preference for infrared 
over optical observations because of the nonavailability of solar illumina- 
tion during the winter months. 

Currents on a global scale can be mapped with relatively low resolu- 
tion. Spatial dimensions of 10 km in extent would be meaningful for ocean 
applications. However, higher resolution may be desired to permit sampling 
between cumulus clouds. In general, the use of ocean color to monitor 
currents and biological and ecological features requires high sun eleva- 
tion angles and a scan which looks away from the sunside of the spacecraft. 
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1.1.2 Data Requirements Survey 


The following sections examine the data requirements of the applica- 
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tion areas discussed in Section 1.1.1. All discussions are keyed to the 
identifying labels introduced in that section. For example, general agri- 
ciJ.tural applications will be referred to as A, while specific areas within 
agriculture will be identified by A followed by a number. 

1.1. 2.1 Number of Spectral Bands 

Figure 1.1. 2(1) summarizes the spectral band requirements of the eight 
general areas discussed in. Section 1.1.1. This figure vividly illustrates 
the differences in needs of these areas. Only in the thermal IR region 
(0.8 to 1,4 microns) does there seem to be a real consistency of needs. 

1.1. 2.2 Repeat Coverage 

Typical minimum rates of coverage for various resolution requirements 
are shown in Figure 1.1. 2(2). The implications of these rates on the vol- 
umes of data generated are discussed in Section 1.1. -2. 4. 

1.1. 2. 3 Spatial Resolution and Field of Coverage 

The unit of resolution used in this study is EIFOV (Effective Instan- 
taneous Field of View) which is defined as "the minimum linear dimension • 
on the surface (at nadir) at which user specified characteristics can be 
discovered." Field of coverage is defined as the swath .width for a nadir 
pointing sensor. Typical resolution and field of coverage requirements 
for .the application areas discussed in Section 1.1.1 are summarized in 
columns two and three of Table 1.1. 2(1), Section 1.1. 2. 4. 


1.1. 2. 4 Data Rates 


The data rate Dj^ (in bits/sec) for each of the application areas dis- 
cussed in Section 1.1.1 may be calculated using the following relation: 


where 


V 


N, 


bp 


swath width (m) 

resolution along a scan line (m) 
satellite ground track velocity (m/sec) 
resolution along scan path (m) 
number of channels 

number of bits per pixel (picture element) (bits) 
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Wavelength (microns) 


KEY 


lO 

to 


u 



Agriculture (A) 
Coastal-Zone Studies (C) 
Forestry (F) 

Geography (G) 

Geology (L) 

Hydrology (H) 

Met ereo logy (M) 

Global Oceanography (0) 


Figure 1.1. 2(1) Spectral Requirements Summary 


EfFOV (meters) 



1 Day 3 Days 10 Days 1 Mo. 3 Mos 6 Mos 1 Yr 


Repeat coverage (minimum) 

Figure 1.1. 2 (2), Tjrpical Rates of Coverage 


Table 1.1. 2(1) shows typical data rate ranges for the application areas 
discussed in Section 1.1.1. All calculations are based on the following 
figures : 

selected from Figures 1.1. 2(1) and 1.1.2(2)‘. 

= Rp - resolution given in Figures 1.1.2(1)-1. 1.2(2). 

V = 6500 meters/sec. 



selected from Figure 1.1. 2(1). 

= 6, which gives a range of 64 gray levels for each picture 

element . 


Table 1.1.2(11) indicates the number of classes for each application 


area. 
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Table 1.1. 2(1) 

Typical Data Rate Ranges 





Field of 


No 4 of 

Data Rates 


Resolution (nt) 

Coverage 

(km) 

Channels 

(M bits/sec) 

Application 

min-max 

min-max 


min-max 

min-max 

A1 

30-50 

185 


4-7 

11.5- 

56.1 

A2 

30-50 

185 


4-7 

11.5- 

56.1 

A3 

30-50 

185 


4-7 

11.5- 

56.1 

A4 

10-30 

50 


4-7 

8.7- 

137.0 

A5 

30-50 

185 


4-7 

11.5- 

56.1 

Cl 

30-50 

200' 


6-20 

18.7- 

173.0 

C2 

30-50 

200 


6-20 

18.7- 

173.0 

C3 

30-50 

200' 


6-20 

18.7- 

173.0 

C4 

50-100 

200 


6-20 

4.7- 

62.4 

C5 

50-100 

200 


6-20 

-4.7- 

62.4 

C6 

3-10 

40 


6-20 

93.6-3470.0 

C7 

30-300 

200 


6-20 

.5- 

173.0 

FI 

50-100 

• 185 


4-7 

2.9- 

20.2 

F2 

5-10 

15-30 


4-7 

23.4- 

328.0 

F3 

10-30 

185 


4-7 • 

32.1- 

505.0 

F4 

30-50 

185 


4-7 

11.5- 

56.1 

61 

30-50 

185 


4 

11.5- 

32.1 

62 

■ 30-50 

185 


4 

11.5- 

32.1 

63 

50-100 

185 


.4 

2.9- 

11.5 

64 

5-10 

15-30 


4 

23.4- 

187.0 

65 

5-10 

15-30 


4 

23.4- 

187.0 

LI 

50-80 

185 


1-5 

1.1- 

14.7 

L2' 

50-80 

185 


1-5' 

1.1- 

14.4 

L3 

50-80 

185 


1-5 

1.1- 

14.4 

L4 

50-80 

185 


1-5 

1.1- 

14.4 

L5 

10-30 

15 


1-5 

.7- 

29.3 

L6 - 

100-200 

185 


1-5 

.2- 

3.6 

HI 

40-60 

200 


1-3 

2,2- 

14.6 

H2 

30-50 

200 


1-3 

3.1- 

26.0 

H3 

10-30 

50 


1-3 

2.2- 

58.5 

H4 

30-70 

200 


1-3 

1.6- 

26.0 

H5 

50-80 

200 


1-3 

1.2- 

9.4 

Ml 

200-400 

800 


2 

.4- 

1,6 

M2 

200-400 

800 


2 

.4- 

1.6 

01 

1-lOkm 

400 


4-20 

0.0* 

0.3 

02 

30-100 

200 


4-20 

3.1- 

173.0 

03 

1-lOkm 

200 


4-20 

0.0** 0.2 


* 624 bits/sec 

*5V3i2 

bits/sec 
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Table 1.1.2(11) List of Applications 

A1 Agricultural Census 

A2 Plant Species Identification 

A3 Plant Stress (Due to Insects, Drought., or Moisture) 

A^■ Soil Conservation Practices 

A5 Crop Yield Estimates 

Cl Mapping of Shorelines 

C2 Mapping of Shoals 

C3 Wetlands Inventory 

C4 Bathymetry Determination 

C5 Bottom Topography Studies 

C6 Mean High/Low Water Line Determination 

C7 Pollution Detection 

FI Forest -Nonforest Delineation 

F2 Forest Typing 

F3 Detection of Forest Fires 

F4 Plant Stress Detection 

G1 Land-Use Change 

G2 Earth Resources Location 

G3 Delineation of Urban/Rural Areas 

G4 Detailed Urban Structure 

G5 Traditional Map Preparation 

LI Structural Geology (Faults, Folds, Lineaments) 

L2 Geoinorphology (Landform Classification) 

L3 Lithologic Mapping 

L4 Geologic - Hazards 

L5 Landslides 

L6 Volcano Studies 

HI Delineation of Land-Water Bovindaries 

H2 Delineation of Hydrologically-Related Terrain Hectares 

H3 Hydrodynamics, Including Floods, Reservoirs, and Estuaries 

H4 • Water Quality Evaluation 

H5 Snow Coyer and Run-Off Evaluation 

Ml Cloud Cover Survey 

M2 Prediction and Assessment of Natural Disasters 

01 . Study of Biological Processes 

02 Sea-Ice Surveillance 

03 Study of Current Patterns , 
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1.1.3 Algorithms Survey 


The response of a multispectral scanner at any sampling time may be 
arranged in the form of a measurement or pattern column vector 


X = 


(1.1. 3-1) 


r-J 

where Xj^ is the amplitude response of the kth channel in the system and n 
is the number of channels. 

Given M categories or classes of desired classification (such as M 
crop types , for example ) , one of the principal approaches in pattern recog- 
nition system design is to determine from representative (training) data M 
decision functions d^(x), d 2 (>c), . . .', with the property that if an 


observation x belongs to the ith category, then 

d. (x) > d. (x) for all j 9^ i (1.1. 3-2) 

i J 

Once the decision functions have been estimated, they are used as the 
basis for automatic data classification Cl3. Thus, given a sample x of 
unknown origin, the sample is assigned to the ith category if d^^(x) yields 
the largest value upon substitution of the observation into all functions. 
Ties are resolved arbitrarily. 

i 

The boundary between the ith and jth classes is given by values of x 
for which d^(>£) = d.(x). In other words, the equation of the boundary 
separating these two classes is given by 

d.(x) - d.( 2 i) = 0 (1.1. 3-3) 
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For M classes, there are M(M - l)/2 such boundaries. The first tw algo- 
rithms discussed below are based on the decision- function concept. 

1.1. 3.1 Maximum-Likelihood Algorithm 

A maximum- likelihood decision rule is one in v?hich the decision func- 
tions are of the form 

dj^(ji) = p(3c/Wj^) p(Wj^) k = 1, 2, ..., M ■ (1.1. 3-4) 

where denotes the kth category, p(x/aij^) is the multivariate probability 
density function of the samples belonging to category and p(^) is the 
probability of occurrence of 

The form of in (1.1. 3-4) is determined by the nature of 

p(x^/tOj^) and p((Oj^). When p(5</w^) is the normal density, it can be shotm [1] 
that Eq. (1.1. 3-4) may be expressed in the form 

d^(x) = £n p (w^) - i ^n|c^| - [(ji - - m^)3 (1.1. 3-5) 

where £n is the natural logarithm, ^ and m^ are the covariance matrix and 

mean vector of the samples of category w. , and |(1 | is the determinant of 

C, ♦ The prime ('") in Eq. (1.1. 3- 5) indicates transposition. 

■He 

The parameters of the decision function shown in Eq. (1.1. 3-5) are 

and m^. Once estimated > these components complete3.y specify the 
decision function of each category, 'studies with multivariate remotely- 
sensed data indicate that the normal density assumption is valid for numer- 
ous classification tasks [2]. 

1.1. 3. 2 Perceptron Algorithifi 

Decision functions based on the perceptron approach are of the form 




(1.1. 3-6) 


where ^ - (w^^. .... 


. . , w, -TJ .1 ^ is the parameter vector, and <()(x) = 

I. — 

((})j^(x), 

valued functions of the patterns jc. Note that Eq. (1.1. 3-6) can represent 
any nonlinear function of finite degree, depending on the choice of ^(x). 
For example, a linear decision function may be constructed by letting N = n 
(see Eq. (1.1. 3-1)), and Iii this case, Eq. (1,1. 3-6) becomes 


d^U) = + 


+ W, X + w, , , 

kn n k ,n+l 


(1.1. 3-7) 
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Fr’om'Eq. ( 1 . 1 . 3 - 3 D, the eqtiation of the boundary between the ith and jth 
pattern classes is given by 

d.(x) - d.Cx) = (w., - w.^) X + ... + (w. - w.^) X 

1 — :i — il ]2 1 in ]n n 

+ Cw. .1 - w. . T ) = 0 ( 1 . 1 . 3-8) 

i,n+l 3 ,n+l 

which is the equation of a hyperplane in n-dimensional space. More com- 
plex functions may be constructed by varying the degree of nonlinearity of 
the functions {(j)^(x)}. 

In applying Eq. (1.1. 3-6) to a classification task, the functions 
{(j)^(x)} are specified. The problem then becomes the estimation of the co- 
efficients for each class. The procedure normally followed is to use 
sample patterns from each class to determine these coefficients. One of 
the simplest ways to accomplish this task is by means of the following basic 
perceptron algorithm. 

Given M pattern classes, (o^^, ( 1)25 * ■ • j assign arbitrary initial 
values to the coefficient vectors w, , w„ , ...» w,, . Then , at the kth itera- 
tive step in the algorithm, if a pattern x(k) belongs to class ( 1 )^ and 

d^C^(k)] > j 5^ i (1.1. 3-9) 

where d^[ 5 ?(k)] = 


w.(k + 1) = w, (k) 
-] -I 


for k=l, 2, ...,M (1.1.3-10) 

On the other hand, if j£(J<) belongs to class and for some H 

d^[x(k)] < d^[x(k)] (1.1.3-11) 

then the following adjustments are made on the ith and £th coefficient 
vectors 


w.(k + 1 ) = w.(k) + c X (k) 

—X —1 — 

Wjj^(k + 1) = ■ c X (k) 


(1.1.3-12) 


The other coefficient vectors remain unchanged. The factor c is a posi- 
tive correction increment. 

Basically, what this algorithm does is to change the parameter vectors 
only when an error in classification occurs. The procedure is said to 
have converged when an entire iteration through all sample patterns pro- 
duces no errors. Several illustrations of this algorithm may be found in 
reference [Ij. 
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1.1. 3. 3 Table Look-Up 


The table look-up approach is based on prestoring in fast, random- 
access core memory the desired answer Ce.g. , crop-tjrpe) for all combina- 
tions of multispectral scanner outputs from selected channels [3-6], Spe- 
cifically, each set of measurements from a given point on the ground is 
interpreted as that address in core memory v/here the answer can be retrived. 
Substituting the simple retrieved operation for the lengthy computations 
required by the conventional approach offers two advantages: 

(1) The processing time is reduced by more than an order of magnitude. 

(2) The multispectral scanner data can be processed by computers 
having minimal sophistication, complexity, and cost. 

These two advantages may inake it possible to use an on-board computer to 
perform the classification function in flight. A general approach to 
table look-up is given below. 

j 

Two Dimensional Case 

All computer based systems for classifying MSS data operate by par- 
titioning the multidimensional measurement space into non-overlapping 
regions associated with each known class. Measurements which are spec- 
trally dissimilar to all of the known classes are regarded as belonging to 
the so-called threshold class. From Figin?e 1.1. 3(1) it is clear that a 
pixel* with the particular measurement vector x = should be as- 

signed to Class 1 provided the following equations are satisfied 

L^(l) 1 (1) (1.1.3-13) 

L^djXj^) 1 X2 < (1,J^^) (1.1.3-14) 

L^(l) and H^^d) are the minimum and maximum values Xj^ can have to be 
associated with Class 1. Similarly, the quantities I. 2 (l>x^) and H 2 (ljXj^) 
are the minimum and maximum values can have for the specific case x^ = 
x^ to be associated with Class 1. 

From this example, it is clear that the procedure defined by Figure 
1.1. 3(2) can be used to decide whether or not a pixel having measurement 
vector X = (x^,X 2 )"* should be assigned to class C. The values of the 


* The term pixel is used often in image processing to denote an image or 
picture element. 
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H2(3) 

4(3) 

Hjd) 

HjCX.Xj) 


Ljd) 

Lj(3) 4<1) Xj Hj^U) Hj(3)' Xj 

Fig.' 1.1. 3(1) Explanation of Classification Algorithm for Tvro Dimensions 

constants L^^CC), H^(C) and L^CCjX^) for all Lj^(C) £ £ Hj^(C) have to be 

prestored in the random -access core memory. 

The values L^(C) and Hj^(C) can be regarded as the contents of two- 
dimensional arrays. Each array xTOuld require 8-bit (assuming 0 < x. < 
255) bytes of the core storage, where is the number of classes. Simi- 
larly, it is possible, to regard 12 ( 0 , x^) and H 2 (C, x^) as two two-dimen- 
sional arrays, each requiring bytes of core storage where is 

given by Eq. (1.1.3-16) below. 

^ From Figure 1.1. 3(1) it can be seen that Nj^(C) is the length which 
results from projecting that measurement space region associated with 
class C 

Ni(C) = Hj^(C) - L^(C) + 1 (1.1. 3-1! . 

'^^Imax ^ CNj^(C)] 1 1 C < (1.1.3-16) 

onto the x axis. In Figure 1.1. 3(1), = 3 and N = N (3) > N (1) > 

* -L ^ JL 7 H 3 X ■ -L JL 

N^(2). 



30 


The core requirements can be minimized by storing only the number of 
values that can actually assume for a given case rather than the fixed 
number dictated by the worst case. This dynainiG core assignment is 

accomplished by storing and retrieving L^CCjXj^) and H^CCjX^) using as a 
pointer in the one-dimensional arrays L^CP^^) and H 2 (P^). The value of P^ 
is computed using Eq. (1.1.3-17) 

Pj^ = 0^(0 + Xj^ (1.1.3-17) 


where 0^(C) is a class-dependent offset given by Eq. (1.1.3-18). 
^ 1 - Lj^(l) f or C = 1 


0^(C) = / (1,1.3-18) 

1 - L^(C) t 2 N^(i) 2 < C < 

^ i=l 

Figure 1.1. 3(3) shows the core arrangement which results from applying Eqs. 
(1.1.3-17) and (1.1.3-18) for the case shown in Figure 1.1. 3(1). 

Once the bo^indary information is prestored in the core memory, the 
classification proceeds. A hypothesis C is formed concerning which class 
gave rise to the measurement vector x = (x^jX^)”". The initial hypothesis 
is that class assigned to the preceding pixel. If this hjrpothesis fails, 
classes are tested in descending order of a priori probability. The class 



PIXEL IS ASSIGNED 
TO CLASS C 


PIXEL HOT ASSIGNED 
TO CLASS C 


PIXEL HOT ASSIGNED 
TO CLASS C 


Fig. 1.1. 3(2) Method for Determining VHiether Pixel is from Class C 
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Figure 1.1. 3(3) Organization of the Prestored Tables of 

Boundary Information for the Case Shown in Figure 1,1. 3(1) 

hypothesis C points to and H^(C), the minimum and maximum values 

can have and still be inside the measurement space region associated with 
class C. If bj^(C) ^ ^ Hj^(C), the hj^pothesis is tested further; other- 

wise a new hypothesis is formed as explained previously. The class h 3 ^o- 
thesis C also points to 0^(C), a class-dependent offset vfhich is added to 
x^. The resulting value points to L 2 (C,Xj^) and H 2 (C,Xj^), the minimum 
and maximum values Xg can have, for the given value of x^, and still be 
inside the measurement space region associated with class C. If L 2 (C,x^) 
— ^2 — pixel is assigned to class C, otherwise a new hypo- 

thesis is formed. If the measurement vector jc = (x^yX^)"^ does not lie 
within the prestored boundaries of any of the known classes, it is as- 
signed to the threshold class. 

Four-Dimensional Case 

The measurement vector for each pixel consists of n ^ 4 measurements, 
and the objective is to assign the pixel to one of the known classes or 
else the threshold class. The first step is to form a class hypothesis C. 
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The initial • hypothesis is that the pixel is from the same class as the 
previous pixel. If this hypothesis fails, the classes ape searched in 
order of decreasing a priori probability. The class hypothesis C is used 
to retrieve the number of those four channels which at'© most effective in 
discriminating class C from all of the other classes. 

The class hypothesis points to Lj^(C), H^(C) and 0^(C). A test is made 
whether or not ^ H^^CC). If not, a new hypothesis is formed. If 

so, the value is added to to derive the pointer which points 

to the core memory address where L^CCjX^), H 2 (C,x^), and 0^(0, x^) are 
stored. The process is continued in the same way to determine whether 
L^CCjXi) 1 ^3 1 H^(C,Xj^,X 2 ) and Lj^(C,x^,X 2 , 

Xg) ^ Xj^ ^ Hj^CCjXj^jX^jXg) . If all these conditions are met the pixel is 
assigned to class C. The first time one of these conditions is not satis- 
fied, a new class hypothesis is formed. If none of the class hypotheses 
satis:^ all four conditions, the pixel is assigned to the threshold class. 

1.1. 3. 4 Clustering 

Clustering is a data analysis technique by which one attempts to 
determine the "natural" or "inherent" relationships in a set of observations 
or data points. It is sometimes referred to as unsupervised classifica- 
tion because the end product is generally a classification of each observa- 
tion into a "class" which has been established by the analysis procedure, 
based on the data, rather than by the person interested in the analysis. 

Presently, the typical multispectral classification experiment is 
conducted as follows [7-163. The data, collected in a single region under 
favorable conditions by airborne or spaceborne sensors, are examined in 
their entirety by the experimenter, who decides which areas are most repre- 
sentative of the region as a whole. The samples from these areas are as- 
sembled to form a training set, which is characterized by ground truth in- 
formation delineating the terrain types of interest. A statistical cate- 
gorizer, or decision box, is constructed on the basis of the statistical 
parameters extracted from the training set. The classification perfor- 
mance is then evaluated on another portion of the data (the test set) for 
which the location and extent of the different types of ground cover are 
also known to the experimenter. 


33 



The details of the various experiments differ vfith respect to the 
sources of data, the method of labeling the training ^nd test seta, the 
terrain types to be identified , the number of spectral bands , the degree of 
the statistical sophistication of the categorizer, and the methods of 
evaluating the results, but the general scheme of classifying samples of 
the test set according to their similarity to a preselected training set is 
the same. 

Extrapolation of the performance levels evaluated in this manner to 
an operational satellite data gathering system is doubtful. In most experi- 
ments both the training sets and the test sets are selected on the basis 
of visual inspection of xhe available data. Terrain classes without large 
uniform representation are frequently deleted from consideration, as are 
regions of abnormally high variability. The data are divided into train- 
ing and test sets in such a way as to enhance the probability of success- 
ful classification. Even in as well conceived an experiment to extend 
the spatial recognition range as that described in [12]], an intruding 
cloud required complete reassignment of the intended training region. In 
view of the amount of information to be collected by the satellite systems, 
it seems unlikely that a considerable fraction of this data can be visually 
screened in time to allow modification of the required decision parameters. 

If, on the other hand, interactive systems are developed to a sufficient 

^ J 

4egree to allow human analysis of much of the imagery, then the whole 
concept of automatic terrain classification becomes superfluous. 

The point of departure from standard statistical classification tech- 
niques can be in the application of an unsupervised learning approach, by 
means of clustering algorithms [173, to circumvent the difficulties of 
collecting representative training sets. 

To get an intuitive idea of what is meant by natural or inherent 
relationships in a set of data consider the example shown in Figure 1.1. 3(4). 
If the spectral reflectance of vegetation in a visible wave band were 
plotted against reflectance in an infrared v/ave band, dry vegetation and 
green vegetation could be expected to form discernible clusters. 

If the data of interest never involved more than two attributes 
Cmeasurements or dimensions), cluster analysis might alvrays be performed 
by visual evaluation of two-dimensional plots such as that in Figure 1.1. 3(4). 
But beyond two or possibly three dimensions, visual analysis is impossible. 
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X2 ('infrared) 



Figure 1.1. 3(4) Example of data clusters 

For such cases, as in multispectral data, it is desirable to develop com- 
puter algorithms to perform cluster analysis. Haralick and Kelly [18] have 
presented two such algorithms : the first partitions the image sequence and 

the second partitions the measurement space. In both, the partition is 
constructed by finding appropriate center sets and chaining to them all 
similar points. The resulting clusters are simply connected. 

The reader interested in the many possible ways of defining clustering 
in quantitative terms may consult references [19] and [20]. Essentially, 
the definition of a clustering algorithm depends on the specification of 
two distance measures: a measure of distance between data points or indi- 

vidual observations j and a measure of distance between groups of observa- 
tions (clusters). Figure 1.1. 3(5) is a block diagram for a typical cluster- 
ing algorithm [21]. The point-to-point distance measure is used in the 
step labeled "Assign each vector to nearest cluster center" . The distance 
between groups of points (clusters) is calculated in the step "compute 
separability information". 
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1.2 MSS SENSOR TECHNOLOGY 
1.2.1 Electromechanical Scanners 

Electromechanical scanners are capable of producing geometrically- 
registered, multispectral images in spectral channels extending from the 
visible to the thermal infrared. In addition, accurate radiometric measure- 
ments can be made. Image formation is now accomplished by mechanical 
scanning. Typically, mechanical motion causes the scene to be sampled in 
the cross-track direction by a detector or array of detectors while satel- 
lite motion provides the orthogonal scan component. 

The information extracted from scanner data is in the spatial, spec- 
tral, and temporal distribution of radiation from a scene. For the most 
part, sensor advancement means improving the spatial resolution for a given 
operating distance. More recently, attention has been given to spectral 
distribution and automatic classification based on the spectral information 
from the scene. Since the spectra of the vegetation features vary with 
their growth cycle or season, the spectral classification task becomes 
easier when temporal variations are included. Thus, the advanced images 
may be viewed as a high resolution multispectral scanner with the ability 
to observe a scene periodically. 

When photomultiplier detectors are employed, the internal gain is 
assumed to be sufficienrly high that amplifier noise can be neglected in 
comparison with photoelectron (shot) noise, in this case, the limiting 
noise can be expressed as the square root of the number of photoelectrons 
released' during an integration period from a resolved elemental area of 
background while the signal is the difference between the corresponding 
numbers’ of photoelectrons associated with target and background. ■ ■ 
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Photomultiplier tubes are . currently used. in multispectral scanner 
systems for the short wayelength channels . For high resolution systems , 
PMT’s are -used for wavelengths shorter than about 800- nm; in lower reso- 
lution systems, PMT's may not be used at all or only at wavelengths shorter 
than 500 nm. 

Silicon planar p-n junction photodetectors represent the most advanced 
semiconductor technology available at this time and provide good resolution 
and parameter control. These units show good high-frequency response, high 
quantum efficiency, low noise, and require no special cooling for most 
applications. Because the photodiode has no internal signal gain, high 
performance channels are usually limited by the noise of the first amplifier 
stage. Silicon photodiodes are normally used in multispectral scanning 
systems for xravelengths between about 0.4 ym and 1.1 ym. Other semicon- 
ductor materials are used in the fabrication of p-n junction photodetectors 
for wavelengths longer than 1.1 ym. In general, infrared photodiodes must 
be operated at low temperature to achieve optimum performance , vrith the 
longer vravelength detectors requiring the lower operating temperature'. 

Photoconductive detectors are most useful at infrared wavelengths 
longer than 3 ym. These devices are essentially variable resistors in 
which the conductivity of the bulk material increases monotonically with 
the magnitude of optical power absorbed in the active volume of the units. 

1.2. 1.1 Detector Cooling Systems 

Future infrared imagers and scanners will utilize cooled quantum- or 
photodetectors. These detectors require cooling to cryogenic temperatures, 
i.e., below approximately 120°K to achieve background-limited performance. 

In general, the long wavelength cutoff and detectivity, as well as 
other detector parameters, are determined by the operating temperature. 
Photodetectors operating in the 8-13 ym atmospheric window require lower 
operating temperatures than those operating in the 3-5 ym region. For • 
remote sensing applications, future imagers and scanners , particularly 
those operating in the 8-13 'ym region, will utilize intrinsic photodetec- 
tors cooled to lOO^K or below to achieve background-limited performance. 

_ Within the next . decade it is anticipated that detector cooling require- 
ment for airborne and space-borne infrared systems will generally lie in 
the 50-120°K region with perhaps a few applications requiring temperatures 
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as lov 7 as 20^K. The cooling capacity at these temperatures \ 7 ill range from 
a few milliwatts for a single photovoltaic detector to perhaps a v/att* for 
large arrays of photocdnductive detectors.* 

The basic design parameters for a cryogenic detector cooling system 
are: the required operating temperature and temperature stability ^ heat 

load at the operating temperature ^ the alignment requirements of the cooled 
detectors relative to the optics 5 and the reliability and operating life of 
the system. 

For spaceborne systems , where weight and power are usually limited , 
the selection of a cooling system involves a tradeoff between the detector 
and cooling system parameters for a given mission duration. 

Three basic types of cooling systems are considered for spaceborne 
systems . 

(a) Passive radiators for spaceborne system vihich cool detectors by 
direct radiation to the low- temperature sink of deep space. 

(b) Open-cycle systems which use fluid or solid cryogens stored in a 
dewarj or stored high pressure gas which provides refrigeration by the 
Joule- Thompson effect. Solid cryogens are only applicable to spaceborne 
syLs terns. 

(c) Closed-cycle systems employing a mechanical refrigerator using 
helium gas as the working fluid or closed-cycle Joule- Thompson systems. 

For airborne application ^ open cycle systems using liquid cryogens 
stored in a dewar, or closed-cycle refrigerators are the only logical sys- 
tem choices. Open-cycle systems using liquid nitrogen or liquid helium 
have been used to cover the temperature range from 42*^to 70^K, These sys- 
tems are relatively simple, low in cost, have good temperature stability 
and do not introduce mechanical vibrations or microphonics to the focal 
plane. They are capable of providing continuous refrigeration for the 
duration of a single aircraft flight without resupply. 

For spaceborne applications, passive radiators are applicable for 
small heat loads down to their lowest temperature limit. Solid-cryogen 
coolers are also applicable for low heat loads, and can be- used to reach 
temperatures below those achievable with radiators. For a specific space- 
borne instrument, the particular constraint of the instrument should be 
considered in detail before selecting the cryogenic cooler. 
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Many of the cooling requirements ; for spaceborne applications can 

probably • be . met with passive. radiators. Current passive radiator systems 

are not large enough to handle the increased heat loads anticipated larger 

detector arrays. Large passive radiators having a cooling capacity of 

o 

approximately 1 watt at 70-80 K temperature level are quite sufficient for 
many spaceborne applications. 

Open-cycle cooling systems may not be applicable for the one-to-two 
year missions with focal plane heat loads which are expected to exceed 100 
milliwatts. 

1.2. 1.2 Optical Systems 

Multispectral electromechanical scanning systems can be classified 
broadly as image-space and object- space scanners. The former requires an 
imaging system that covers the full field of view with the required resolu- 
tion. Object-plane scanners take the load off the optical system and place 
it on the scanning system sequencing a narrow optical field of view. 

Narrow- field systems usually include simple spheres, parabolas, 
Cassegrain, Newtonian, Ritchey-Chretien, and Dall-Kirkham systems. Using 
the most expensive optical system (the Ritchey-Chretien), one can obtain 
resolution of 0.02 mr for about 1 degree full field (17 mr). 

Wide- field systems (10 degrees or more) are much more difficult to 
design and build. They have scanning mechanisms that are simpler to imple- 
ment, hox-iever, and are therefore of considerable interest. The two, main 
candidates -for such wide-angle optical systems are the Schmidt and the 
Bonx'fers-Maksutov system.- Both provide resolutions on the order of 0.1 mr 
out to angles of about 15 degrees. 

1. 2.1.3 Hadamard System 

A newer optical technique that makes use of the multiplex advantages 
in _ spectroscopy is an image-space scanner called the Hadamard system. 

Where the normal spectrometer uses' one exit slit and rotates a prism for 
multiple imaging, the Hadamard system uses a mask that provides multiple 
exit slits. A system of this type requires that the entire system be 
imaged. Any variations in the radiation from the opaque elements causes 
variations in the apparent scene radiation. The nonuniformities cause 
scene noise and uniform rad-iation can increase photon noise. This system 
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is of great . advantage for sampling -a acene with a few variarions, compared 
to. a single -.detector initially. Scanning -.noises', however , .i^. a serious pro- 
blem. 

A comparison can be made with a system' employing a linear array of 

detectors. For an array of n detectors to cover the n resolution elements 

1/2 

of the scan line, the gain in signal'-to-noise ratio is n . The detectors — 
probably cost more, having a varying responsivity from element to element 
and require a bias supply, cold shielding, preamplifiers, etc. The Hada- 
mard scanner uses one detector which reduces cost. The number of pre- 
amps and bias supplies, weight, and cooling requirements are lower, but 
the system requires a scanning mask, sampling system, and data processor 
for inverting the matrices. 

1.2.2 Solid-State Scanners 

Solid state scanners are in the early stages of development at the 
present time, but they offer important advantages for future earth-orbit 
satellite missions. Advantages include: no mechanical scanning, built-in 

geometric accuracy, the high quantum efficiency of silicon for the visible 
range, high resolution, good stability, low voltage operation, and improved 
signal- to-noise ratio performance. 

The fundamental performance criteria for integrated, self-scanned 
solid-state arrays closely parallel those for alternative sensors. .Among, 
the important criteria are: 

CD Geometric resolution which is determined by element-size spacing 
and number. 

C2) Quantum efficiency and its variation with length, as determined by 
reflection and absorption losses in layers overlying the silicon. 

C3) Scan rates using integration and, together with element-size, 
quantum efficiency, and available light, determine available signal. 

These factors set limits on d 3 raamic range and boundary conditions on 
signal-to-noise ratio of the sensor itself. 

Of importance in practical -systems is device stability, particularly in 
short term and, although of lesser consequence, power' supply requirements. 

Scan rate, taken with the number of elements. Imposes bandwidth require- 
ments per monolithic chip and for the total system. Dynamic range require- 
ments and signal compression and encoding requirements are dicta-ced by 
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system constraints . on the- one handj on the other,. they are limited by 
available transmission bandijidth. and. by precpmpression ^nd precoding noise • 
associated with. the' signal i 

As with other sensor s ,• solid state self-scanned arrays require geo- 
metric, radiometric, and electrical calibration. ‘Intrinsically, geometric 
calibration is straightforward and, once established, is. subject only to 
mechanical distortions of the sensor structure. Element-to-element varia- 
tions in photosensitivity or dark current must; as with other sensors, 
be -calibrated out.' Depending on the particular sensor design, additional 
calibration must be provided when additional per-element variations exist; 
for example, photo'transistors require gain or linearity calibration, usually 
at two or more points, plus additional calibration for temperature varia- 
tions. In general, determination of photosensor geometry by the photo- 
masking process during fabrication and the chemical stability of the sili- 
con/silicon dioxide system make geometric, radiometric and electrical 
characteristics of solid-state self -scanned arrays very stable. 

Silicon-detector arrays represent an Imaging technique that makes use 
of linear arrays of solid-state detectors operating in what is termed a 
Vpushbroom scan" .mode. In such a system, a detector array is used to 
image the scene in the cross-track direction and spacecraft motion is- used 
.to provide the orthogonal scan component. The primary- advantages resulting 
from the use of these arrays are: 

(1) With the array oriented in a cross- track configuration, continuous 
coverage of a wide swath width of terrain can be obtained. Mechanical 
scanning is eliminated, since the satellite subpoint motion along the 
ground track provides the single-axis scanning motion that is required, and 
the detector elements are interrogated electronically. 

(2) Large arrays can be formed containing several thousand detector 
elements. The use of high-density arrays offers high resolution capability. 

C3) The precise geometric alignment of the detector elements resulting 
.from the use of micro circuit fabrication techniques, in addition to the 
precision of alignment within the optical system, offers an advantage in 
accuracy of ground reconstruction of 'images, 

-C4) The detectors are operated in an integration mode, thus providing 
high scan efficiency. The exposure time of each detector element is 
limited only by the permissible image motion at the optical focal plane. 
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The. high scanning efficiency -results in, a. long exposure time, .increasing 
the' signal-to^noise ratio and.resoiution. > 

There are three. candidate sensors:' Photodiodes, .phot otransistors and 
charge-coupled devices (CCD).' 'All three are -similar in mechanism, limita- 
tions i and po-tential of the sensing-.portions, but differ significantly 'in 
the arrangements for signal amplification and scanning. In photodiode and 
pHototransist or • arrays , scanning is 'accomplished by' switching, while, in 
charge-coupled devices, scanning is accomplished by movement of potential 
wells . 

In general, there are two self-scanning modes: digital multiplexing 

and analog charge-transfer. The photodiodes and phototransistors ordinarily 
use digital multiplexing. The CCD sensors have analog charge-transfer 
readout. Since' all of these sensors are silicon, their operation is re- 
stricted to the visible and near infrared' (1.1 ym). The readout circuitry 
is generally on the same silicon chip as the array of detector elements 
(usually 100 or more sensors per chip). - 

At present, silicon photodetector elements’ can be assembled into line 
arrays of several-thousands elements with a center-to-center element spac- 
ing of 15 ym for a, high-resolution, push-broom mode- of scanning. 

Charge- coupled devices are in the early developmental stage, but CCD 
scanners possess better performance capabilities for multispectral earth 

imaging from orbit. Using the' basic CCD principle,, it is possible 'to 

transport a photon-generated signal charge over long distances within a 
.chip and to sense the charge- with a preamplifier having potentially an 
extremely small input capacitance. This in turn leads to amplification 
with very low input noise levels . Using the buried-channel CCD concept , 
it is possible to achieve very- high transfer- efficiencies over wide d 5 mamic 
ranges, which minimizes image distortion. Using.. transparent gates over 
the photosensing regions, it is possible to achieve high net quantum effi- 
ciencies . 

CCD is a new class of semiconductor sti’uctures normally operating in 
Cthermal) non-equilibrium and utilizing, as the signal carriers, minority 
charge transported by moving- potential wells. In essence, therefore, a 
.CCD is a nearly ideal semiconductor analog shift -register . The CCD con- 
cept permits the design of highly complex functional devices at potentially 
low cost. In addition CCD has the attributes of: 
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-(1) silicon fabrication simplicity 

(2) high packing density 

(3) high reliability 

(4) low power requirements 

(5) potential low-noise analog signal processing 

Self-scanned arrays also appear worth exploring for infrared applica- 
tions » although the advantages in the infrared range appear not to be as 
significant as in the visible spectrum. The limitations result from the 
fact that the devices are still in the early developmental stage and infra- 
red detector development has not • reached the technological maturity of 
silicon detectors and arrays. Major differences are contained in the fol- 
lowing: (a) present HgCdTe detectors operate generally in the photoconduc- 

tive mode; (b) small scale (approximately 100-element) IR arrays have been 
manufactured only through the physical assembling of individual detectors; 
and (c) 1/f noise in the infrared detectors is considerable higher than in 
silicon detectors. The implications of these differences for the IR systems 
are as follows; (a) present systems require a drop in the incoming radia- 
tion; (b) detector bias (1 mW per detector) is required; and (c) the low 
impedance of the detectors (25-2505^) requires high-power preamplifier 
stages. 

It is c\arrently possible to make high-density arrays of IR detectors 
(photo conductive, photovoltaic) out of narrow band gap semiconductors for 
use in mechanical scanners. However, it is not possible yet to inte- 
grate these detectors with, their associate electronics onto one chip; this is 
because the technology for fabricating transistors in materials such as 
InAs, InSb, HgPdTe, PbSnTe, etc., has not been established. It is 
possible to consider hybrid devices consisting of IR detector arrays ce- 
mented to a silicon -chip that contains the necessary readout circuitry. 
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1.3 PREPROCESSING ALGORITHMS 


1.3.1 Sensor Corrections 

In order to illustrate the kind of MSS sensor corrections that will be 
required on-board a satellite we briefly review the ERTS system which is 
typical of most multispectral scanners. 

The purpose of the MSS system is to accurately produce images of the 
earth from a low altitude, satellite. The energy received from the earth’s 
surface, both reflected and emitted, is sensed and digitized by the system. 
The digitized data are transmitted to the ground where they are recorded 

% 

for later processing* The recorded data are either analyzed directly by a 
computer or transferred first to a photographic transparency* 

The recorded data on the magnetic tapes includes timing and radio-- 
•metric calibration information that can be converted by computer programs 
into spectral signatiires for direct analysis* In photographic transparent- 
cies, the density of the film is proportional to the scene energy levels. 
There are twenty- four channels in the system and six channels are involved 
in making a one band picture . ‘ To produce a perfect pictiire of the scene 
imaged, the channel gain and the offset of each of the six channels in a * 
band must be known exactly. This is necessary so that the correct film 
density can be assigned to the digital words generated by the multiplexer. 
If this condition is met, individual lines in the picture at almost all 
density levels will not be discernible. However, if the gains and offsets 
are in error in one or more channels, the individual adjacent lines which 
should have nearly the same scene Information will have different average 
levels of density in each channel and, consequently, will be displayed as 
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stripes in the picture. Also, for signature analysis the channel gain and 
offset are equally important to maintain relative response between the four 
bands. For these reasons, a calibration system is required in the MSS 
system to report periodically on the status of the channel gain and offset. 

1. 3.1.1 Radiometric Response 

The integrating sphere is used as a light source to establish absolute 
gain/sensitivity for the scanner. The integrating sphere is calibrated 
using a primary standard source. Table 1. 3.1(1) gives a listing of the 
radiance required for each channel to produce full scale signals at 
the output of the multiplexer. The integrating sphere used to determine 
the radiance required was calibrated at 6SFC. The details are given in the 
Acceptance Test Report (HS 324-5196). 

Table 1.3. 1(1) Radiance Necessary to Produce Full Scale at the 

2 

Multiplexer Output (mw/cm • ster) 


CHANNEL 


BAND 



.1 

2 

3 

4 

A 

2’. 45 • 

1.86 

1,84 

4.81 

B 

2.43 

1.95 

1.68 

4.80 

C 

. 2.34 

1.34 

1.73 

4.76 

D 

2.33 

1.87 

1.80 

4.84 

E 

2.39 

1.92 

1.74 

4.70 

F 

2.35 

1.96 

' 1.78 

4.70 


1.3.1. 2 Determination of Gain and Offset 

Over a wide range of shades of gray, the eye responds to ratios, in- 
stead of absolute levels. The eye can detect a sharp-edged junction be- 
tween two fields that differ by little more than 2 percent in level. As a 
result the gain determination procedure should be designed to compute the 
gains of all channels in a band to within an uncertainty of 2 percent peak- 
to-peak. Furthermore, based on visual impressions of a test transparency, 
it has been determined that an acceptable offset spread in a band for strip- 
free pictures is 30 mv. The gain and offset are defined in Figure 1.3. 1(1). 
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Figure 1.3. 1(1) Determination of Gain and Offset 

The gain and offset of each, channel in a band are determined by pro- 
cessing the gray wedge in conjunction with, the sun calibration pulse.. The 
sun- pulse is used to modify the gain and offset computed from the gray 
wedge. 

The reference (initial) gain and wedge have to be known a priori be- 
cause the calibration gray wedge .(different for each scanner channel) can 
only provide information to compute the change in gain and offset from a 
previously recorded reference gray wedge. The integrating sphere is all 
that is needed to obtain the reference calibration data. Two different 
radiance outputs are used to establish the initial offset. One of these 
two settings, or an average of several settings, is used to establish the 
gain. 

1,3.2 Data Bulk Reduction 

The extremely large volume of data generated by a MSS imposes a 
severe computational burden on the on-board processor. The application of 
appropriate data compression and/or feature selection techniques can some- 
times reduce the severity of this problem. The algorithms usbd should not 
destroy more than. the maximum acceptable amount of Information and should 
be capable of efficient compression of different kinds of data (i.,e. , 
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vegetation, desert, mountain, etc.) that a multispectrhl sensor might 
encounter over changing terrain. 

1.3. 2.1 Coding 

Data compression algorithms for multispectral scanner data, (e.g., 

ERTS data), have been developed, analyzed and tested to the point where 
these algorithms and their performances are reasonably well understood. 

From the results of these studies and experiments there appear to be two 
basic types of data compression algorithms that are applicable: (1) trans- 

form coding; (2) coding by BLOB. 

1.3. 2. 1.1. Transform Coding 

Transform coders perform a sequence of two operations. The first opera- 
tion is a linear transformation that transforms the set -of statistical 
dependent data elements into a set of more independent coefficients. The 
second operation is to individually quantize and code each of the coeffi- 
cients. A variety of linear transformations is available to implement 
transform coding or MSS data, but the eigenvector transformation is the 
optimum linear transformation in two senses: The mean square error between 

the original and reconstructed data is less than for any other linear trans- 
formation; also, it eliminates all correlations in the data. 

In this method an N-vector of data samples 

x= (Xj^,X2, ..., Xj^)^ , (-1..3.-2-1) 

is transformed into an n-vector of coefficients 

y = CyjLs y2> (1.3. 2-2) 

by the transformation 

y = T (x - m) (1.3. 2-3) ‘ 

where m is the mean vector of » 

m = E Cx D (1.3. 2-4) 

and ^ is a n x N matrix whose rows are the eigenvectors of the covariance 
matrix of 3^: That is, if 

^ = E C(j£ - m) (^i ~ (1.3. 2-5) 

Then the rows ^ of X ^ solutions to the equation 

Ct=Xt (1.3. 2-6) 

corresponding to the n largest eigenvalues A = > • • > X^. A replica 

^ of the data is reconstructed from the coefficient y by transformation 
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X = T" (y f £) (1.3. 2-7) 

For n = N the reconstructed data x is identical to the original data 
i.e., ^ = X. For n < N some error (x - x^) is incurred. However , the 
eigenvector transformation results in the least mean square error 

= {E I I X - xi 1^} ' ■ (1.3. 2-8) 

of ail linear transformations. The elements Y 2 * ♦ • • s of ^ 

correlated and have variances given by the eigenvalues , . . . , X^ . It 

can also be shown that xhe mean square error Eq. (1.3. 2-8) is given by 
„ N - n N 

e = 2 X. - Z X. = 2 X. . (1.3. 2-9) 

1=1 1=1 i=n+l 

so that retaining only the first n of the N coefficients results in a mean 
square error given by the sum of the variances (eigenvalues) of the dis- 
carded coefficients. 

4. 

In any particular application the number of coefficients *n that must 
be retained depends on how fast the eigenvalues > . .* . > X^ 

decrease. If the data samples are not correlated, all N eigenvalues have 
significant values and we must choose n - N for negligible distortion. On 
the other hand if the data samples are highly correlated thePx the eigen- 
values decrease rapidly so that only the first few have significant values 
and all but the first few can be discarded. In this case n « N and signi- 
ficant sample compression ratios 

^ (1.3.2-10) 

s n 

can be obtained with negligible distortion. Hence the compression ratio 
that can be achieved depends on the amount of redundancy in the data and 
the amount of error that can be tolerated. 


For correlated data the n retained coefficients have different rms 
values > . . . > so that a different number of bits should be 

used to code each coefficient. The mean square quantization error is mini- 
mized by choosing m^ 'v log X^. This is called block quantization. If the’ 
original N data samples have m bits each, then the bit compression ratio 


achieved by the eigenvector transformation and block quantization is 

^ ~ (1.3.2-11) 

m, +m^+...+m 
12 n 

2 2 2 
It can be shown that R (e ) .can be made less, than R (e ) (for some e ). 

b s 
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The eigenvector transformation and the block quantization method have 
been used to encode waveforms [2,33 and pictures [4,5J. Monochromatic 
pictures were encoded and reconstructed with no -noticeable distortion at 
2.5 bits per picture element. Ready and Wintz [6] implemented the same 
technique for MSS data. They considered two additional transforms (Fourier 
and Hadamard) also applicable to the MSS data. Both of these transforma- 
tions are non- source-dependent in that the set of orthonormal basis vectors 
is fixed regardless of the source characteristics. The rows of the Fourier 
transformation matrix are the sampled, harmonically-related sine and cosine 
functions. The rows of the Hadamard matrix are the discrete version of ortho- 
gonal Walsh functions. 

1.3. 2.1.2 Coding by BLOB 

BLOB first subdivides the multi-imagery into 2x2 sxabsets (pixel 
groups) of picture elements (pixels). A hypothesis testing algorithm com- 
putes the mean (gray level) and variance (texture) of each 2x2 array and 
compares the first and second order statistics of adjacent subsets. Adja- 
cent subsets having similar first and second order statistics are merged 
into blobs. In this manner the entire imagery is partitioned into blobs 
such that the picture elements within each blob have similar gray levels 
and/or variances. By varying either of two parameters, the amount of con- 
sideration given to gray level and texture can be adjusted [8]. 

The boundaries may be weak in some spectral channels and strong in 
others. To take this consideration into account, the variance (F-test) and 
mean (student t-test) tests are implemented separately in each channel. 

This has been called a ’’multiple-univariate approach." 

Finally, BLOB transforms the boundaries into a sequence of binary 
digits with the help of contour tracing algorithm discussed in [33. The 
output of the contour tracer is in the form of directionals which are coded. 

1.3. 2. 2 Feature Selection 

Use. of coding procedures reduces the data volume that must 
be processed by the data analysis algorithm. This in turn reduces both the 
time required to read the data into the computer and the time required to 
do the mathematical operations on the data. For example, classification 
algorithms are usually designed to operate on fewer number of spectral 
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channels because including more channels requires excessive computer time. 

In that case it is desirable to know the relative importance of the indivi- 
dual features from the classification viewpoint. This suggests the study 
of feature selection; that is, the selection of subsets of feature measure- 
ments from the complete set. 

The eigenvector transformation discussed in Section 1.3. 2. 1.1 is also 
one of the most popular techniques for feature selection. The procedure is 
usually referred to as the principal components method [1] and consists of 
the following steps. 

(1) Compute the covariance matrix C of the given data. This matrix is 
n X n. 

(2) Obtain the n eigenvectors and associated n eigenvalues of The 
eigenvectors are n-dimensional. 

(3) Choose the m eigenvectors associated with the m largest eigenvalues 
of where m < n. 

(4) Form an m x n transformation matrix whose rows are the m eigen- 
vectors selected in Step (3). 

(5) Reduce all original vectors into a set of y vectors by means of 
the transformation 

y = A”^x 

The resulting vectors are of lower dimensionality but, according to Eq. 

1. 3.2-8, this reduction results in the least mean square error. 
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1.4 CANDIDATE FOR DETAILED STUDY 


In this section a detailed study of a particular application area is 
carried out. The area which has been selected for this study is the remote 
detection and classification of agricultural crops. This choice is based 
on the importance and relevance of this problem in terms of human as well as 
economic considerations. 

There are two principal items of information relative to agriculture 
remote sensing systems - (1) total crop acreage and (2) total expected 
yield. On a world-wide scale, the first item requires classification of 
ground crop information into one of approximately twelve different crops 
(pattern classes). Once the total acreage of each crop has been determined. 
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the yield can be calculated by assessing the state of health’ of each crop . 

It has been - suggested that a minimum of three states of health ranging from 
normal to abnormal should be adequate for a close approximation of the 
expected yield from a particular crop El]. Therefore, the crop classifica- 
tion problem maj be divided into two stages. The first stage is used for 
classifying ground information into one of the crop categories. The second 
stage then classifies each into one of the health states. Assuming twelve 
crops and three states of health this is equivalent to a thirty-six-class 
pattern recognition problem. 

1.4.1 Baseline Data Format 

From Table 1.1. 2(1), entry AS, the average resolution for crop classi- 
fication is 40m, while the scanner swath width (field of coverage) is 185 km. 
From Figure 1.1. 2(1), the desired number of bands (channels) for this appli- 
cation is seven. These bands consist of the four ERTS bands, plus two 
bands in the raid-infrared range and a band in the thermal region. 

Using the data rate formula given in Section 1.1. 2. 4 with 

S = 185 km 
w 

R^ = R = 40 m 
L p 

V = 6500 m/sec 



we find that the data rate for this application is 

D = 3*16 X 10*^ bits/sec 
K 

This data format will be referred to as the "baseline data format" in 
subsequent sections. 

1.4.2 Preprocessing Algorithms 

1.4. 2.1 Coding 

From the results of data compression algorithms studies, there appear' 
to be two candidate algorithms applicable to on board data compression. 

They are: (1) transform coding (Section 1.3. 2. 1.1); and (2) Coding by BLOB 

(Section 1.3. 2. 1.2). Both of these algorithms are non- informat ion pre- 
serving in the sense that an exact replica of the original data cannot be 
constructed from the coded data. 

The studies further indicate that coding by BLOB achieves a higher 
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compression ratio than that obtained by transform coders, and this suggests 
that this approach is a better candidate algorithm for on board processing 
purposes, though it involves more computation. 

To determine the performa,nce of this method of data compression experi- 
mentally two sets of data were considered [2], One was aircraft imagery 
(Laboratory for Application of Remote Sensing CLARS) Run. #71053900) and the 
other was ERTS-1 imagery (LARS Run #73041801). 

BLOBS are larger and fewer at high significance' levels. The' contour 
tracer output contains fewer initial point locations and gray levels that 
must be coded resulting in a higher compression ratio. Figure 1.4. 2(1) 
illustrates the effect of various significance levels on the compression 
ratio . 
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Figlire 1.4* 2(1) Data Compress^ion Ratio for BLOB'-Contour Encoding Scheme* 
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% CLASSIFICATION ACCURACY 


To determine the effect of data compression on classification accuracy ’ 
in both these images , the blobs of the decoded data were classified using a • 
minimum distance classification algorithm C2], The classifier was imple- 
mented to distinguish corn, forage, soybean, forest and water in aircraft 
imagery, and wheat, pasture and other in ERTS-1 imagery. The results are 
shown in Figure 1.4. 2(2). Presently, machine classification of multispectral 
data images is most commonly done on a point-by-point (spectral signature) 
basis. Both data .sets were also classified according to their spectral sig- 
nature analysis. It is evident from Figure 1.4. 2(2) that compression ratios 
of the order of 12 and 32 for ERTS-1 and aircraft imageries, respectively, 
can be obtained without degrading the classification accuracy. 


'PER POINT' (SPEC. SIGN.) 



Figure 1.4, 2(2) Compression Ratio vs Percent Classification Accuracy 


1.4. 2. 2 Feature Selection 

The • candidate algorithm for feature selection is proposed to be based 
on eigenvector transformation discussed in greater detail in Section 1.3. 2.2. 
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To illustrate' the effect of this transformation on ERTS-1 data, an estimate 
of the' 4x4 covariance matrix of the data was determined. The' data was 
taken from an agricultural scene in Ogle, Lee and DeKalb Counties in the 
midwestern United States (LARS Run #72032806), The eigenvalues, matrix of 
eigenvectors and the variance in the individual channels is given in 
Table 1.4. 2(1). 

Table 1.4. 2(1) Variance Eigenvalues and Eigenvectors of ERTS-1 Data 


Variance 

Eigen values 

, Matrix of' Eigenvectors 

Channel 1-14.3 

124.39 

0.0491 

-0.0327 

0.8255 

0.5612 

Channel 2~27,6 

41.97 

0.5602 

0.8024 

0.1143 

-0.1705 

Channel 3-85.9 

1.97 

; 0.1927 

-0.1147 

-0.5526 

0.8027 

Channel 4-41.7 

1.19 

‘ 0.8040 

-0.5846 

0.0023 

0,1079 


Clearly the first two principal components of the transformed data 
contained 98.14% of the total variance contained in all four spectral bands 
of the original data. A set of training fields was selected to train the 
classifier on the data from the first tv/o principal components to classify 
the data into three classes (corn, soybean and other). A scatter plot of 
the first principal component versus the second principal component is 
given in Figure 1.4. 2(3). The decision boimdaries shown yield a classifi- 
cation accuracy of 98% determined on the training fields. 

This indicates that the eigenvector transformation of the ERTS-1 data 
reduces the data bulk by two orders of magnitude without any loss of infor- 
mation. In addition, since the most information is contained only in two 
dimensions, handling and processing of the transformed data becomes an 
easier task. 

1 .'4 . 2 . 3 Conclusion 

Although preprocessing plays a central role in data transmission, we 
have concluded that its usefulness is limited for on-board processing. The 
increase in hardware and software required to carry out the coding opera- 
tions is a necessity if data is to be transmitted from the host system to 
some processing station. In most applications, especially in agricultirre , 
the actual classification can usually be best carried out directly on the 
input data. For this reason, attention will be focused on the classifica- 
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• Figure 1.4.2 (3) A Scatter Plot of First Principal Component vs the 
Second Principal Component of ERTS-1 Data 

tion aspects of an on-board system for agricultural applications. 

1.4.3 Analysis Algorithms 

The analysis algorithms which are most suitable for on-board agricul- 
tural data processing are the maximum-likelihood and table look-up approaches 
discussed in Section 1.1.3. 

Attention will be focused in Section 2 on the hardware and software 
requirements of these two algorithms; as will be .seen, the two approaches 
pose quite distinct implementation, problems. 
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2 ON-BOARD PROCESSOR REQUIREMENTS 

This section contains the computational requirements of the data analysis 
algorithms, a review of some computer architectures and organizations, a 
design of some computer architectures capable of handling the algorithm 
computational requirements, and a discussion of the on-board processor en- 
vironmental effects. 

In Section 2.1 we survey some earth resources user requirements and the 
.data analysis algorithms for implementing them, in addition, we- analyze 
in detail the basic computational requirements; i.e., the number of multi- 
plications, additions, etc., required to implement these algorithms for the 
baseline data format. 

In Section 2.2 we present some computer architectures and organizations 
with particular emphasis on pipeline, array processors and multiprocessors, 
since it is apparent that some sort of parallel processor will be required 
to keep up with the high data rates required by the users. Memory and soft- 
ware requirements are also discussed. A number of on-board processors were 
then designed to efficiently implement the maximum likelihood and table look- 
up algorithms at the required rates. 

Finally, the environmental effects on the on-board processor for both 
the earth- synchronous and the sun-synchronous orbits are discussed in 
Section 2.3. 

2.1- ALGORITHM COMPUTATIONAL REQUIREMENTS 

This section deals with the basic computational requirements of pre- 
processing and analysis algorithms. It will be assumed throughout the 
following discussion that the estimation of the parameters for the eigen- 
vector transformation of Section 2. 1.1. 2 and the training of the maximum- 
likelihood and perceptron algorithms of Section 2 . 1 . 2 . 1 are carried out on 
the ground . 

The following notation is used throughout this section: 
n: number of features Cchannels) 

2 ^ = (Xj^, ..., pattern vector 

M: number of pattern classes 

Npi total number of picture elements Cpixels) in a data frame 

N, : total number of additions 

A 

Njj! total number of multiplications 
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N^; total nurabey of comparisons- 
Ng : total number of .square roots ■ 

2.1.1 Preprocessing Algorithms 

The preprocessing algorithms considered in this section are the data 
bulk reduction procedures' discussed in. Sect ion 1.3.2, The two principal 
approaches discussed in that section are based either on an eigenvector 
transformation or BLOB coding. The computational requirements for these 
two appraoches are discussed below. 

2.1.2 Basic Computational Requirements 
Cl) Eigenvector Transformation 

The reduction of a pattern vector with q^ components to n components by 
means of an eigenvector transformation requires nq additions and nq multi- 
plications. For a total of Np pattern vectors the following totals apply: 

= nqNp (2.1. 2-1) 

= nqNp ' (2.1. 2-2) 

where N^ and are the total number- of additions and- multiplications, res- 
pectively. 

(2) BLOB Coding 

The BLOB coding algorithm operates on two rows of an image at a time. 

V 

Letting W represent the number of columns per row we have the following 
processing requirements : 

(a) Additions in calculating the mean - 4nW 

(b) Additions in calculating F-test = 0 

(c) Additions in calculating t-test = 5nW 

(d) Multiplications in calculating mean and variance - n(n+l)W 

(e) Multiplications in calculating F-test - 2nW 

(f) Multiplications in calculating t-test - 6nW 

(g) Additions in storing mean and. variance - nW t 

Then, the total number of additions required to process two rows of an 
image is given by 

additions = ClOn + ^ 

X 

and the number of multiplications is given by 
multiplications = [8n + n^n+l)J W. 
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If we let.V equal the numhej? of rows in an image th,en we have. the following 


total figures : 

N^.= ClOn + (2. 1,2-3) 

= [8n + nCn+1)] ~W (2. 1.2-4) 

However, VK = N^, the total number of elements in the digital image, so that 
" I C2. 1.2-5) 

i n(n+l)3 Np (2. 1.2-6) 

2.1.3 Storage Considerations 


The storage requirements of the eigenvector transformation are simply 
those associated with the q x n transformation matrix. The transformation 
program itself consists of nothing more than the instructions required to 
multiply -the transformation matrix by each input pattern vector. Therefore, 
the total storage requirements are qn words for the matrix plus that re- 
quired to store the short multiplication program. 

- Experiments conducted at LARS with the BLOB algorithm indicate that 
approximately 35,000 32-bit words, including data storage, are required for 
coding with this approach. 

2.1.4 Analysis Algorithms 

2.1.4. 1 Basic Computational Requirements 


The operations required to implement the maximum-lilcelihood , percep- 
tron, clustering, BLOB, and table look-up algorithms are examined in this 
section. Expressions for the number of arithmetic and other required opera- 
tions such as comparisons and square roots are derived, The results are 
then compared in Section 2.1. 2. 2 using typical parameter values. 

(1) Maximum-Likelihood Approach - Classification of a pattern vector 
= (Xj^, ..., into one of M classes by the maximum-likelihood 

approach requires implementation of Eq. (1.1. 3-5) 

dj^Cx) = An p(Wj.) - /^n j^l - j C(x-I\)'' (2. 1.2-7) 


k-=.l, 2i 


M 


The terms ' lii p Co), ) , ™ Jin | C [ ; » ^ ' are' estimated during ‘ training and . 

• X ^ iC — K 

consequently, need not be' calculated 'during the classification process. 
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Therefore, the principal computation for each, inyolyes multiplications 
and additions, in .the last term-.of ’Eq. . (.2.1.2-'7) aboye. 

For n. channels, x and m haye.n components and ^ ^ is an n x n matrix. 
Thusj the number of arithmetic operations is determined as follows: 

(a.) formation of (x-m, ); n additions 

-1 2 

(b) multiplication of by Cx-n^): n(n-l) additions; n multipli- 

• cations • 

(c) multiplication of Cx-nu ) by the' results of (b): Cn-1) additions; 

n multiplications 

(.d) multiplication of 0.5 by the results of Cc); 1 multiplication 
(e) addition of all three terms in Eq. C2. 1.2-7): 3 additions 

'The results are, therefore, 

2 

number of additions = n + n + 2 

2 

nvraiber of multiplications = n + n + 1 

These results are for one class and one pattern. Letting M equal the number 

of classes and Np the total number of patterns to be classified yields: 

= (n^ + n + 2) M Np C2. 1.2-8) 

H„ = (n^ + n + 1) M -N„ (2.1. 2-9) 

M P 

In order to classify each observation x into one of M classes it is 
required that the inaxiinum of M decision functions be determined as was 
indicated in Section 1.1^3*1, Therefore ^ in addition to the total number 

/ 

of operations given in Eqs. (2.1. 2-8) and (2. 1.2-9), it is also necessary 
to perform comparisons, where 

= Np(M-l) (2.1.2-10) 


This relationship follows from comparing M decision function values for Np 
patterns . 

(2) Perceptron Approach - Classification of a- pattern vector into one 
of M classes by the perceptron approach requires implementation of Eq. 
(1.1. 3-6) 


= "k" i (2.1.2-11) 

k = 1, 2, .. ., M 

where the coefficients aiid functions {^^(x)} are determined during 

training. 
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In. order. to establish,. the number of arithmetic operations required to 
implement the.perceptron approach.it is. necessary to choose a particular set 
of functions One of. -the most common approaches is to express the 

above decision functions in'. the < form of an rth degree pol 3 niomial [1] by 
means of the following recursive formula; 


d'^Cx) = 


n 

n 

Z 

2 

^P^=l P 

Vs. 

2"?1 

+d^-^ 

(x) 


n 
S 

P =P- , 
r r-1 


w. 


P P . 
12 


.P 


X. 


P ■ 
1 ^2 


.X, 


r 


C2. 1.2-12) 


where the first term d°(x) ’is given by d^(x) = this equation the 

superscript indicates the degree of the polynomial. The class subscript has 
been dropped from the equation-'to simplify the notation. It is understood 
that there are M decision functions, one for each class. 

The number of additions in Eq. (2.1.2-12) is given by: 


additions = 


_ (n + r)I 


nJ r I 


- 1 


and the number of multiplications by: 


multiplications = 


.'r 

2 (q+1) 

q=l 


' (n+q) ! 
nl ql 


(n-hq-1) ! 
nl (q-1)! 


-2 


These results are^ for one pattern and one class. Letting M equal the num- 
ber of classes and Np the total number of patterns to be classified results 
in the following total number of arithmetic operations: 


N, 


fCntr) ! _ 1 

= l^nl rl J 


MN,. 


-2 


MN,. 


(2.1.2-13) 

(2.1.2-14) 


('n+ 1 'k r (ntq)! (ntq-l)I 

njq. n!(q-l)I 

As in the case of the maximum-likelihood approach, it is also necessary to 
perform comparisons, where 

= Np(M-l) (2.1.2-15) 

(3) Clustering Approach - The clustering approach discussed in Section 
1.1. 3. 4 operates directly on the complete set of data to be classified. 

This is in contrast with the two approaches discussed above which classify 
one pattern vector at a time. The number of operations required to imple- 
ment the clustering approach is determined as follows: 


(a) sample mean for each class M, p. ^ 


L 2 X. ; j = 1, 2, 


i=l 


— 1 


n 
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.additions:' n L M 
multiplications ; > n ' 

where L = Np/M. ' It is noted. that M, the -class or mode -of a .cluster, is not 
specified' as in. the case of the two approaches discussed above, but is in- 
stead determined by the algorithm based. on the properties of the data being 

analyzed. ^ ■ 

2 1 ^ ’2 

Cb) sample variance O. = 2 (x. y.) ; j = 1, 2, ..., n 

3 13 3 


additions: 2 n L M; multiplications: n(L+l) M 


(c) mode Cclass center ra 


k3 


, -f2Ck-lf 

- Uj + cT 

U- ^ 


1; j = 1, 2, 


n 


additions: (3M+1) n; multiplications: (3M+L) n: k = 1, 2, 

n ' 2 

C = 2 (x . -nw .) k = l, 2, ..,,M 

• T X KX 
X=1 


M 


Cd) Euclidean distance 


additions: n Np M; multiplications: n Np M 


Let I be the number of iterations allowed for convergence of the algorithm. 

Then, the total number of computations involved in assigning all the data 

vectors to one of M classes is: 

additions = (n L M + 2n'‘"li M + (3M + l)n + n Np M) I 

= (n Np + 2n Np + 3n M + n + n Np M) I 

= n itM + 3)N_ + 3M + 1] 

ir 

multiplications = (nM+nNp+nM+ 3Mn + n + n Np M) I 
= n IC5 M + CM + 1) N„ + 1] 


to determine the distinctness of classes we have the following additional 
computations: , 

Ce ) Swain-Fu distance : V = 






where = Tr ^(m^-m 2 ) (ni-mg . Calculating 


-1 3 

C, involves n additions 

3 

and n''' multiplications. Then, the calculation of A, involves 2n + 1 addi- 

3 2.. ' ^ 

tions and 2n + n multiplications. Finally, the calculation of y for two 

classes involves: 

3 3 2 

additions: 4n +3; multiplications: 4n + 2n +2; 

square roots : 3 

Considering next every pair of classes . yields : 

additions: t 3]; multiplications: [4 un 4 - 2n^ t 2] 


2.. 


square roots: 


'smCm-D 
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The -total number of computations .required for clustering is. then:' 


= n I -E(M + -3)Np t 3M + 1 t [4n.^ + 33 

= n I E5M + Qi+1)N„' + 13 +'H|^ E4n^ + 2r? + 23 

n ‘ “ P 2 , 

M = 3MCM-1) 

S . 2 

where Ng is the number of square roots. 

In addition it is also necessary to perforin 
= NpCM - 1) 


•C2. 1.2-16) 
(2.1.-2-17) 
C2', 1.2-18) 


(2.1.2-19) 


comparisons . 

C4) Blob Classification - After the data have been compressed by the 
BLOB preprocessing algorithm, it may be desired to classify each blob into 

one of M classes. Two principal approaches are used in the BLOB classifi- 
cation algorithm: Cl) minimum-distance classification, and C2) maximum- 

likelihood classification. Letting represent the number of blobs, we 
have the following figures for the minimum-distance approach: 

=(2n^ + n^ + n + 2) MNg (2.1.2-20) 

= (2n^ + n^ + n + 5) MN^ - (2.1.2-21) 

= Kp(M - 1) (2.1.2-22) 


In addition, there are N^m log operations. 

The maximum- likelihood figures are obtained from Eqs. (2. 1.2-8) 
through (2.1.2-10) with = Ng. 

N. = (n^ + n + 2) MN_ (2.1.2-23) 

A 13 

= (n^ + n + 1) (2.1.2-24) 

M B 

Nq = Ng(M - 1) (2.1.2-25) 

(5) Table Look-Up - The principal operation in the table look-up 
approach discussed in Section 1.1. 3.3 consists simply of 

= (11 + M) HNp (2.1.2-26) 

additions. Although this technique is attractive from a computational 
point of view, it must be kept in mind that the actual program and storage 
requirements are considerably more complex than the approaches discussed 
above . 

2. 1.4. 2 A Numerical Comparison 

The results obtained in the previous section are best appreciated by 
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comparing the number of require^i operations on an absolute basis. Table 
2. 1.1 CD presents; such, a comparis.on for the .following, typical parameter 
values': ‘ 

n = 4 . 

0 

Np = (3300 rows) x (2500- columns) = 8.25 x 10 pixels 
M = 5 

I = 15 • • 

Number of points/blob = 10.' Then Ng = N^/IO = 8.25 x 10 
These values are typical of those used in. connection with ERTS-1 MSS data. 


’"Table 2. 1.1(1) A comparison of Analysis Algorithms 


for One Data Frame 


Equations Number 
Used of 

Approach (Section 2.1.2) Additions 


Maximum-Likelihood 

(8), (9), (10) 

9. 

075 

X 


Perceptron 

(13), (14), (8), 

5. 

775 

X 

10*^ 

[r=2], (15) 




9 

Clustering 

(16,(17), 

3. 

960 

X 

10^ 


[1=15] „(18), 






(19) 




8 

BLOB [Min. Dist.] 

(5), (6), (20), 

8. 

251 

X 

10 


(21), (22) 




9 

BLOB [Max-Likeli- 

•(5), (6), (23), 

1. 

114 

X 

10 

hood] 

(24), (25) 




8 

.Table Look-Up 

(26) 

6. 

,600 

X 

10° 


Number of 
Multipli- 
cations 

8.663 X 10 
1.485 X 10^ 

2.970 X 10^ 

8.044 X 10^ 
1.040 X 10^ 


Number 
of Other 
Operations 


N : 3.30 X 10^ 
N^: 3.30 X 10 


Ng: 30 

Nj,: 3.30 X 10 


7 


Log: 4.13 X 10 
N„: 3.30 X 10^ 
Ng: 3.30 X 10^ 


.The results obtained in Table 2. 1.1 (I) for the perceptron approach 
assume decision boundaries of the second degree (r=2). This assumption is 
made to arrive at a meaningful comparison between the perceptron and maxi- 
mum'- likelihood approaches 5 since the latter implements second degree 
boundaries when the data are assumed to have Gaussian properties. This is 
by far the most common assumption followed in processing remotely- sensed 
data by the maximum-likelihood approach. 

The results shown for the BLOB approach take into account Eqs. (2. 1.2-5) 
and C2. 1.2-6). Although these equations are associated with preprocessing 
functions, they are also required -in classification. In other words, pre- 
processing is a necessary step' prior to classification by the blob approach. 

\ 
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2.x. 4. 3’ Storage Considerations 


Implementation of the maximum-likelihood algorithm for M classes re- 
quires M storage locations for to pCw^^), M storage .locations* for ^ » 

Mn storage locations for the mean vectors m^^ and Mn^ storage locations for 
the' matrices ^ The actual classification algorithm can be implemented 
in less than 2500 storage words. In addition, it is necessary to have a 
buffer for the input data. The size, of the buffer depends on the processor 
speed. 

The perceptron algorithm requires M(n+r)I/nlr! storage locations for 
the coefficients of the decision functions. The algorithm itself can be 
stored in less than 2000 words. A buffer is also required for this 
approach. ■ 

Experiments conducted at LARS with the clustering algorithm indicate 
that 57,000 32-bit words are required to process 40,000 data vectors (with 
n=4). This includes program and data storage. In -addition, a buffer whose 
size depends on processor speed is required'. 

The storage requirements for the classification aspect of the BLOB 
algorithm depends on the approach taken. For the minimum-distance approach 
it IS necessary to store the. mode centers along with the covariance matrix 
of each blob. Since the number of 'blobs varies with the application the 
storage requirements must be determined by actual experimentation . • The 
classification algorithm itself can be stored in less than 2000 words and 
a buffer is also required. 

The maximum-likelihood approach to blob classification has the same 
requirements as the general maximum-likelihood algorithm discussed above. 

Experiments with the table look-up approach indicate that this approach 
requires in the order of 31,000 words of storage [5,6]. This figure 
applies to a maximum of 24 classes with n=4. A buffer whose size depends 
on processor speed is also required. 

The above comparisons are intended only as rough guidelines of algo- 
rithm storage requirements. The actual figures for each algorithm depend 
on a number of factors such as processor instruction set and speed, program 
optimization, and the application area. 

2. 1,4. 4 • Determination of the' Algorithm Constants 

The constants in the preprocessing and analysis algorithms depend on 
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the particular application. For example, rhe etc., in the maxi- 

mum likelihood algorithm- decision functions given by equation (2. 1. 1-7 )depend 
on the statistics of the. data classes to.be classified. The on-board pro- 
cessor, must implement this equation, but the processor design, speed of 
operation, etc., are invariant to the numerical values of the constants. 

There are two possibilities for determining the numerical values for 
'these constants: (1) the spectral signature bank approach; and (2) training 

The' signature bank approach assumes that the required numbers for each 
application are known and stored so that for any particular application the 
required set of numbers can be looked up in the bank. This approach is not . 
feasible at this time because the constants for any particular application 
depend on many different parameters in a complicated way that is not well 
understood at this time. For example, the spectral signature of corn de- 
pends on lighting conditions, soil characteristics, the amount of moisture 
in the leaves, the point in the growing cycle, etc. When, if ever, the 
. dependence of the constants on these parameters becomes well enough under- 
stood to allow application of this approach is uncertain. 

„jrhe'. second approach is, called training .and. requires^ that a sample of 
each type of data to be classified be obtained and the required numerical 
values measured from this data sample. For example, the in the maximum 
likelihood algorithm are the mean intensity values in each spectral band. 

The idea is to choose a typical data sample (called a training sample ) , 
make the measurements from this sample, and use them in the decision func- 
tion. 

If the training approach is used there are two possibilities^ for 
training the classifier ; i . e . , obtaining the numerical values of the con- 
stants. The first approach is to collect the data sample and compute these 
parameters on board the satellite. The second’ approach is to transmit the 
raw data sample to a ground computer, compute the parameters on the ground, 
arid transmit these back to the satellite. The latter approach appears to 
be the most feasible at this time . The entire problem of how to go about 
collecting training samples, verifying ground truth, etc. , is not well under 
stood at this time. NASA is presently expending considerable resources in 
an attempt to solve these problems. The results of these studies will com- 
plement the results of oin? study. 
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2.2 COMPUTER ARCHITECTURE AND ORGANIZATION 
2.2.1 Historical Data 

Historically, computer organization has developed along the classical 
lines of the Von Neuman machine Cl]. When the arithmetic logic unit (ALU) 
is designed to operate on n-bit operands in a time-sequential mode, the 
operation of the computer is termed serial mode. Traditionally, a computer 
whose ALU is serial has been classified as a serial computer. If the ALU 
execution on an n-bit operand is performed separately and simultaneously, 
the mode of operation is parallel. Approximately n times as much hardware 
is required in a parallel organization as compared to a serial organization. 
This is somewhat offset by the fact that the serial control unit is more 
complicated than that of a parallel organization. However, the parallel 
structure executes operations nearly n times faster [2,3]. 

First generation computers were serial. The perforraance/cost ratio 
was relatively small. To increase the speed of operation, second genera- 
tion machines were designed with a parallel structure. Transistor and 
better random access memory technologies significantly increased the per- 
formance/cost ratio. The user base broadened with the development of 
higher' level programming languages 'which', generated a demand for more com- 
puting power. 


71 



Third generation computers . use integrated circuit technology and im- 
proved . peripheral and memqry.techhdlogies.to,continue .the improvement in the 
performance /cost ratio. . .The. organization of these computer systems still 
uses' a parallel'-organization. ' System' resources are utilized. more efficiently 
in a few of these systems by . adding a second processing unit. In this 
organization, independent processors, each with identical physical organiza- 
tion, performs relatively independent functions. The memory unit is a 
resource shared' by each processor . Although the processor cost has in- 
creased two-fold, the performance has increased by a lesser factor. Even in 
third generation machines' the' Von Neuman influence is evident in most com- 
puter’ organizations. • Flynn E4] and Higbie [6] have presented new concepts 
which attempt to avoid the ambiguous term "parallelism" in categorizing 
computer. organizations. 

Four major classifications for computer processors have been defined: 

Cl)' Single Instruction Stream-Single Data Stream (SISD) processors; 

C2) Multiple Instruction Stream-Single Data Stream (MISD) processors; 

(3) Single Instruction Stream- Multiple Data Stream (SIMD) processors; 

(4) ,. Multiple- Ins.truction, Stream-Multiple Data Stream- (MIMD) processors. 

The Instruction Stream is the sequence of instructions as executed by 

the processor while the Data Stream is the sequence of data called for by 
the instruction stream . Computer organizations are described by the multi- 
plicity of hardware subsystems utilized to service the Instruction and Data 
Streams. No longer is a system structure identified by the physical con- 
struction of its ALU, but rather by the number of instruction streams, the 
number of data streams, and their interactions. 

2.2.2 Parallel Processor 

SISD Computer . This system has a single processor. -Several computer 
systems using this basic organization achieve computing power by overlap- 
ping the execution of sequential instructions. That is, the execution of 
the next instruction begins as soon as the first phase of the current in- 
struction is completed. When conditional instructions are encountered, 
"dummy" instructions are executed' until the decision or conditional pro- 
cess is complete. The CDC 6600 series [5] and the IBM '360/90 series [7j 
are -examples of computers falling into this class. 
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2.2.3 Pipeline > Processor ' 

MISD CoTTiPuter .. This structure has. become known as a. pipeline computer 
The whole philosophy of pipelining' is to' implement a series of 
independent processors, ' or processor subsystems j, in the form of a long 
pipeline. As a computation flo>7S through the collection of processors, each 
processor performs a specific sub-operation.' Several computations may be 
in process in distinct processors at the' same time. The overlapping of 
computations is made possible by locating output registers with each segment 
of the pipeline. The' CDC STAR [1Q3 and the TI-ASC £113 are xhe most signi- 
ficant examples of computers in this class. 

2-2.4 Array Processor 

SIMP Computer . Several computers (SPAC, [12] SOLOMON, [13] and 
ILLIAC IV [14]) have been proposed to process multiple data streams. ILLIAC 
IV is the only one of consequence to have achieved an operational status . 
Computers falling in this class of structure are described alternately as 
array- processors, associative array processors, parallel processors, vec- 
tor processors, or othogonal processors depending upon their organizational 
differences. Communication between processors is predetermined and instruc- 
tion stream execution operates simultaneously on all data streams. No 
overlapping of instruction execution is presently used in this organization* 

2.2,5 Multiprocessor 

MIMD Computer . This is an example of a multiprocessor system with 
either shared or multiple memories in which more than one program is 
executing on its own data multiprogramming). The Carnegie-Mellon 

system* [15] -is an example of such an organization and employs a number of 
independent minicomputers with a shared memory. 

Looking at these organizations with respect to the data acquisition 
‘ from a multispectral scanner, an obvious observation is that data acquired 
from the multispectral scanner represents a multiple data stream. A SIMD 
organization is one possible consideration. However, to achieve the data 
rate-constraints, a MIMD structure represents an alternate approach. 

2.2.6'* Memory- 

In a small general-purpose digital computer the memory is the most 
expensive of the various computer subsystems. The memory size must be in 
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excess of the -requirejnents .set .forth. by. the average job that is anticipated, 
must '.be fast to prevent a limiting of the , cojnputer performance, and must be' 
low povter'.to reduce .power .supply .costs. . Small’ computers today are turning, 
almost . exclusively , to semiconductor .memories ..which are out-performing core 
in both speed and. cost for small to. raedium-'Sized. memories Cl6,17!]. 

. The architecture of contemporary machines has. included a single memory 
which is used for both, program storage and temporary storage . ¥ith the 
advent of the. microprocessor and the semiconductory memory, the memory has 
been partitioned into two separate memories, program memory and working 
memory- Because of the volatility of semiconductor memories, the program 
memory is a read-only memory CROM) and the working memory is a random- 
access memory CRAM). 

2-; 2 . 6 . 1 Program Memory 

In examining the algorithms for earth resources applications the pro- 
grams are not lengthy and complex since speed is the predominating factor j 
therefore, the program memory can be relatively small. Projections 
indicate that the on-board computer will be of the multi-instruction 
multiple-data- input type in the 1980-1990 time span. To implement this 
type of architecture, one method is to interconnect a number of micropro- 
cessors; one or more would be in each stage of the pipeline, and to accom- 
plish the throughput requirements, pipelines could be parallelled. Each 
microprocessor or group of microprocessors could be assigned a small program 
memory. At’ present the access time for a semiconductor memory is much 
faster than the cycletime of the microprocessor so that something like five 
microprocessors could share a memory. Since both microprocessor speed and 
memory speed will increase, this ratio may still stand in the 1980-1990 time 
frame . 

The program memory must be non-volatile and there must be a convenient 
way to change the program. Today, this is accomplished in programmable 
read-only memories (PROM). Electrically-alterable read-only memories (EAROM) 
will be available soon. The common type .of PROM today requires an ultra- 
violet light source to erase the raeroroy. By 1980 the bits’ of the EAROM can 
be' selected randomly for writing into either state. 

2 . 2 . 6 . 2 ' Working Memory ’ 

The working memory criteria can be satisfied by a random-access 
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volatile memor>y (RAM). P;resently these -are either, bipolar or .jnetal-oxide. 
semiconductor CMOS). The", bipolar -.RAMs enjoy on the' average, a speed advan- 
tage and the.MOS RAMs are. superior in their efficient use of power. Currently 

2 . 

the' bipolar chips contain 4000 bits, .require 6 mil /bit and an access rime 
of 50-100 ns.' This performance- is expected. to improve. at a .substantial rate 
through the 1980's. 

2. 2. 6. 3 Bulk Memory 

The on-board earth resources computer will not need an on-board bulk 
memory to store the MSS signals for later_ processing. An on-board computer 
is being designed to operate in .real time and to immedi3.tely process all the 
data as.it is produced' by' the sensors. 

The output of the computer contains processed data in the form of 
classification data for each pixel. For the candidate format discussed in 
Section 1.4 there are four spectral bands with 6 bits of gray-level quanti- 
zation -for each pixel. There are 1.3 ys per pixel so that the input data 
rate is 31.2 MHz. 

The ou-tput data is classified into 12 classifications or four bits per 
pixel. This- gives an output data rate of 5.2 megabits per second. At this 
point the on-board processor has compressed' the data by a factor of six to 
one. If these output data are to be preserved, an on-board magnetic tape 
recorder would be required or a wide band data link would be necessary 
between' the spacecraft and the ground. Another approach is to sort the 
classification data into twelve registers aboard the spacecraft. The regi- 
sters could be read out at the end of each scan line indicating the combined 
'type and quality of crop along the scan line. This greatly reduces the 
output data rate but sacrifices the knowledge of the exact location and 
quality of a particular crop. 

2. 2. 6. 4 New Memory Technologies 

A number of new technologies TlTj have been developed in the memory 

area. One of the most promising is. the silicon-on-sapphire complementing 

MOS memory element. These circuits offer one very important advantage. . 

The insulation resistance between the components is extremely high and the 

> 

decrease in . s-uBstrate parasitic capacitance is so great that it appears 
that MOS circuits built on sapphire could be as fast as or faster than 
bipolar integrated circuits'. 
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The metal-alumlna-silicon field-effect transistor is especially well- 
adapted for use as an electronically-alterable read-only memory. This 
transistor is the same as a conventional MOS except that the alumina is used 
in place of oxide as the gate dielectric material. The characteristics of 
the alumina dielectric material can be changed' by applying a voltage to the 
gate. There is a critical gate voltage level which, when reached, alters 
the transistor’s threshold voltage. Shifting the threshold voltage from 
one level to another corresponds to the storage of one bit of information. 
Fortunately, the change in threshold voltage is non-volatile allowing the 
device to be especially suifed as a program memory. ■ 

Charge- coupled devices appear to have their’ greatest potential in re- 
placing disk and drum storage since the device works an excellent shift 
register. Presently there appears to be no need for high-capacity shift ' 
registers for on-board earth- Resources computers. 

For severe environment applications’ the Metal-Nitride Oxide Semiconduc- 
tor (MNOS) has shown great promise coupled with the fact that it is a non- 
volatile memory element. The devices have been applied in a large memory 
which replaces disk or drum Cl93. An MHOS memory could serve either as a 
program memory or a working memory or both. 

2.2.7 Software 

The requirements on the software for the on-board processing of earth- 
resources data are basically different from the requirements of ground-, 
based processing software. This stems from the fact that the ground com- 
puters have traditionally been large-scale general-purpose digital computers. 
These computers are used for many tasks and may even be assigned to com- 
pletely different roles at the conclusion of the earth-resources mission. 

The on-board processor is dedicated to a specific task and more of the 
system implemented in hardware than in software. This is complemented by 
the speed advantage of hardware over software. 

A little over a decade ago the software costs were less than the costs 
of the hardware. With the introduction of Large Scale Integration (LSI), 
the cost of the computer has declined considerably; however, the labor 
costs for software design has continued to rise. Some estimates put the 
cost of a single instruction in a complicated program in the neighborhood - 
of ten dollars; another figure is that -one hour of labor is required for 
each two instructions in a program [20], The trend, then, is to put more 
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software functions -into hardware. For the* earth-resources on-board compu- 
ter this trend is emphasized even more strongly. Since speed is the pre- 
dominant factor, as many of the computer functions as possible will have to 
be done- in hardware. In early computers binary addition was the only 
function implemented in hardware. To increase speed, more and more functions' 
have been implemented in hardware. 

.The on-board earth-resources > computer does not 'need all^the features 
of a large general-purpose machine^, but it does need speed, !t^ seems practi- 
cal, then, to do the time consuming jobs by hardware and *use software where 
speed is not essential so that flexibility can^ be maintained. As an 
example , the maximum- likelihood algorithm involves primarily additions^ and 
multiplications. It was very evident from the start that multiplication 
could not be done in software and that either hardware' multiplication 'or a 
table look-up multiplication method would be employed. The addition func- 
tion, of*' course, is -readily available in hardware form in all CPU chips or * 
arithmetic logic units.* 

' Since high speed is a prime requirement of the on-board computer, the 
programs must' be short, simple and efficient. Efficiency dictates th'e use 
of an assembly language for the programming. In the maximum- likelihood 
algorithm the total number of operations required is small; further, the 
pipelining technique breaks down the computation into individual tasks. 

The microprocessors are assigned to each task, and each task has its own 
program. These programs are very short since 'the total operation has been 
broken down* into parts (tasks). The individual programs would reside in an 
electrically-alterable read-only memory (EAROM) . By ground command to the 
spacecraft, coefficients could be changed in the program to match ground- 
truth “measurements . * ’ ’ ' . ‘ 

An advantage of utilizing identical microprocessors in the on-board 
computer is that programming becomes fairly standard. A cross-assembler 
would be used to develop each program for each task for loadinginto the 
EAROMs. ’The cross-assembler program would be run on any available computer 
of adequate capacity. A higher- level language would be of little use 
since it is not envisioned that a large’ amount of program will be needed; 
also, a higher-level language is not as efficient. 

* • m * 

2.2.8 Computer Architecture' Examples 

Attention, will be focused on implementations of the maximum-likelihood 
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and table look-up algorithms. Section 2. 2. 8.1 deals with a microprocessor- 
implementation of the maximum-likelihood algorithm. This is followed in 
Section 2.2. 8.2 with a microprocessor implementation of the table look-up 
approach. Section 2. 2.8.3 deals with a microprocessor implementation of an 
expanded form of the maximum-likelihood algorithm. This is followed in 
Section 2. 2.8.4 by a hardware implementation of this approach. Finally, 
Section 2. 2. 8. 5 presents a hardware implementation of the table look-up 
algorithm. All examples are based on the baseline data format defined in 
Section 1.4. 

2. 2. 8.1 Microprocessor Organization for the Maximum-Likelihood Algorithm 

Four processing elements are pipelined to execute the maximum likeli- 
hood function for one class as shown in Figure 2.2. 8(1). Each processing 
element (PE) is organized about a microprocessor; the performance measure 
is the execution time. A PE utilizing -the INTEL 8080 is shown in Figure 
2 . 2 . 8 ( 2 ). 



= < 5 .- 1 ^)' 

- 1/2 (x-;^) 

Figure 2.2. 8(1) Pipelined Organization for Executing 
the Maximum Likelihood Algorithm 

Input • Output 

Registers Registers 



Figure 2. 2. 8(2) Microprocessor System Organized 
About an INTEL 8080 Microprocessor. 
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The vector x. is input to PEI. PEI forms. x~j\ outputs the result to 
PE2 which generates (jc-n^) ^ and. outputs it to PE3 where is compute'd 

PE4 determines the classification category, k. 

Processing Element PEI . ’ 

The n input registers and nM output registers are required for n 
spectral bands and M classification categories. In order to determine the 
time required to compute x-mj^ the required code was written using the instruc 
tion set of the 8080 as illustrated in Table 2.2.8{I). This computation 
must be repeated for each of the m categories, k=l, 2, ..., M. The toral 
execution time is ' . 

TE = (13. 5M + 7.,5)n ys. 


Table 2. 2. 8(1) 
MNEMONIC OPERAND 


Program for Generating x-^ 


BYTES 


in PEI 
COMMENT 


LXI 

H, ARRAY 

3 

IN 

XI 

2 

MOV 

D,A 

1 

SUB 

Mil • 

1 

.OUT ’ ■ 

Xll 

2 

INX 

H 

1 

MOV 

A,D 

■ 1 

SUB 

M21 

1 

OUT 

X21 

2 

INX 

m 

• 

• 

• 

K 

1 

5(M-3) 

• 

MOV 

A,D 

1 

SUB 

MMl 

1 

OUT 

XMl . 

. 2 


initialize mj^j. a^ray 
input component of Ji 

generate. Xj^ - m^^^ 
output Xj^ - m^^ component 
increment HL 

generate x^^ - 

output Xj^ - component 

increment HL 

generate and output x^^ - m^^ 
component i = 3, . . , M - 1 

generate Xj^ - m^^^^ 
output Xj^ - component 


This segment of program is written for the x^ component of Jt, but 
must be executed for each of the remaining components x^, . . . j x^. 


Processing Element PE2 . • . 
PE2 computes the. product 




nM input registers and 2nm output registers are required. Table 2. "2. 8 (II) 

shows the 8080 code required to use table look-rup to form the products. 

-1 

and the coefficients of • are used to address a multiplication 
ROM from which the product is read. This code must- be executed for each 
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Table 2.2.8CII) 


Program for Generating 


-1 


in PE2j^ 


MNEMONIC OPERAND BYTES 


COMMENT 


LXI 


H, TEMP 

3 

LXI 


B, CROWl 

3 

IN 


Xll 

2 

OUT 


Xll 

2 

MOV 


D,A 

1 

LDAX 


B 

1 

MOV 


E,A • 

1 

LDAX 


D 

1 

MOV 


Tl, A 

1 

IfJX 


B 

1 

INX 

« 

« 

« 

H 

1 

6(n-2) 

LDAX 

« 

B 

1 

MOV 


E,A 

. 1 

LDAX 


D 

1 

MOV 


Tn, A 

1 

LXI 


H, TEMP 

3 

LXI 


B, CROWi 

3 

IN 


Xli 

2 

OUT 


Xli 

2 

MOV 


D, A 

1 

LDAX 


B 

1 

.MOV 


E, A 

1 

LDAX 


D 

1 

ADD 


T1 

1 

MOV 


Tl, A 

1 

INX 


B 

1 

INX 

• 

* 

* 

H ■■ ' 

1 

■ 7(n-2) 

LDAX 

* 

B 

1 

MOV 


E, A. 

1 

LDAX 


D 

1 

ADD 


1 

1 

MOV 


1, A 

1 

LXI 


, TEMP 

3 

LXI 


j CROWn . 

3 

IN 


Xln 

2 

OUT 


Xln 

2 

MOV 


D, A 

1 

LDAX 


B 

1 

MOV 


•E., A 

1 

LDAX 


D 

1 

ADD 


Tl 

1 

OUT 


All 

2 

INX 


B 

1 

INX 


H 

1 


initialize location of n tem- 
porary storage locations 
initialize row 1 in C^^ matrix 
input x^ - m^^^ component 
output Xj^ - m-j^ to PE 3 


generate product (x^ - 

increment BC 
increment HL 

generates products (x^ -m^ ^ ) . C^ . 
r— 2 y « • • ^ n— 1 


generates product 
(x‘ -ra, , ).Cl 


TnV U) 


In 


initialize row i in C, matrix 
input x.-m^ . component 
-output x.-m^^ to PE 3 


generates product (x.-m . ).C^ 
accumulate sum in Tl^ ^ ^ 


product (x. 

j = 2. ...t 


-m. . ) . C7 • 

n-J’- ^3 


fJroduct 


(x.-m, . 
1 li 



initialize row n in C,^ matrix 

input X -m, component 
1 . n In r.T^o 

output X -m. TO PE 3 
n in 


product (x -m, ).C^, 
^ n In nl 

output A~^ of 
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k=l, 2, ...jM. The program execution time in PE2 is 
TE = (21. 5n + 15.5)nM ys. 

Processing Element PE3~ . 

PE3 accepts the 2nm outputs of PE2 and generates 
“ 1/2 

with 2nm input registers and m output registers. A program is presented in 
Table 2.2.8(111). The execution time is 
TE = (23. On + 15.5)M ys. 


Table 2.2.8(111) Program For Generating ^ 
MNEMONIC OPERAND BYTES COMMENT 


H, CONS 


COMMENT 

initialize location of constant 


IN 

Xll 

2 

MOV 

D, .A 

1 

IN 

All 

2 

MOV 

E, A 

1 

ODAX 

D 

1 

MOV 

B,.A 

1 

IN 

,X12 

2 

MOV 

D, A 

1 

IN 

A12 

2 

MOV 

E, A 

1 

LDAX 

D 

1 

ADD 

B 

1 

MOV 

■Bj A 

1 

9(n-3) 

IN 

Xln 

2 

MOV 

D, A 

1 

IN 

Ain 

2 

MOV 

E, A 

.1 

LDAX 

D 

1 

ADD 

B 

1 

RAR. 


1 

MOV 

B, A 

1 

MOV 

A, Cl 

1 

SUB 

B 

1 

OUT 

D1 

2 


input x^ - component 

input A^ ^ ^ component of A^ ^ 
E (A)^’^ ^ 

product A~^. (x, -m, , ) 

B (A) ^ 


product A2.(x2-mi,) 

accum^llate sum in B 

AT^(x.-m . ) 
li._i li , 

1—0 ^ ^ 


product a 7 (x ) 

^ In n In 

divide by 2 

l/2(x-mj^)"C^-^ (x-ULl) 


generates d- (}^) 
outputs dj^ Ik) 


Processing Element PE4. 


PE4 compares the M values of dj^(jc) to determine the classification. 
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This requires m input registers and a single output register. Table 2.2.8 
(IV) describes the program for PE4. The time required is 
TE = (19. 5M - 2.5)ps 


Table 2. 2. 8 (IV) Program for Determining Classification k in PE4 


MNEMONIC 

OPERAND 

BYTES 

COMMENT 

MVI 

C, 1 

2 

c ■«- 1 

MVI 

E, 2 

2 

E ■«- 2 

IN 

D1 

2 

input d- (x) 

MOV 

B, A 

1 

B contains largest 

IN 

D2 

.2 . 

input (x) 

CMP 

B 

1 

compare d^ (x), dj (x 

JC 

ALPHA 

3 


MOV 

B, A 

. 1 

replace d^ (x) , 1 

MOV 

C, E 

1 

with d 2 (x) 9 2 

ALPHA: INR . 

E 

1 


IN 

D3 

2 


* 

LAST: MOV 

A,C 

i 


OUT 

k *: 

f ^ 

2 

I 

Output' classification 



- 



Table 2.2.8(V) presents the comparison of current and projected pro- 
gram execution times for each processing section for n and M eq^ual to 4 and 
12 respectively. The 1974 column and 1983 column result from data taken 
from the median curve of Figure 3. 2. 2(1) for projected microprocessor add 
times.. 


Table 2.2.8(V) Comparison of Current and Projected Program 
Execution Times for the Processing System of Figure 2.2. 8(1) 



Number of 
Microprocessors 

8080 

Microprocessor 
3.5 ys ADD 

1974 

Microprocessor 
1.0 ys ADD 

1983 

Microprocessor 
100 ns ADD 

PEI 

1 

678.0 ys 

193.7 ys 

19.37 ys • 

PE2 

1 

4872.0 

1392.0 

139.20 

PE3 

1 

1290.0 

368.6 

36.86 

PE4 

1 

231.5 

66.1 

6.61 


The bottleneck in 

the pipeline is PE 2. 

PE2 can be split 

into M parallel 


processing elements as in Figure 2. 2.8(3). Each (J^ = Ij •••» M) is 

input to a separate PE2, and the computation for each category proceeds in 
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Figure 2. 2. 8(3) Organization of a Processing System for 
Executing the Maximum Likelihood Algorithm. 

parallel. Table 2. 2. 8 (VI) gives the execution times for M=12. 


Table 2. 2. 8 (VI) Comparison of Current and Projected Program 
Execution Times for the Processing System of Figure 2. 2. 8(3) 



Number of 
Microprocessors 

8080 

. Microprocessor 
3.5 ys ADD 

1974 

Microprocessor 
1.0 ys-ADD 

1983 

Microprocessor 
100 ns ADD 

PEI 

1 ■ ■ ■ 

678.0 ys 

193.7 ys 

19.37 ys 

PE2 

12 

406.0 

116.0 

11.60 

PE3 

1 

1290.0 

368.6 

36.86. 

PE4 

1 

231.5 

66.1 

6.61 • 


• ' The bottleneck is now PE3. Suppose that PE3 is organized as 2 pro- 
cessing elements in parallel. Each processing element handles only 6 of 
the 12 categories. The program execution time in PE3 is halved, requiring 
645.0 ys in the INTEL 8080 microprocessor. Now PEI is the constraining 
section of the pipeline . Proceeding in this manner we arrive at the pro- 
cessor shown in Figure 2. 2. 8(4). The resulting number of processors and 
their execution times for each of the processing elements are summarized 
in Table 2.2.8(VII). 

The maximum likelihood processor of Figure 2.2. 8(4) requires 5,814 
•Intel 8080 microprocessors for implementation with a projection of 171 
microprocessors in 1983. It also requires the sixteen 16K multiplication 
ROMs required in PE2 and PE3 . Another method for implementing the multi- 
plication is with a peripheral, hardware multiplier as illustrated in 
Figure 2. 2. 8(5). The microprocessor outputs the multiplicand and the multi- 
plier to output registers MCAND and MLIER, respectively, and the multiplier • 
generates the product (_PR0D) which is input to the 'microprocessor. The 
resulting program execution times were determined and found to be 
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PE4 

k 




d^(x) 


Figtire 2. 2. 8 (4) A More Efficient Organization for 
Executing the Maximum Likelihood Algorithm. 


Table 2.2. 8 (VII) Comparison of Current and Projected Program 
Execution Times for the Processing System of Figure 2. 2. 8 (4) 


Number of 
Microprocessors 


8080 1974 1983 

Microprocessor Microprocessor Microprocessor 

3.5 ys ADD 1.0 ys ADD 100 ns ADD 


PEI 

PE2 

PE3 

PE4 


2 

12 

4 

1 


354.0 ys 

406.0 

322.5 

231.5 


101.1 ys 
116.0 

92.1 

66.1 


lO'.l'l ys 
11.60 
9.21 
' 6.61 


' PROD 



Figure 2. 2. 8(5) Microprocessor System Organized About 
an INTEL 8080 Microprocessor and a Hardware Multiplier. 
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PE2: T = 22. 5n^ + 16n + 2 ys 

PE3: T = (29. Sri + 10.5)|- ys. 

The program execution times for PE2 and PE3 when organized with a ROM and a 

hardware multiplier are listed below. These differences are primarily due 

to the faster execution times of the instructions for addressing memory as 

compared to that of the INPUT/OUTPUT instructions. 

PROCESSING ROM HARDWARE ’ 

ELEMENT MULTIPLIER MULTIPLIER 

PE2 406.0 426.0 

PE3 322.5 385.5 

Use of the hardware multiplier affects the previous results very little. 
However, since PE2 was the critical section in the pipeline organization, 
the increased execution time in PE2 from 406 ys to 426 ys will cause the 
number of microprocessors to increase . This is reflected in the’ following 
comparison . 


ROM HARDWARE 

MULTIPLIER MULTIPLIER 


no. 

multiplexed 

total no. no. 

multiplexed 

total' no. 


systems 

microprocessors 

systems 

microprocessors 

INTEL 8080 

306 

5,814 

321 

6,099 

1974 Microprocessor 

88 

1,672 

92 

1,748 

1983 Microprocessor 

9 

171 

10 

190 


As long as the multiplication time is much less than the execution time of 
the INPUT instruction, these results are valid. ^ 


2. 2. 8.2 Microprocessor Organizat-ion Using Table Look-Up (TLU) 


Most of the processing in the TLU algorithm involves calculating the 
new address to use as a .pointer to recover the boimdaries for a particular 
dimension. Little processing is required to compare the data with the 
boundaries. If no microprogramming feature is available, the single micro- 
processor approach is inviting. A typical search loop for the TLU imple- 
mentation with the Intel 8080 is: 

Instructions Execution 

' Cycles ’ 


. INB 

CMP M 2 
OP OUT 3 
INXH 1 


85 


/Jump out of Loop 



CMPM 2 
JM OUT 3 
INXH 1 
ADD M 2 
MOV E,A 1 
INX H ■ 1 
ADC M 2 
MOV L,E 1 
MOV H, A 1 
JMP START 3 


23 


/Jump out of Loop 


In the worst case 90% of the processing time is spent in this compari- 
son loop . Since the processor spends 46 ys in this loop ahds in the worst 
case must execute the loop 48 times (n=4, M=12) for each pixel, the maximum 

pixel rate is about 500 pixel/second. To achieve the baseline format data 
0 

rate of 1.3 x 10 pixels/second, the processing load must be distributed 
over 2600 identical systems. 

If no microprogramming feature is available , the pointer calcula- 
tion must be done using the last eight steps in the programming. The over- 
head due to this processing amounts to 12 machine cycles or about 25 
microseconds per pointer. For the n=4, M=12 machine a data throughput of 
10,000 pixels/second is achieved using 36 microprocessors and 300 K bytes 
of memory. The same data throughput may be achieved by paralleling about 
20 single microprocessor units with about the same total memory require- 
ments. The fully parallel system is faster but requires more memory to- 
handle the worst case boundary size. 


2. 2. 8. 2.1 On-Board Computer Organization for the Table Look-Up Pattern 
Classification 

In this section, we present an on-board computer organization for the 
4-dimensional table look-up algorithm described in Section 1.1. 3. 3. These 
include the general specifications, memory planning for "the table entries, 
and the multi-microprocessor system design. 

We consider n=4 channels per scan line so that each pixel is a 4-di- 
mensional vector X = (Xj^,X2,Xg,Xj^)'", each element requiring 6 bits. If we 
include the sign bit and a parity-check bit, we allow one 8-bit byte per 
measurement Thus, the computer can be either an 8 bit or a 16 bit 

machine. We consider M=12 classes and a data rate of 1 MHz. There are 
six multispectral scan lines. We present the computer architecture per 
scan line; the processors for all scan lines are identical. 
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Eppler’s dynamic memory assignment technique discussed in Section 
1.1. 3. 3 requires a memory of roughly 32K 16 bit memory words. 

Figure 1.1. 3 Cl) suggests that the region in measurement space asso- 
ciated with a particular class can be one of two types. Class 3 is com- 
pletely separate from all other classes and shares a boundary with only the 
Threshold Class ; this type of class will be- referred to as disjoint . The 
regions associated with Classes 1 and 2 share a common boundary. If Class 1 
were not present, the boundary for Class 2 would be extended into the 

region now assigned to Class 1 and vice versa. Classes 1 and 2 will be 
referred to as overlapping. • 

Boundary information is stored in main memory for* only "those measure- 
ment space points actually assigned to the class. For this reason over- 
lapping classes require less storage than disjoint classes having the same 
statistics. Therefore the upper- limit on the memory requirement can be 
derived by considering each class separately -( i . e . , the disjoint case). It 
is shovm in [25] that N^(C), the number of Vcilues x^ can have and still be 
inside the measurement space region assigned to Class C, is given' by 

N^(C) = 2 1k^(C)[^'^^ (2. 2.8-1), ■ 

In this equation Q is a user-specified threshold parameter indicating 'the 
maximum Mahalanobis distance 'from the. mean that a measurement vector can 
have and still' be assigned' to that class. For a four-dimensipnal classifier 
using Q=12'. 0 excludes fewer than 2 percent of samples taken from a Normal 
distribution and K^^C'C) is the one- dimensional covariance matrix (Involving 
terms for only the first channel) for Class C. The value N^(C) represents 
the length of the line which results from projecting Region C onto the x^- 
axis as in Figure 1.1. 3Cl). Similarly N^CC), the .number of (x^.x^) -combina- 
tions inside the' region assigned -to Class C, is given by ‘ 

N2(C) = ttQ[K2(C)|^/2 • (2. 2.8-2) ' 

■ ' 

. In this equation K 2 (C) is that part of the covariance matrix for Class C ’ 
which involves x^ and X 2 terms. The value N 2 (C) represents the area which 
results from projecting Region C onto the (Xj^,X 2 )-plane as in Figure 1.1.3 
(1). Similarly Ng(C), the number of (x^,X 2 ,Xg)-combinations inside the 
region assigned to Class C; is given by 

^3^^^ = (2. 2.8-3) 
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In this equation is that part of the covariance matrix for Class C 

which involves the (x^,X 2 ;,Xg)-subspace, Finally Njj^(C), the number of 
(Xj^ 5 X 25 Xg,x^)-combinations inside the region assigned to Class C, is 
given by 


Nj,(C) = I iK^(C)-| 


1/2 


C2.2.8-4) 


In this equation K^(C) is the full four-dimensional covariance matrix and 
N^(C) represents the volume inside the four-dimensional ellipsoid for a dis- 
joint class. 

Equations (2. 2. 8-1) through (2. 2. 8-4) give the n-dimensional volume 
for each class separately^ It is useful to define the total and the average 
n-dimensional volume for an application according to Eqs. (2. 2. 8-5) and 
(2. 2 . 8 - 6) 5 respectively. (Notations were defined in Section 1 . 1 . 3. 3). 


N 


n 


N 

/ N (C) 
^ ,n 
c=l 

N 


N = 
n 


n 


( 2 . 2 . 8-5) 


( 2 . 2 . 8 - 6 ) 


Nc 

The quantity very important in that it is the number of 

different values pointer can have; see Figure 1.1. 3(3). For each pos- 
sible value of main memory must be provided to store and O 2 - The 

quantities and can each be stored in one 8-bit byte , but 0^ requires 
2 bytes because it can exceed 255. As shown in Figure 1.1. 3(3), the total 

memory required to store all and 0^ values for a complete applica- 

. . ^ ^ ^ ^ 

tion is bytes. Similarly 4 N 2 = bytes are required to 

store all values of L^, and Og. The memory required to store all • 

values of L., and H„ is 2N„ = 2N N_. It can be seen from Figure 1.1. 3(5) 
that 4N^' bytes of main memory must be provided' to store all values of 

and 0^^. The total memory required to store all of the boundary infor- 
mation for a complete four- dimensional Table Look-Up classification is 
given by 

• N_ = 4N + 4N- + 4^f„ + 2^„ = 4N (1 +N. + ?„ + 0.5N.,) (2.2. 8-7) 

nCLzoC±Z d 

The main memory requirement was calculated using statistics from a 
study involving a wide variety of sensors, platforms, and applications. 
Flight Line Cl [28] is a 6.4 km by 1.6 km agricultural test site near 
Purdue University; data was collected on June 28, 1966, using the Michigan 
multispectral scanner at an altitude of 800 meters. 
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An upper-limit for main memory requirement (i.e., assuming disjoint 
classes) was calculated using Eqa. (2. 2. 8-1) through (2. 2. 8-7) and the 
covariance matrix for each of the classes. Table 2. 2. 8 (VIII) gives the 


Table 2,2.8CVIII) Computed Upper-Limit Estimates for Main 
Storage and Table-Making Calculations for the Case Q=12 


• Study Class 

Optimum 

Search 

Sequence 


"2 

“3 


"b 

Storage 
Requi red 
(Bytes) 

Number of 
Points 
Tested to 
Computer 
Tables 


Soybean 

1.6,9,11 

18 

230 

2,530 

46,300 

6,530 

46,300 


Corn 

6,9,1.12. 

17 

193 

2,090. 

32,800 

5,020 

32,800 


Oats 

1,6,9,11 

21 

2^j2 

4,600 

90,500 

10,252 

90,500 

C 

Wheat 1 

6,8,9,12 

2C 

284 

2,160 

25,100 

5,560 

25,100 


Red Clover 

1,9,6,12 

15 

158 

1 ,740 

43,200 

4,172 

43,200 


Alfalfa 

1,6,9,11 

■ 13 

126 

1,300 

37,900 

3,156 

37,900 

o> 

Rye 

1,6,9,12 

16 

163 

1,670 

15.700 

4,056 

15,700 


Bare Soil 

9,12,1,11 

11 ' 

93 

760 

8,150 

1,936 

8,150 

to 

Wheat 1 ! 

1,6,9,11 

19 ^ 

614 

11,950 

232,000 

26,432 

232,000 

Uu 

< 

Total : 5^ __ 

— 

156 

2,103 

28,800 

531,650 ^ 

66,636 

531,850 . 


Average: F 


17 

233 

3,200 

59,072 

7,404 

59,094 


search sequence which minimizes the memory requirement and , the. number 
of 8-bit bytes of main memory required to store the necessary tables. This 
table shows that to store the complete representation of the classification 
regions of the measurement space requires at most 66,636 bjrtes in the case 
of Flight Line Cl. For this case 531,850 different (x, ,x_,x ,x )-combina- 

J. ^ v) T" 

tions must be tested to compute the prestored tables. 

In addition to the disjoint case, in which each class is treated 
separately , computations were performed to determine the actual memory re- 
quirement taking into account the overlap between classes. This was 
accomplished by finding the various volumes using the maximum likelihood 
classifier for assumed normal statistics as in [28]. Table 2.2.8(IX) shows 
that the actual core requirement is 92 percent of the worst-case require- 
ment (see Table 2.2.8(VIII)) for the case of Flight Line Cl. 

The results presented in Table 2. 2. 8 (IX) support the following impor- 
tant conclusions: 

(1) The main memory requirements for a four- channel Table Look-Up 
classifier are easily satisfied. Flight Line Cl can be run on a minicom- 
puter with 32K 16-bit words or 64K 8-bit words. 
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Table 2,2. 8 (IX) Actual Requirements for Main Storage and 
Table-Making Calculations for the Case Q=12 


Study Class 

Optimum 

Search 

Sequence 


•‘a 


"« : 

Storage 

Required 

(Bytes) 

Number of 
Points 
Tested to 
Computer 
Tables 


Soybean 

1,6.9,11 

18 

236 

2,2l»6 

150,096 

5,506 

46,300 

— 

Corn 

6,9,1,12 

17 

178 

1,895 

28,569 

4 , 570 . 

32,800 


Oats 

1,6,9,11 

21 

251 

^,172 

78,290 

9,432 

90,500 

<u 

c 

Wheat 1 

6,8,9,12 

26 

302 

1 2,136 

24,074 

5,584 

25,100 


Red Clover 

1,9,6,12 

15 

15 ^ 

1 1,670 

37,934 

4,012 

43,200 


Alfalfa 

1.6,9,11 

13 

■ 129 

1 1,153 ! 

32,243 

2,874 

37,900 

sz 

O) 

Rye 

1,6,9,12 

16 

161 

; 1,509 

13,805 

3,718 

15,700 


Bare Sol 1 

9,12,1,11 

11 

90 

! 762 

8,150 - 

1,924 

8,150 

<r> 

V/heat It' 

1,6,9,11 

19 

587 

, 10,720 

202,556 

23,860 

232,000 

3 

Total : ' 


156 

2,088 

26,263 

465.717 

6l-,480 

531,850 


Average: N 


17 

232 

j 

2,918 

51,746 

6,831 

59,094 




(2) Use of the Table Look-Up algorithm in five dimensional measurement 
space is probably not practical* For five dimensions the memory require- 
ments are 15045,440 bytes for Flight Line Cl* If more than four channels 
are required, it is more practical to compute the four best linear combina- 
tions of channels and use the result in the four- dimensional Table Look-Up 
algorithm'. 

(3) Using Eqs. (2. 2. 8-1) through (2>2.8-7) gives moderately good esti- 
mates of the main RAM memory requirements without having to actually carry 
out the maximum likelihood computations. 


2. 2. 8.2. 2 Multi- Microprocessor Computer System for On-Board Table Look-Up 
Pattern Classification 


For n = 4 channels and M = 12 pattern classes, we need a multiprocessor 
which allows parallel retrieving and comparing boundary information as 
well as concurrent address pointer calculations. A modified version of 
the table look-up algorithm emphasizing the inherent parallelism was de- 


veloped; the required multiprocessop/multi-memory module computer system 
per each class is presented ’in Figure 2.2. 8(7). The lower boundary infor- 
mation L. (?•), the upper boundary information H.(P,), and the pointers 0.(P.), 
Dt 31 ji 
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are stored in memory module M ^ M ^ and respectively. At each memory 

j- 2. o 

cycle, these are simultaneously fetched. The processors are assigned 


special missions; i.e., performs the comparison Lj(P^) 


P^ perforins 


the comparison ^ Hj(P^) and is responsible for pointer updating = 

O.(P^) + X.. Each measurement in the p.ixel x may be input to all three 

processors at the same time so that concurrent executions are enhanced. 

Since each processor is required to perforin only simple computations 
and logical decisions, we use microprocessor approach with microprogrammed 
control. Table 2.2.8(X) shows the system characteristics in the Intel 8080 
microprocessors . 


X. measurement 
1/0 area 


‘t 



^2 1 

pn 

Processors 
(1 ysec per 
instruction cycle) 

(P. )<x. 

f"' 

x,<H.(P.) 

P.=0.(P.)+x. 

1111 

«1 

ll' . 
*^2 

«3 



L.CPp 

: f 

' 9 

f 


H.(Pi)i 

% 1 
. I 


O.(Pi) 

• 

i . 

1 • 


Parallel 

memory 

modules 


Figure 2. 2, 8 (7) A 3 -Processor and 3-Memory Module Subsystem 
for Table Look-Up Classification (Per Class). 


Table 2.2,8CX) Major System Characteristics of the Multiprocessor . 


System for *t-Channel Table Look-Up 


# of u-processors used 

# of 4Kxl RAM chips required 

memory cycle 

processing speed per each 
pixel 

' # of registers or buffers 
required 

Power required 

Total # of IC chips 


36 Clntel 8080 n-MOS) 

-- 

ifSO n sec (n-MOS RAM’s) 

4 psecs or 250 KC data rate - (I' ljsec per addi- 
tion or comparison) 

36 16-bit address buffers; 72 8-bit data 
buffers; (need 2 8-bit buffers, 1 16-bit 
buffer per p-processor) 

36 watts + 100 watts = 136 watts (1 watt/ 
processor and 0.5 watts/4Kxl RAM) 

36 + 300 + 108 = 450 
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2 . 2 . 8 . 3 Microprocessor Organization for the Expanded Maximum Likelihood 
Algorithm • 

•By expanding Eq. (1.1. 3-5) and merging coefficients we generate the 
following .expression for the maximum likelihood decision function: ' 

dj^Cx) = ^4^4 . ‘5' 

■ + \8’‘2=‘3 * \9’^2’‘4 *k,10='3=‘4 + \,lA 


\,12='2 ^■\,13’'3 *lc,14’'4 * *k.l5 


(2. 2. 8-8) 


Figure 2. 2. 8(8) depicts a system archi-tecture for implementing Eq. (2. 2. 8-8) 
in' which each PE is organized about a microprocessor . PE3 is the same as 
PE4 in fhe previous microprocessor system organizations. PEI and PE2 are 
organized for efficient' execution of Eq. (2. 2. 8-8) rather than the matrix 
format presented in Eq. (l. 1.3-5). PE1-. of PEI generates the output P, 

*. K K 

which consists of the partial products formed from the pattern vector input 
X, and a partial sum. 



Figure 2,2,8C8) A Microprocessor System Organized Similarly to the 
Optimal Pipelined Arithmetic Unit 


2 2 2-2 

”^k “ ^1 *^.2 »^3 »^l’^2’^l‘^3»^l’^4»^2’^3»^2'^4’^3*%’ 


PARTIALj^ A]^,ll'^l \,l-2''2 \,13'"3 '^■*k,14''4 

for k = 1 , 2 , . ; 12 . 


^k input to PE2j^ of PE2 which forms dj^(jc). i® input to PE 3 (k = 1, 

...j’l2) for classification. ' ’ . 

Table 2.'2.8(XI) presents the'program for PE2. The program for PEI is 
‘ listed -in Progress Report No. 8. Each PE in PEI and PE2 has a peripheral 
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Table 2.2.8(XI) 

Program for 

Generating in PE2 

MNEMONIC OPERAND 

BYTES 

COMMENT 


LXI 

HjCOEFFA 

3 

initialize coefficient location 


IN 

• PARTIAL 

2 



MOV 

■TEMP, A 

1 

initialize TEMP with partial sum 


IF 

XISQ 

2 

input X 2 


OUT 

MCAND 

2 



MOV 

AjCOEFFA 

1 

input Aj^ 


OUT 

MLIER 

2 

2 


IN 

PROD 

2 

input product Aj^ x*^l 


ADD 

TEMP 

1 

accvunulate sum 


MOV 

TEMP, A 

1 



INX 

H 

1 




• 


multiply each A, . by' appropriate 

2’ 2 2 



« 

« 

‘96 

partial product, X2 5 Xg , x^ , 


. 

• 

• 


Xi*X2, Xj^.^-Xg, x^'Xg, x^'Xj^ 





• ±or. - j • ^ 2; ...» 9 


IN . 

X3X4 

2 

input product Xg-x^^ 


OUT 

MCAND 

2 



MOV 

AjCOEFFA 

1 

input 


OUT 

MLIER 

2 



IN 

- PROD 

2 

input product A^^ io"^3*^4 


ADD 

TEMP 

1 

accumulate final sum 


OUT 

Dk 

2 

output dj^(x) 

hardware multiplier - Table 2.2.8 (XII ) 

gives the comparison of current and 

projected execution times 

for each processing section in Figure 2. 2. 8(8). 


Table 

2.2,8CXII) Comparison of Current and Projected Execution 



Times. for the Processing System of Figvire 2,2,8C8) 



Number of 

8080 

1974 1983 


Microprocessors Microprocessor 

Microprocessor Microprocc 




3,5ys ADD 

l,0ys ADD 100 ns ADD 

PEI 


12 . 

338.0 

96,57 9.66 

PE2 


12 

342,5 

97.86 9.79 

PE 3 


1 

231,5 

66.1 6.61 


In Table 2.2. 8 (XIII) the number of multiplexed , parallel systems and the 
total number of microprocessors required to process data at the acquisition 
rate of 1*33 ps per pixel for the baseline data format defined in Section 
'1.4 is summarized. 
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Table 2 <2,8 (XIII) Conipai*ison of the Numbei? of Micfoprocessoi's 
Required for the Processing System of Figure 2. 2. 8 (8) 

Nuaiber of Total Number of 

Multiplexed Microprocessors 

Parallel Systems 

INTEL 8080 258 6,450 (258 x 25) 

1974 Microprocessor 74 1,850 

1983 Microprocessor 8 200 


Each PEI, duplicates the generation of the partial products formed from 

K 

the input vector. These products are not functions of the classification 
categories and need be generated only once . Figure - 2 - 2 . 8 ( 9 ) presents a 
much more efficient organization. PEIA generates the partial products P 

T. 2 2 

P = ^1 , X2 , Xg-Xj^- 



Figure 2. 2. 8 (9) A More Efficient Microprocessor Organization 
For Implementing Eq. (2 , 2 . 8-8 ) 


in parallel, PEIB generates the partial sum, PARTIALj^, k = 1, ..., 12. 

PEIB consists of four (4) parallel microprocessor processing elements. Each 
processing element generates three (3) partial sums. For instance, PEIB. 
generates 


PARTIAL = A- - . X,, + + A. . _ 

1 1,11 1 1,15 

PARTIAL^ = + + ^2,15 ^ 

PARTIAL^ = + ^ S,15 

PE 2 and PE 3 are the same as in the Figure 2. 2. 8(8) organization. 
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The program can be split between PEIA and PEIB. However, each pro- 
cessing element in PEIB must execute the segment of code related to the 
partial sums three times, one for each of the classification categories. 
Table 2. 2. 8 (XIV) and Table 2.2. 8 (XV) summarize the results for this system 
architecture. Notice that this organization requires the fewest number of 
microprocessors, but not a spectacular decrease. Eq. (2. 2. 8-8) requires 
twenty- four (24) multiplications while the matrix form 

requires only twenty (20) for n = 4. The efficient use of PEIA provided 
the slight improvement. 


Table 2.2.8CXIV) Comparison of Current and Projected Execution 
.Times for the Processing System of Figure 2. 2. 8(9) 



Number of 
Microprocessors 

8080 

Microprocessor 
3 . 5ps ADD 

1974 • 

Microprocessor 
l.Ops ADD 

1983 

Microprocessor 
100 ns ADD 

PEIA 

1 

200 

57.1 

5,71 

PEIB 

4 

384 

109.7 

10,97 

PE2 

12' 

342 

97,7 

■ 9.77 

PE3 

1 

231 is 

66.1- 

6.61 


Table 2. 2. 8 (XV) Comparison of the Number of Microprocessors Required 
For the Processing System of Figure ’2.2.8(9) • 


INTEL 8080 
1974 Microprocessor 
1983 Microprocessor 
^ ■■ 


Number of Total N\amber of 

Multiplexed • Microprocessors 

Parallel System 

289 5,202 

83 1,494 

8 162 


A summary of all the microprocessor organizations is provided in Table 
2.2.8(XVI). The effort to identify the hardware requirements for these 
various architectures culminated in Tables 2. 2. 8 (XVI I) through 2. 2. 8 (XXI). 
Several considerations were taken into account and are: 

1. In calculating the required RAM memory each processing element in 
a particular processing section has the same amount of MM. How- 
ever, the program store was assumed to be in only one of the 
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Table 2.2,8CXYI) Compaplson of Number of Microprocessors for Different 
Function and Multiplication Implementations 


Hardware Multiplier 

ROM Multiplier System of Figure 2,2.8C5) 

System of Figure 2.2.8C5) Using PE in Fig, 2. 2.8(6) 



no. 

no. 



multiplexed' total no. 

multiplexed 

total no. 


systems microprocessors 

systems 

microprocessors 

INTEL 8080 

306 5,814 

321 

6,099 

1974- 




Microprocessor 

88 1,672 

92 

1,748 

1983- 




Microprocessor 

9 171 

10 

190 


Ca) Eq, Cl. 1.3-5) Implementation 



Hardware Multiplier 

Hardware Multiplier 


System of Figure 2.2.8C8) 

System of Figiare 2; 2. 8 (9) 


no. 

nOt 



multiplexed total no . 

multiplexed. 

total no. 


systems microprocessors 

systems 

microprocessors 

INTEL 8080 

. 258 6,450 

289 

5,202 

1974- 




Microprocessor 

74 1,850 

83 

1,494 

1983- 




Microprocessor 

8 200 

9 

162 


(b) Eq. C2.2.8-8) Implementation 


processing elements. Since each processing element is working in 
synchronization in a processing section, the instruction to be 
executed is output to all by the master element. 

2. Different types of multiplication schemes were utilized to access • 
the impact on system program .execution times . ROM multipliers , 

PLA multipliers, and hardware multipliers were considered. 

calculating power, a power standard was developed by identify- 
ing the major components in the system and associating a typical 
power dissipation with each. In future projections, the power 
dissipation was assumed to remain constant. For instance, al- 
though the projection gave fewer and assumed more powerful micro- 
processors, the power dissipation per yP chip remained 1.0 watt. 

The same was true for RAMs, ROMs, PLAs hardware multipliers, and 
I/O registers. A power standard is provided in the bottom row of 
Table 2.2.8(XVII). 
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Table 2.2.8CXVII) System Hardwate 'Requirements - for the Microprocessor 
Processing System Us,ing ROM Multipliers 

Total ROM 


Number 

Required Multiplier 

Niomber 



Volume 

Micros 

Memory 

Memory 

of I/O 

■ *Power 

Weight’ 

(Cubic Cost 

processors (Bytes ) 

(K Bytes) 

Registers 

CKwatts ) 

(Kgraras) Meters) (K$) 

INTEL 5,814 

78,697 

78,336 . 

97,002 

16y286 

978 

1,96 6,300 

8080 

1974 Micro- 





0.57 1,870 

4.684 

1,672 

22,889 

22,528 

27,896 

281 ■ 

processor 
1983 Micro- 




* 




2,665 

2,304 

2,853 ■ 

0.479 

29 

0,06 186 

processor 






Power INTEL 

INTEL 

INTEL 

2-Tl 




Standard 8080 

8102 

8316 

SN74L95 




1.0 

RAM 1024 

ROM 2048 

38 mwatt 




watt 

words X 1 

words X 8 






bit 150 

bits 175 






mwatt 

mwatt 





* , based on Low-power TTL MSI 

and Silicon Gate n-channel MOS technology. ^ Power 

standard is Provided as bottom line of table. 




Table 2.2. 8 (XVIII) 

Same as 

Table 2'. 2. 8 (XVII) Except the ROM ' Multipliers 


Are Replaced by PLA Multipliers 



Total 

Number 

Required 

Number.,,* 

Number 



Volume 

Micro- 

Memory 

of PLA 

of I/O 

*Power 

Weight 

CCtibic Cost 

processors (Bytes) 

Multipliers Registers 

(Kwatts)(Kgrams ) Meters) (K$) 

INTEL 5,814 

78,697 

4,896 

97,002 

12.285 

845 

1,69 5,854 

8080 




« - 



1974 Micro- 1,672 

22;889 

1,408 

27,896 

3,533 

243 

0.49 1,684 

Processor 
1983 Micro- 171 

Processor 

2,665 

144 
DM 7575 

2,853 - 

0.361 

25 

0.05 172 

• 

550 mwatts 






* Same power standard as Table 2. 2. 8 (XVII) 

** National Semiconductor DM 7575/DM 8575 programmable logic array (PLA). A 
16K X 8-bit memory would be required to provide equivalent function. Power 
Dissipation is 550 mwatts . 


4. Realistic estimates were made concerning 'weight , volume, and cost. 

In each table the number of chips required was calculated. Weight 

was generated on the basis of 250 integrated circuits/kg and 125 

pP circuits/kg. • Integrated circuit and circuit board volume was 

3 

calculated on the basis of 8 cm for integrated circuit chips and 
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Table 2,2,8CXIX) System Hardwai’e Reg^uirements for the Microprocessor 
processing System Using Hgrd^fare Multipliers 

Total 

Number Required Number** Number Volume 

Micro- Memory “of ’ of I/O *Power Weight (Cidjic Cost 
processors (Bytes) Multipliers Registers CKwatts)(Kgrams) Meters) (K$) 


INTEL 8080 

6,099 

82,584 

5il36 

101,757 

14,983 

1070 

2.14 

65911 

1974 Micro- 
processor 

1,748 

23,960 

1,472 

29,164 

4.204 

307 

0.62 

1,967 

1983 Micro- 
processor 

190 

2,968 

160 

3,170 

0.467 

33 

0.07 

215 


* Power standard same as Table 2, 2. 8 (XVII) except for multipliers 


** Iterative array multiplier using 8-2 bit x 4 bit binary parallel multipliers 
such as TI SN74LS261, F 93S43, or AM 2505. Includes multiplicand and multi- 
plier registers. Power dissipation 958 mwatts. 


Table 2. 2. 8 (XX) System Hardware Requirements for the Microprocessor 

Processing System of Figure 2. 2. 8 (8) Using Hardware 

Multipliers and Eq. (2, 2, 8-8) Implementation 
/ 

Total 

Number Required Number Number Volume 

Micro- Memory ' of of I/O *Power Weight (Cubic Cost 


processors (Bytes) Multipliers Registers (Kwatts ) (Kgrams ) Meters) (K$) 


INTEL 8080 6,450 

52,991 6,192 

18,834 

13.160 451 

0.90 

3,936 

1974 Micro- 1,850 

15,455 1,776 

5,402 

3,774 129 

0.26 

1,129 

processor 
1983 Micro- 200 

1,991 192 

584 

0.408 14 

0.03 

122 

processor 






* Power standard same as Table 2. 2. 8 (XVII), Table 

2.2.8(XIX). 



Table 2.2.8(XXI) 

System Hardware Requirements for the Microprocessor 


Processing System of Figure 

2. 2. 8 (9) Using Hardware 



Multipliers and Eq. (2.2. 8-8) Implementation 



Total 

Number 

Required Number 

Number 


Volume 


Micro- 

Memory of 

of I/O 

*Power Weight 

(Cubic 

Cost 

Processors 

(Bytes) Multipliers Registers 

(Kwatts ) (Kgrams ) Meters ) 

(K$) 

INTEL 8080 5,202 

57,098 4,913 

22,253 

8.817 418 

0.84 

3,441 

1974 Micro- 1,494 

16,722 1,411 

6,391 

2.533 120 

0.24 

987 

processor 
1983 Micro- 162 

2,218 153 

693 

0.275 13 

0.03 

107 

processor 







* Power standard same as Table 2. 2, 8 (XVII), Table 2. 2. 8 (XIX). 
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16 cm for yP chips. Quantity costs were considered to be $50/yP5 
$10/RAM, $5/R0M, $4/RE6ISTER, etc. 

2. 2. 8.4 Hardware Organization for the Expanded Maximum Likelihood Algorithm 

In this section we design some processors for computing the maximum 
likelihoods given by Eq. (1.1. 3-5) and the expanded version given by Eq. 

(2. 2. 8-8) that are based on hardware implementations as opposed to software 
implementations using microprocessors. This amounts to using hardware multi- 
pliers and adders in an efficient arrangement for performing the required 
confutations. The first two organizations are -based on Eq. (1.1.3 t 5) and 
the third organization is based on Eq. (2. 2. 8-8). 

2:2. 8. 4.1 Serial Organization 

In a serial computer system, we simply compute each of the 16 terms 
g, (ij) of Eq. (1.1. 3-5 for 1 < 1, j < 4 one at a time. The following serial , 

K 

pipelined- arithmetic unit (SPAU) is designed to serve the purpose. Two 
adders and two multipliers 'are required in this SPAU as sho™ in Figure 
2 . 2 . 8 ( 10 ). 



synchronizing clock 


Figure 2.2.8C10) The Serial Pipelined Arithmetic Unit (SPAU) 

The whole serial computer, consisting of three processing elements- 
pipelined in cascade is shown in Figure 2.2.8(11). RAM stands for random- 
access memory and EAPROM means electrically-alterable programmable read- 
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PE 


] 





Figure 2.2,8(11) The Serial Computer Organization 


Table 2. 2. 8 (XXII) Time and Hardware Requirements for the 


Serial Computer Organization 

No. of Arithmetic Units No. of. Memory for Data Total Processing 

. or p-processors used Registers & Coefficients Time Required 


2 adders, 2 multipliers & 
1 multiplexer for PE^ ; 

1 adder of PE^ and 
1 comparator for PE^ 

• 

9 for P,E^ 

2 for PE2 

3 for PE^ 

8-bIt RAM words 
for X 

252 “8-blt EAPROM 
words 

for 

{ Cj^bj^and 

{a }^2 
^ k-^1 

Let t be the 
pipeline period. 
We need 16 t^ 

per class and 
16-tg X 12 = 

192 tg for 12 

classes 

✓ 

Total 7 devices 

Total 1^1 
registers . 




only memory. 

The hardware and time requirements are listed above in Table 2.2.8 
(XXII). 
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In Table 2.2.8(XXIII) = time delay of register 

t^ = time delay of adder; 

t^ = time delay of multiplier 

t^ = time delay of comparator. 

The system speed t = t + max* (t ,t ,t ). 

s r a m’ c 

2. 2. 8. 4. 2 Parallel Organization 

All of the 192 terms in Eq. (1.1. 3-5) for gj^(i,j) for 1 £ i, j £ 4 and 
1 £ k £ 12 can be calculated in parallel . Also , the summing stage in PE« 
and the comparing stage in PEg of Figure 2.2.8(13)' can be done in log 2 (n^+l) 
~ ■ ~ ^ log 2 ^ ” log2l2 = 4 steps instead of 15 or 11 steps , respec- 

tively. A parallel system is shown in Figure 2.2.8(14). Each PE , for 
1 £ ^ £ 12 in Figure 2.2.8(14) contains n = 16 copies of the SPAU shown in 
Figure 2.2.8(10). All of the for 1 £ k = 12 are identical as shown in 

Figure 2.2.8(l3a). The detailed diagram for PEg in Figure 2.2.8(12) is 
shown in Figure 2, 2. 8 (13b). The memory requirement is the same as in the 
serial system. The parallel system is about 192 times faster than the 

97X + 2477 

serial system unit requires ^ — = 164 times more hardware cost. 


Table 2, 2. 8 (XXIII) Time and Hardware Requirements for the 
Parallel Computer Organization 

No. of computing No. of Total processing time 

devices used registers required 


2 X 16 X 12 = 384 adders 

9x16x12 = 1728 in 

One 

and 384 multipliers in 
computing stage 

computing stage 


16 X 12 « 192 adders 

33x12 « 726 In 

for 12 classes 

used In summing stage 

summing stage 


11 comparators used in 

23 In comparing- 


comparing stage 

stage 


Total 971 devices 

Total 2477 
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. . 

(Fig. 

2 . 2.8 

(13b) 


d,2(x) 


Max 


Figure 2.2.8(12) The Parallel Computer Organization 


2. 2. 8. 4. 3 The Optimal Organization 

All of the coefficients in Eq. (2. 2. 8-8) can be computed in 
advance. , This implies a gain in speed of the on-board processor and a 
reduction in hardware. 

In Figure. 2.2. 8(14) the first pipeline stage PE^ contains m=12 identi 
cal optimal pipelined arithmetic units (OPAU) as in Figure 2.2.8(15). The 
comparing stage PE„ is identical with that shown in Figure 2. 2. 8 (13b). In 







adders 



Figure 2. 2. 8 (13 a) The Parallel Sunmiing Stage for 1 k ^ 12 



Figiire 2. 2, 8 (13b) The Parallel Comparing Stage PE^ 


Table 2. 2. 8 (XXIV), we summarize the complexity associated with such an 
optimal system- 













Table 2,2,8CXXiy) Time and Hardware Requirements for the 
Optimal Computer Organization 


No. of Computing 
dev I ces 

No, of Registers 
used 

Memory Require- 
ment for coeff's 
and data 

Total Process- 
ing Time req'd. 

24xl2«288 
multi pi lers 

Uxl2=168 

adders 

11 comparators 

59 X 12 + 23 « 
731 

latches 

A 8-blt RAM 
words for x^ 

15x12+12=192 
8-blt EAPROH 
words for the 
coefficients 

One unit 
t =t 

s register 

+Max(t , , 

adder, 

^multipl ler, 
^comparator) * 

Total 

A67 devices 

• 

I 

: 


Table 2. 2. 8 (XXV) Comparison of the Three Computer Organizations 


Organization 

Total 

Computing 

Devices 

(a) 

Total 

Processing 

Time* 

(3) 

Total 
Memory 
Requl remen t 

Performance 

Efficiency 

Measure 

(a* 3) 

Serial Computer 
(cheap but 
slow) 

7 

192 

256 8-blt 
words 

7 X 192 = 107A 

Optimal Computer 
(med r urn 
cost and 
fast) 

A67 

•' ‘s 

196 8-blt 
words 

A67 X 1 = A67 

’artliel Computer 
(expensive 
but fast) 

971 


256 8-blt 
words 

971 X 1 = 971 


* t = t^ , + Max. (t t . . . , . t ) 

s latches adder, raultxplier, comparator 
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PE^CFig, Z;2.8a5)) 



2. 2. 8.4*4 Comparison of the Three Computer Organizations 

In Table 2. 2, 8 ( XXV) ^ we summarize the processing time and hardware re- 
quirements for the three computer organizations. A performance efficiency 
measure is defined as the product of the number of time units and the num- 
ber of hardware devices required according to the current level of techno- 
logy. From Table 2.2.8CXXV)j we conclude with the following remarks: 

(1) The serial organization would cost the least in hardware at the 
expense of intolerable slowness. 

(2) The parallel system is very costly because of excessive hardware. 

C3) The memory requirement for storing input data and constant coeffi- 
cients is relatively the same for all the three organizations. The S3m- 
chronizing registers (latches) required in each organization is roughly 
proportional to the required device numbers. 

(4) The optimal organization offers an improvement over the other two 
organizations. Its performance measure is roughly 2 times better than the 
strictly parallel system and 2.3 times better than the serial system. 

2.2.8. 5 A Hardware Table Look-Up Processor 

A hardware implementation of the TLU algorithm using currently avail- 
able TTL and MOS technologies is shown in Figure 2.2.8(16). The timing 
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Figiire 2.2,8(15) The Optimal Pipelined Arithmetic Unit (OPAU). 

generator is shared between parallel units (one for each pattern class) and 
contributes insignificantly to the' overall IC count and power dissipation. 
Two TTL IC's are required for each of the digital comparators; the total 
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dissipation in each of the processors would be approximately one watt. The 
processing delay encountered in the comparator stage amounts to no more 'than 
six gate delays. A twelve-bit ^ full carry , look-ahead adder may be imple- 
mented in three TTL IC’s requiring 1.6 watts and performing a sixteen bit 
addition in eight gate delays or approximately 60 nanoseconds. The memory 
address and data input latches may comprise three TTL IC*s which dissipate 
a total of one watt and introduce an additional 15 nanoseconds processing 
delay* 

The time required to generate and latch the new address is 75 nano- 
seconds; the time required to process each feature dimension is this time 
plus the memory access time. For currently available NMOS RAM’s, the 
memory access time is in the neighborhood of 500 nanoseconds, so the real 
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processing delay in this type of processor is caused by memory access delays.. 
An estimated 4096 words of length 24 bits are required for each of the 
classes the processor must recognize. Using current memory density and 
dissipation figures, the number of memory IC's required is about 24 per 
processor and these IC's dissipate about 12 watts. 

The. total IC count for a twelve-class, six-bit-resolution computer is 

approximately 400, and the total dissipated power is roughly 200 watts. 

The volume of the unit is approximately 1/2 cubic foot. 
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2.3 ON-BOARD PROCESSOR ENVIRONMENTAL EFFECTS 


The components of an on-board processor can be divided into several 
classes. Each class will contain a particular type of component, and each 
class will be affected differently by the space environment. The majority 
of the components will be made up of four classes: (1) the electronically- 

alterable programmable read-only memory (EAPROM)j (2) the random-access 
memory (RAM); (3) the microprocessor chip (yP); and (4) the hardwired logic 
chip. The environmental effects on each class of components is heavily 
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dependent upon rhe spacecrcift orbit and the spacecraft environment. 

2.3.1 Orbit 

For earth-resources missions there are two types of orbits. One is the 
earth- synchronous orbit in which the spacecraft continuously vievfs a given 
area of the earth's surface, and the other is the sun-synchronous orbit *in 
which the spacecraft crosses the equator at the same time of day. 

2. 3.1.1 Earth- Synchronous Orbit 

In order that the spacecraft remain at a fixed point with respect to 
the earth's surface, the spacecraft orbit must be equatorial, circular, and 
at an altitude of 35,870. km. Such an orbit would provide opportunities for 
viewing of short-lived events while they occur, or at the first cloud- free 
opportunity, and/or monitoring time- variable environmental phenomena. 

The. major potentially-damaging effect at synchronous altitude is space 

radiation due to the trapped electrons in the Van Allen belts. The electron 

flux can be enhanced by a magnetic storm which raises the quiet-day flux 

0 

more than an order of magnitude. The quiet-day flux is around 10 electrons 
per square cm per sec. [l]. Of more importance is the dose received by the 
spacecraft during the transfer orbit when the radiation belts are pene- 
trated. The major spacecraft damage is to the solar cells since they are 
on the outer surface. The spacecraft shell provides some shielding for the 
electronics. Metal-oxide silicon field-effect transistors (MOSFET) are more 
susceptible to the effects of radiation than are bi-polar transistors. An 
evaluation of the effects of radiation received during a penetration through 
the Van Allen belts can be made by observing the radiation data taken on 
MOSFETs with various thicknesses of shielding. In previous NASA tests, ten 
orbits of penetrating the Van Allen belts each orbit (four days) produced 
a 30-mV shift in the gate threshold for one-quarter gram per square cm of 
aluminum shielding (one thickness of spacecraft skin). One-gram and two- 
grams-per-square-cm shielded MOSFETs showed no deterioration. For the 
types of components projected at present, there appears to be no problem 
with space radiation at synchronous altitude when normal shielding pre- 
cautions are taken. 

Another area of concern in spacecraft environmental system design is 
the temperature control. At s 3 mchronous orbit altitudes the spacecraft is 
almost always in the direct sunlight except for brief periods when the sun 
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is eclipsed by the earth. The temperature can be controlled fairly accurately 
and maintained eft a value which provides good' safety margins. Based on the 
operating temperature range of the existing components in the four classes 
above, no temperature problems are anticipated at synchronous orbit. 

At the time of writing this report, ATS-6 had experienced a year's 
exposure [23 to the synchronous orbit environment. The spacecraft carries 
two small general-purpose digital computers for use in the altitude-control 
system. These computers show no effect to this environment. The tempera- 
ture of the Service Module has varied from 22°C to 26°C which is well within 
the 5°C to 35°C temperature specification. 

2. 3.1.2 Sun- Synchronous Orbit 

The typical earth-resources mission has been flown at -SOO-km altitude. 

The spacecraft crosses the equator at the same local time each revolution; 
LANDSAT-1 erbssed the equator at approximately 9:40 a.m. each orbit. The 
relationship of the ERTS-orbital plane to the centers of the earth and the 
sun remains constant while the earth rotates beneath the observatory [2]. 

At this altitude the earth is completely covered in 18 days. 

The 900-km earth-resources orbit is just below the proton and the elec- 
tron trapped-radiation belts. The flux from magnetic storms does not pene- 
trate to this altitude. Since considerable experience has already been 
gained for operation of field-effect transistors at this altitude, the 
present indication is that normal chassis covers- provide sufficient shield- 
ing against space radiation for MOSFET devices. 

Currently, NASA is experiencing problems with CMOS devices at higher 
altitudes. While early CMOS devices were sufficiently i mmune to radiation 
effects, the process line was changed and later devices show a greater sus- 
ceptibility to radiation. This is a temporary problem, and based on past 
performance, it is safe to predict that during the 1980 's the CMOS devices 
will once again be sufficiently immune to radiation effects. 

2.3.2 Spacecraft Environment 

During the 1980 ’ s the method of launch of. an earth-resources spacecraft 
will be by the space shuttle. The shuttle can place a 6000-kg payload in a 
620-km sun- synchronous orbit; however, it can only place a 600-kg payload 
into a 900-km sun- synchronous orbit. A small decrease in orbital altitude 
results in a large increase in orbited payload weight. 
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To be cost-effective, the shuttle requires a large number of users and 
a large weight-carrying capability. This can be achieved- at the lower alti- 
tudes. Utilization of a higher weight payload usually results in a lower 
cost for the payload. Redundancy can be incorporated plus additional 
shielding so that lower cost components can be used. The additional weight 
can also be used to control the payload environment so that environmental 
testing can be reduced. Based on present-day technology and the successes 
in flying payloads in both earth-synchronous and sun-synchronous orbits, 
the on-board processor designer should have no trouble in meeting the 
spacecraft environmental specifications for either sun- synchronous or earth- 
s 5 mchronous orbits. 
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3 TECHNOLOGY FORECAST AND ASSESSMENT 


The ability of the on-board processors designed in Section 2 to imple- 
ment the algorithms described in Section 1 in real time for the required 
throughput data rates depends on the components that will be available at 
the time of system design. The lead-time required for design, procurement, 
fabrication, checkout, and launch is about 5 years, so that 1S80-1990 
launches will utilize 1975-1985 technology. Consequently, we require ac- 
curate component and system technology forecasts for the next 10 years. 

Section 3.1 deals with performance measurement criteria. 

Section 3.2 contains a survey of the electronic component technology 
available in 1975. Future improvements in component technology from 1975 
to 1985 are projected. 

In Section 3.3 we review the computer system technology available in 
1975 and we forecast future system technology using both manufacturers' 
estimates and a technology forecasting model. 

Section 3.4 contains a survey of existing satellite on-board computers 
and a discussion of future on-board processor technology. 

In Section 3.5 we develop a forecast feedback system that allows the 
incorporation of the projected component and system technologies into the 
on-board. processor architectures. The results are then used to obtain a 
better estimate of projected performances. 

Finally, Section 3.6 contains a discussion of a component and system 
technology that is in its infancy, but which, with sufficient stimulation, 
could make a significant impact on on-board processors toward the end of the 
1980-1990 decade. 

3.1 PERFORMANCE MEASUREMENT CRITERIA 

/ " ~ 

Many measurement performance criteria have been proposed to evaluate 
hardware structures, algorithmic processes, and computer systems. See, 
e.g., bibliographies on performance evaluation by Buchholz [2] and Miller 
£33. An interesting trend can be discerned from the historical literarure 
on performance measures. 

Skalansky £43 suggested a measure to evaluate the hardware performance 
characteristics of several binary adders. The evaluation criteria was 
defined in terms ‘of the adder's efficiency n = R/I. The value of R was 
taken as the number of bits in each operand. Several formulas for generating 
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values of I were proposed using the gate - normalized addition time and the 
number of two- input AND gates and OR gates. 

In an effort to compare the efficiency of algorithms in processing 
matrices and linear equations. Householder [53 indicates that the number of 
numerical operations plus the number of recordings of intermediate results 
generate an overall measure of the efficiency of the computational process. 
If a benchmark such as time were used as the measure for evaluating each 
system's performance, the affect of the system's instruction set and archi- 
tectural hierarchy would be integrated by the execution of the computational 
task and provide a significant indication of system performance. 

A method for evaluating the computing power of various systems known 'as 
throughput [6-83 is defined in terms of the specific task to be performed. 
The task might be the compilation of a FORTRAN statement or the execution of 
an I/O statement. In this context, throughput is measured in jobs per unit 
time. An interesting aspect of this measure is that the operating system., 
compiler, and data management affect the measure as well as the hardware 
architecture and memory hierarchy. This method of evaluation has been pro- 
posed and used for large, complex digital systems in multiprocessing, multi- 
programming, and real-time environments. 

Other meaningful measures for performance measurement and evaluation 
of complex computing systems are based upon workload characteristics [93, 
utility functions [10 3, hardware monitoring [113, and data base of perfor- 
mance [12,133. 

Anacker and Wang [143 presented a method for performance evaluation of 
computing systems with memory hierarchies in which the degradation of per- 
formance due to system architecture and data flow paths in the CPU and 
memory were taken into consideration by calculating the sum of instructions 
plus data words processed or manipulated per unit time by the system in 
carrying-out a given task. 

Leis [153 introduced a parameter called "setup" and explored its rela- 
tionship to processing power and instruction set in evaluating hardware 
design. A computer architecture was Identified as the set of constraints 
the system imposes on the user in executing his computational task; setup 
is a measure of the amount of data manipulation performed prior to the 
system's execution of an instruction identified as a step in the algorithm 
being executed. For Instance, if a system architecture is organized 
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around a single address instruction format^ setup is required to access one 
of the operands prior to the execution of an arithmetic operation. A system 
organized around a two- or three-address instruction format does not require 
the same amount of setup. Changes in setup are measured in terms of over- 
head. The more setup required in a given application points to a lower per- 
formance . 

Trends are toward the development of performance measures which inte- 
grate all the operating characteristics of the computing system. This study 
incorporates these considerations by taking into account both the hardware 
organization (system architecture) and the software overhead (number of 
microprocessor cycle times) required to optimize and implement the algorithms 
discussed in Section 2. 

3.2 COMPONENT TECHNOLOGY 

3.2.1 1975 Component Technology 

Some 1975 micro- electronic technology families are listed in Table • 

3, 2. 1 ( 1)5 ^ qualitative comparison /of some MOS circuit techniques are 

listed in Table 3.2.1(11). 


TaJ>le 3,2,lCXl Sojne 1975 Component Technologies 


LSI ON-CHIP 



DENSITY 

PONER-DELAY PRODUCT, pj 

SMALLEST' 

FAMILY 

GATES/mm'^ 

15V 

5Y 

IV • 

DELAY, ns 

SCHOTTKY BIPOLAR 

30-40 


5 


2 

CMOS 

30-40 

50 

.5 

- 

10 

STATIC NMOS 

80-120 

50 

5 


20 ■ 

CMOS/SOS 

80-120 

25 

3 

- 

3 

I2l BIPOLAR 

100-120 


5 

1 

10 


The reliability of these components is illustrated by the statistics. 

3 

listed in Table 3.2.1(111). The mean time between failure (MTBF) of 10 - 

10^ hours for a 1000 component system is marginal for an on-board processor. 
Present LSI technology allows up to several thousand gates per chip. 

The maximum power dissipation per chip ranges from 0. 5-2.0 W; very few 
dissipate more than 1 W. 

Two present day memory technologies are illustrated in Table 3.2.1(IV). 
Both are 8^000 bits. The charge coupled device (CCD) has a faster transfer 
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Table 3.2,1CII) Qualitative Comparison of MOS Circuit Techniques 
Circuit 

Technique Speed Power Density Advantages Disadvantages 


Static 

Logic 

0.-2 

MH 

z 

5 mW 

Per 

Gate 

200-300 
Gates 
Per Chip 
Readily 
Obtained 

Runs at any frequency d.c. 
to max. No external clocks 
required. Easiest of PMOS 
types to interface with bi- 
polar. Simplest to use for 
system design. 

Slowest . Highest 
Power, Higher cost 
than 2i> or 45^ 
(larger chip area 
per function). 

2 i> 

Dynamic 

10 Kc- 
5 HH 

z 

Lower 
Power 
than 
Static 
^ Logic 
)on 
Pulse 

Higher 
Density 
than 
Static 
Logic . 

Lower Power. Higher Speed. 
Better Noise Immunity than 
Static Logic, 

Minimum Frequency 
of Operation. Ex- 
ternal clocks must 
be supplied, there- 
fore more system 
noise and complica- 
tion. 


(depending uf 
Frequency anc 
Width . 

4 0 ■ • 

Dsmamic 
Logic 

10 Kc- 

10 MH„ 
z 

Only 
Reac- 
tive 
Power, 
)wer is 
static 
lie 

’owatt 

dissi- 

Higher 
Den sity 
than 2^ j 
Dynamic j 
Logic 

Still Lower Power. Higher 
speed than 2 <f>. Better Noise 
Immunity than Static Logic 
r 2<l> 

Minimum Frequency 
of operation . Ex- 
ternal Clocks re- 
quired, therefore 
more system noise 
and complication. 


so total pc 
less than s 
or 2i> Dynan 
Logic. Mici 
d . c . power 
pation 

Comple- 
mentary 
MOS Logic 

0-25 . 

MH 

z 

Micro- 

watt 

id.c, 

power 

dissi- 

pation 

Lower I 

Density 
than - 1 

Static 1 
Logic , j 

Approx, i 
300 Gates 
Per Chip.] 

Lowest System Power. Highest 
'speed. Good Noise Immunity, 

;S ingle Power Supply. Low 
putput Impedance in both 
directions . Wide Power 

Supply Variations. 

[ 

! 

i 

1 

Largest MOS Chip 
Area* More pro- 
cessing steps re- 
quired. Highest 
MOS die cost. 


Table 3,2,1CIII) 1975 Component Reliability 


Cbi^s passing die sort electrical test 
Chips passing visual inspection of lead bonds 
Components passing all tests after packaging 
Components operable at first system test 
Components operable after 18 hour burn at 
maximum temperature . 

Subsequent failure rate of svirviving components in 
normal service 

Lab or computer room severe environment 
MTBF for 1,000 component system 


100 

98 

80 

78 

76 

.001— ,01 per 1,000 hours 
.01-0.1% per 1,000 hours 
lO^-lO^ Hours 
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Table 3. 2.1 (IV) Comparison of RAM and CCD for Bulk Storage 


Die Size (MILS) 
Access 

Bit Transfer Rate 
Active Power (MW/8K) 
Standby Power (MW/8K) 
Technology 


8K RAM (IBM, 1973) 

145 X 201 (29K MIL^) 
Random to 64 b Blocks 
1.8 X 106 b/s 

22.5 

2.5 

Std. P-SI Gate 


8KCCD (Bell-Northern, 1974) • 

168 X 178 (30K MIL^) 

Random to 256 b Blocks ' 

10^ b/s 
15 
1 

2-Level N-SI Gate 


rate but slower access time because it is a d 5 maTnic shift register device so 
that 5 like a' disc memory ^ readout is delayed until the data comes past the 

readout electronics^ 

\ 

3.2*2 1975-1985 Component Technology 

The equivalent number of gates per IC package (ChipO increased two 
orders of magnitude from 1960 to 1970 (from about 5 to about 500), and 
another order of magnitude between 1970 and 1975 (about 5,000 by the end of 
1975). Three more orders of magnitude remain before the optical diffraction 
limit is reached. Figure 3. 2. 2(1) illustrates this evolution from SSI 
(small scale integration) through MSI (medium scale integration) to LSI 
(large scale integration) [203. All indications are that the number of 
gates/chip will increase by another two orders of magnitude between 1975 
and 1985. 

For some IC’s the number of pins on the package must increase in order 
to effectively utilize the increased number of gates. The pin count per 
package has increased exponentially over the past 15 years. A less than 
exponential increase is predicted for the next 10 years because of problems 
in manufacturing, assembly, seal, and testing. The insertion yield and 

repairability will also tend to slow the growth. A fast 16-bit CPU would 

/ 

require 100-150 pins. This appears feasible as illustrated in Figure 
3.2, 2(2). 

Future memory technology is projected in Table 3. 2. 2(1). For the 7 
year span that is shown it is expected that cell area, povjer/bit, and MTBF/ 
bit will improve by about one order of magnitude. Storage capacity and 
cost increases approach one-and-a-half orders of magnitude. Some experts 
are predicting a 128K bit chip by 1980. All indications suggest between 
one and two orders of magnitude improvement between 1975 and 1985. Power 
requirements are directly proportional to cycle time, but only to the 
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COMPONENTS/CHIP 



1960 1970 1980 


Figure 3. 2. 2(1) Maximum Components/Chip vs Time (Projection) 



" "I960 ' 1970 1980 

Figure 3. 2.2(2) i,SI Package Pin Count Limitations 


-Table 3*2. 2(1) Possible Trends in I^namic MOS RAM 500 ns Cycle Time 


Storage Capacity 

IK* 

4K 

16K 

64K 

Date Introduced 

1972 

1974 

1976 

1979 

Cell Area, Mil^ 

6 

2 

1 

0.5 

Chip Area, In^ 
Power/Bit, pW 

.02 

.03 

.04 

.06 

Active 

300 

100 

20 

4 

Standby • 
MTBF/Bit, Hr. 

4 

lOlO 

1 

0.5 

0.2 

loii 

Component Mfg. Cost/Bit, 2 yr. After Intro. 

0,3 

0.1 

.03 

.01 

Memory System Price/Bit to OEM 

1.2 

0.4 

.12 

.04 


square root, of storage capacity (Power /# bits * access 'time) and are not 
expected to pose any problems • ' 

Present prices of microprocessor components reflect an attempt by sup- 
pliers to recover some of their development costs. The factory cost of a 
microprocessor (or any manufacturable LSI component) is under ten dollars. 
Rapid decreases in the price of microprocessor components are likely over 
the next fev7 years. Memory , input-output devices » and power supplies may 
well be more important factors in the total cost of a system. ' 

The engineering cost of a microprocessor chip is substantial ^ perhaps 
on the order of 0.5 to 1 million dollars. This is a significant barrier to 
the development of custom microprocessors for single users. 

The performance, cost, and reliability of many microcomputer systems 
is determined by memory. Continued rapid decrease in the cost of memory is 
likely. Charge-coupled device memories and bubble memories are not likely 
to challenge the dominance of conventional semiconductor memory techniques. 
Eventually, read-write memory and mask-programmed read-only memory will 
approach the same cost, while electrically programmable read-only memory will 
continue to be more costly. 

It should be possible to build highly reliable microcomputer systems 
through use of well-planned but relatively cheap testing, screening, and 
burn-in procedures. 

Some 1980 LSI costs projections are presented in Table 3.2.2(11) alo: 
with some of the parameters governing costs. 

3.3 SYSTEM TECHNOLOGY 


System technology is progressing at a rate comparable to component 
technology. Some 1973 and projected 1978 minicomputer and microcomputer 
system characteristics are illustrated in Table 3.3(1). 
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Table 3.2.2(11) Projected 1980 LSI Component Prices 


Chip size for 30% yield 
Cost/good Chip 
Number of devices/chip 
Number of gates /chip 
Factory cost of package 
Factory cost/gate 
Selling price/gate 


MOS 

BIPOLAR 

, ^ 2 

2 

1 . 5 cm 

■ 0.6 cm 

$9.00 

$3.60 

10^ 

10 5 

10^ 

10^ 

$28.00 

$17.00 

.03<|: 

.17« 

.06<:- 

.35<|: 


Table 3.3(1) Comparison of Minicomputer /Microcomputer Characteristics 


1973 ■ 1978 


* 

Minicomputer 

Microcomputer 

Minicomputer 

Microcomputer 

Execution time (ps) 

0. 5-2.0 

2.0-25.0 

0.1-2. 0 

0.1-10.0 

Word Length (bits) 

8-32 

4-8 

8-32 

4-16 

Number of Instructions 

100-200 

20-60 

150-250 

100-200 

Technology 

Core 

MOS 

Bipolar /MOS 

Bipolar /MOS 


3.3.1 > Technology Forecasting .Model 

Under the Apollo program El6ll NASA developed a system to forecast space 
. - ' *• * 

vehicle weight .and performance based upon trend forecasting. Because of The 
success of this technique and the excellent documentation available, it was 
adapted to the forecasting of computer performance.’ 

As illustrated in the preceeding sections, component and system para- 
meters are progressing at an exponential rate. Consequently, our model is 
based on the logarithm of the NASA linear maximum likelihood model.- 

The computation and plotting routines of the linear maximum-likelihood 
model were modified to accept the logarithm of the input data. To make a 
projection that would be good to the year 1990, input data over the range 
of 1960 to 1975 was needed. The periodical. Computers and Automation , 
publishes a monthly survey of computers along with their first date of in- 
stallation. Every few years, beginning in 1960, the periodical has published 
a digital computer specification survey. This survey includes cycle time, 
add time, multiply time, etc. Using the monthly survey and the specifica- 
tion survey, data versus .time can be extracted for use in the prediction 
program. 

Figure 3.3(1) resulted, from using input data taken from minicomputer 
surveys for the years 1970 to 1975 to predict computer cycle time.. A 
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Figure 3, 3.1(1) Computer Cycle Time Forecast 


more meaningful result was obtained by taking minicomputer add times ^ for 
the years 1970 to 1975 and combining it with the add times of computers from 
1960 to 1970. ' The resulting prediction is shown in Figure 3.3(2). The 
expected value for the computer add time in 1983 will he 100 ns. The small 
circles and asterisks represent the 80% confidence limits. This curve is 
based upon expected values. Even today there are machines with add times 
of 100 ns» but the average or expected value today is 800 ns. These com- 
puter generated projections are in close agreement to the predictions made 
by experts in the electronics industries discussed in Section 3.2 and 3.3. 

-3.4 OM-BOARD 'PROCESSOR TECHNOLOGY 

NASA is at present in a program for standardizing its on-board com- 
puters. ' The -Goddard space Flight Center has developed the Advanced On- 
Board Processor (AOP) [2l3 which had its first flight on OAO-3 21 August 
1972. There is an on-going program to reduce size and power consumption 
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Figure 3. 3.1(2) Computer Add Time Forecast 


of this computer so that it can be used on smaller spacecraft. The machine 
has an 18-bit word length and a 4-microsecond add time. 

The Marshall Space Flight Center is sponsoring several versions of on- 
board computers: The Hybrid Technology Computer (HTC) by IBM, the Space 

Ultra-reliable Modular Computer (SUMC) by RCA, and the Hughes Aircraft Long- 
Life Fault- Tolerant System. The RCA SUMC is being made of Silicon-on-Sap- 
phire CMOS. 

Technology requirements for the on-board processor differ from computer 
system technology because of space, power, weight, reliability, environment 
and cost constraints. Many of these parameters, e.g. , space qualification 
and check out time and costs, increase as the square of the number of com- 
ponents. Experience has shown that, while the number of transistors or 
gates in a logic system. has increased by orders of magnitude, the number of 
integrated circuit packages has remained essentially constant. To support 
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this theory 5 rhe advice of an expert in digital logic design 5 D. R* Lokerson 
of NASA/GSFC furnished the data shown in Figure 3,4CD* This figm?e shows 
the number of IC packages and the nuinber of transistors or gates for the tele 
metry encoding system for the Interplanetary Monitoring Platform (IMP) space- 
craft over the ten year period 1965-1975 • The IMP telemetry encoder in- 
creased in computational power over the past decade. It used the newest 
technological developments; it was the first spacecraft to use integrated 
circuits and the first spacecraft to use MOS transistors. Note that the 
number of IC^s remained essentially constant while the equivalent number of 
gates increased by two orders of magnitude over the 10 year period. It is 
the opinion of Lokerson that systems with more than 1,000 packages become 
increasingly difficult to test. While the number of gates in the packages 
may continue to rise, testings checkout, and reliability of pin connections 
will dictate a maximum system size of approximately a 1,000 packages per 
system for the next ten years. 


o 



Figure 3.‘4(l) Numbers of IC Packages and Gates for the IMP Spacecraft 
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The constraint of 1,000- chips is based mainly on cost. The technology 
is available for building larger processors; however, the cost to design, 
test, and check out a processor goes up with the number of packages much 
faster' than linearly. We will, of course ,' investigate the sensitivity of 
our results to variations in this data point. 

In summary, we conclude the following for the on-^board processor for 
1975-1985: 

(1) The maximum allowable number of IC packages will stay constant at 

1 , 000 . 

(2) The add time will increase by one order of magnitude. 

C3) The number of gates per IC package will increase by one-to-two 
orders of -magnitude. 

(4) The data throughput capability will increase by two-to- three orders 
of magnitude. (This follows from items 1, 2, and 3.) 

3.5 FORECAST FEEDBACK SYSTEM 

Because of the complex interrelations between on-board processor re- 
quirements, architectures,* etc., a systematic* approach based on a feedback 
model v/as developed. The object of the forecast feedback system is to help 
give the best estimate of on-board processor architecture organization and 
performance. The forecast feedback system block diagram is shown in 
Figure 3.5(1). 
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Figure 3,5(1) Forecast Feedback System Model 
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3.5.1 Input Data 

Data is £ed into the loop in tuo inputs* The p;?iniary input is the on- 
board computer requirements data. These' are the data requirements generated 
in Section 2.1 of this study. To generate the first set of data, typical 
applications were selected in the fields of agriculture, geography, meteoro- 
logy, and oceanography. On-board processor requirements were generated to 
satisfy the spatial and spectral requirements. They are fed into the loop 
by means of the block in Figure 3, SCI), labelled Algorithm Appraisals, Here 
the appropriate algorithm is selected or indicated. In the first iteration 
of the loop there is no feedback to perturb the algorithm appraisal. The 
algorithms are then passed on to the block .labelled Management Decision 
Strategy Trade-Offs. The second set of inputs occur in this block. These 
inputs are the program constraints, the component history data and any sti- 
mulus to industry which may affect the rate with which technological ad- 
vances are made. NASA management, along with the members of the study team, 
made the decisions and trade-offs required to select the appropriate algo-' 
rithms consistent with the program constraints. The algorithms are then 
passed on to the computer architecture designer Cbox labelled Computer 
Architecture and Organization). 

3.5.2 Output Data 

The' output from the feedback system is taken from the box of Figure 
3. 5 Cl) labelled Computer Architecture and Organization. The output is two 
computer designs, one a general-purpose computer and the other a hardwired 
special-purpose computer. Each is based upon using present-day technology 
without benefit of a technology forecast. 

3.5.3 Feedback Loop 

The component parameters required for the implementation of the tvjo 
designs are fed into the technology forecasting program for a projection of • 
the expected performance in the 1980’ s. The parameters were ranked and 
weighted j processor add time was deemed the most meaningful and useful 
'parameter. This information was passed to the Algorithm Appraisal box 
where the selection between the maximum- likelihood and the table look-up 
algorithm occurred. Two new designs resulted, one based upon using a set 
of general-purpose computers and the other on the use of hardwired logic. 
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A second and third iteration around the’ loop produced the final processor ' 
architectures described in. Section 2, 

3.6 OTHER TECHNOLOGIES 

There are some new approaches that show promise for drastically in- 
creasing the computational power of an on-board processor, e.g. , the tse 
computer and the Josephson tunneling circuits. Both would require consi- 
derable stimulation to develop them to their full potential. 

3.6.1 The’ tse Computer 

Image information is two-dimensional. Two-dimensional data processing 
is usually executed as a sequence of serial computations by a word-oriented 
machine. Several parallel-processing architectures [22-24] have been pro- 
posed in which the processor organization is an array of processing elements 
capable of operating on a number of words simultaneously. Most of these 
systems have not reached an operational status because of the prohibitive 
cost involved in their construction. 

Some forms of .optical computing operate at speeds exceeding those of 
digital techniques. Optical computing has generally been limited to optical 
analog devices based upon the mathematical concepts of coherent or Fourier 
optics and holography. Optical analog computing derives power from the fact 
that the computing systems are relatively simple and that parallel pro- 
cessing is a natural property of the lens [25]. The area of optical digital 
devices is not as well-developed. Strong[26,27] has explored the feasibility 
of performing -logical operations on binary images using coherent and non-co- 
herent optical techniques. This research indicates that a two-dimensional 
computer architecture using coherent and non-coherent optical devices has 
attractive image processing capabilities. 

A program at the Goddard Space Flight Center has removed the ground 
rules with which conventional computers are designed and a ‘new design philo- 
sophy and concept has emerged. Schaefer and Strong [28,29], Earth Observa- 
tion Systems Division, NASA/6SFC, have proposed a family of two-dimensional 
logic devices capable of performing simple, parallel logical operations 
simultaneously on one or two binary images. The following summarizes a' 
small part of their research. 


126 



Consider an- image composed .of a 512 x 512 rectangular 
array of picture elements’ in vjhich the gray level pf each element 
is quantized to six bits. One v/ay to visualize the digitized 
image is as six binary image planes, each plane containing 512 
. X 512 bits. The binary image plane or bit plane is a two- 
dimensional binary data array called a "tse" which comes from the 
word for the Chinese writing characters. Just as one tse repre- 
sents many English letters, one binary tse represents many binary 
bits. 

A family of tse . logic devices which utilize electro-optical 
technology to perform simple, parallel logical operations simul- 
taneously on one or two tses has been proposed. Figure 3. 6. 1(1) 
illustrates a tse gate capable of AHDing two binary image planes. 
The interleaver is a passive element which combines corresponding 
positional elements in the image A and the Image B inputs to the 




s LOGICAL 1 


1 I = LOGICAL 0 


Figure 3. 6. 1(1) .A Two tse Input, Digital AND Gate 
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same elemental position at the. interface of the electro-optical 
device.' The electro-optical device is an active integrated 
circuit which converts the optical inputs -into electrical sig- 
nals which are logically ANDed in a conventional manner. The 
output electrical signals are converted "to an optical output by an 
electro-luminescense' process; Since the fan-out of the active circuit 
is one, an interleaver must be used in reverse to duplicate the output 
image and to increase the effective fan-out. Each output from the” ■ 
duplicator must interface to a REfORMATOR which is an active tse buffer 
device to restore the proper optic signal levels. Figure 3. 6.1(2) 
demonstrates how the effective fan-out from the AND gate can he 
increased to four. 
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tse REFORMATOR 


Figure 3. 6. 1(2) Use of DUPLICATORS to Increase' Effective Fan-Out 

Of a tse Device to Four 


Other elementary tse operations are OR, Exclusive-OR , 

NEGATE, and SLIDE operations. These primitive operations and 
the results of performing them on typical images are depicted 
in Figure 3. 6. 1(3). All of these operations, with the exception 
of the SLIDE operation, represent conventional logical operations. 
The SLIDE operation performs an image translation • either up, 
down, right, left, or combinations of these directions. Con- 
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Figur>e 3. 6, 1(^3) An Example of Elementany tse 
Operations on Typical Images 


ceptuallyj this operation is generated by interfacing two fiber 
optic bundles with a physical offset between the bundles. There- 
fore, when a tse is transmitted along the optic data path of one 
bundle, the tse at the output of the second bundle is transmitted 
with the same physical offset. Positions in the' output path 
which are not interfaced to an input because of the offset are 
masked in such a way that they contribute logical zero in the 
output image. 

A tse "contractor" device is a control device to indicate 
the presence of a 1-element in any position of the tse. If 
there are no 1-elements in any position of the binary image, 
the' output of the device is logic-0; otherwise, the output is 
logic-i. This device is different in that the input is a data 
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plane, but the. output is a single logic signal. A device. of 
this type is very. necessary, for implementing conditional image 
operations . 

Tse devices can.be interconnected in much the same manner 
as conventional digital logic gates to implement special purpose 
architectures for executing specific algorithms. As they are 
interconnected to organize more complicated structiores, methods 
for controlling the devices are essential . All ' active tse de- 
vices are assumed to have a one-bit control line, for turning 
them on or off. In the off state the output tse is a zero-tse 
Call elements in the array are logic-0). 

An example of a system organization with tse components is the LOGICAL 
OPERATIONS UNIT of Figure 3. 6. 1(4). This unit is a subsystem in a larger 
tse computer organization CSOj' for executing Strong’s parallel counting 
algorithm C28j. 

Although tse components are still in an early developmental stage, the 
parallelism produced by these devices readily projects a system structure 
which allows for the parallel, concurrent, or simultaneous execution of 
.processing tasks. Future advances in integrated' circuit and optical tech- 
nologies in the next 5-10 years are expected to provide effective computer 
architectures for processing two-dimensional data. 

3.6.2 Josephson Tunneling Devices 

For over twenty years, attempts have been made to utilize the principles 
of superconductivity in the building of computers. If the extremely low 
losses of the superconductive state could be utilized in a logic element, 
then a very low-power computer could be built. In 1956, D. A. Buck [SI] 
demonstrated the "cryotron," which was capable of performing both logic and 
memory functions. This device was not competitive with semiconductor com- 
puting elements in speed; therefore, it has not made an impact in the com- 
puting field. In 1962, B. D. Josephson predicted [323 that supercurrent 
could flow through an insulator sandwiched between two superconductors. 

The predicti’on was later confirmed, and his work spawned a whole host of 
investigations into devices utilizing this principle. The promise of ex- 
tremely fast switching speeds at very low powers is now being realized 
with the implementation of the Josephson Tunneling Logic (JTL). 
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Figure 3. 6. 1(4) 


Orgamza'fcion of the tse Logical Operations Unit 


The first measurements [33j of the speed-power product of this device 
showed an improvement of three orders of magnitude over that of the fastest 
semiconductor logic. This represents a speed capability of a single 
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processor which,' if it were built es an on-board processor, could do the 
multispectral scanner' processing job serially in real time. 

As reported by W. Anacker' C34D, the basic gate is made iip of two super- 
conductors, a and b, with an oxide j c, between them as shown in Figure 3,3.2 
(1). Control lines, d, e,'and f, insulated from the gate, supply the input 



Figure 3,6 .'^Cl-*) Basic Josephsbn Tunneling Gate 

control signals. Each gate is supplied with a gate current, Ig, in and out 
through strip lines, h and i, respectively. For the JTL to be operative, 
this current must be present. With Ig present, and with no control-line 
currents present, maximum current flows through the gate, and the gate vol- 
tage is zero. With the application of a control current in one of the con- 
trol lines and in the direction of the gate cin?rent, the gate current will 
switch to a low value and the voltage will rise to a quantum level. This 
sends a pulse down the output strip lines, j and k, which are terminated in 
their characteristic impedance with a resistance 2R^. The current in one of 
the output strip lines is used to control the next Josephson tunneling gate. 
Depending upon the ratio of the gate line width to the control line width, 
the gain, Vig/Vic can be made greater than unity. This allows the direct 
interconnection of units with no need for amplification of signals. The OR 
operation is formed by a current on any of the three control lines, d, e, 
and f. The AND operation is accomplished by using the opposite-sense cur- 
rents in the control lines. The opposite-sense currents must be present* in 
all three control lines for the devices to switch. Since the gate output 
strip lines are the control for the next gate, the use of one of the lines 
will set up an OR operation, and the use of the other will set up the AND 
operation. A NOT operation' is also available, but it must be performed in 
a time sequence. Once a gate has switched, it can only be brought back to 
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the initial state by removing the gate current. Memory cells have also been 
built using three gates per. cell; The direction of the current in the stor- 
age loop is a measure of the stored bit value. The memory is non-destructive 
readout.’ A switching time of *600' ps* has been achieved [353 ^ which indicates 
access rates' of 1.5 GHz. The write cycle ^ approached 1 ns. 

The performance [33] of the JTL device is depicted in Figure 3. 6,2(2). 
The measured speed-power product of several gates in series showed an average 





flight use, since the space vacuum can .he used to prevent heat losses into 
the supercooled . system . Likewise ,, the space shuttle provides the capability 
of insuring th,at the cryogenic systems can be orbited. at .low cost. 

A rule-of-thumb by . J. K. .Hulm I36l] states that the' improvement using 
cryogenics, to.be worthwhile,' must generally be a factor of at least 1,000. 
Using this rule', there will be no sudden 'rush to the use of Josephson Tunnel- 
ing devices; however, 'the results' of Fi*gure 3. -6.202) were 'ob-fained using 
non-optimum devices. If improvements in JTL technology are made faster than 
in semiconductor technology, then this S 3 rtem may have a' practical applica- 
tion for space flight use . Because of. the handicap of the low operating 
temperature, the devices will probably never enjoy mass production. There- 
fore, there will be substantial economic reasons to continue to use semi- 
conductor logic circuits. 
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4 FEASIBILITY, TRADE-OFF, AND SENSITIVITY ANALYSIS 

This section identifies the pertinent performance parameters, isolates 
the independent and necessary parameters, and relates these parameters to 
the system requirements for each of the user requirements discussed in the 
preceding sections. This will allow us to determine the feasibility of on- 
board processing for each user type in the 1980-1990 time frame and to per- 
form a tradeoff analysis to determine the sensitivity of our results to each 
of the important system parameters. 

The significant parameters related to the performance of the on-board 
processor are tabulated in Table 4(1). 


Table 4(1) 


On-Board Processor Performance Parameters 



data bit rate (bits/sec) 

number of spectral bands (channels) 

number of bits per resolution element per spectral band (bits) 
swath width (meters ) 

resolution along a scan line (meters) 

satellite ground track velocity (meter s/sec) 
resolution along scan path (meters) 

number of classes 
system cost (dollars) 
system power (watts) 
system weight (kilograms) 
system volume (meters^) 
time (years ) 
number of components 

number of additions 

number of multiplications 

pixel rate (resolution elements/sec) 
system constants (i=0,l) 

system constants (i=0,2,3,4,5) 

System Complexity function (i=l,2,3,4) 


From Section 1,1. 2. 4, the incoming data bit rate is: 

R - S /R, X v/R xnxb = rxnxb 
w L p 

where 

r = S /R-x v/R 
w L p 


(4-1) 


(4-2) 
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The pixel rate r^ the number of ground resolution elements per second ob- 
served by the MSS 5 has been determined for each user requirements One could 
choose any three independent parameters from the relation in equation (4-1). 

chose the set n^ and b to specify system performance. The number of 
classes M is also an independent parameter. 

I 

since each of the processors investigated achieve their speed by dis- 
tributing the processing load between many similar sub-processors, the cost, 
power, weight, and volume are all essentially proportional to the number of 
units required to satisfy the design. • For example, a doubling of the data 
rate R requires doubling the number of parallel processors which corresponds ' 
to a doubling of the power, weight, etc. We now define a system complexity 
function for each of the four processors i = 1,2, 3, 4 that relates the 
effects of the various parameter requirements. By proper choice of the 
system scale constants the functions are synonymous with "cost", 
"power",' "weight", "volume", or "number of components". 

4.1 ' COMPLEXITY FUNCTION DEPENDENCE ON M, n, r, b, and R . 

In this section we define a complexity-function for each of the four pro- 
cessor implementations .described'- in Section 2 in terms of the parameters de- 
fined in Table 4 Cl). The . dependence on the time T is discussed in Section 4.2 

■'4.1.1 ' ' The' Microprocessor Maximum Likelihood Method (yPML) 

The circuit complexity required for both the matrix equation (Sectiori 
2. 2. 8.1) and the reduced form equation (Section 2. 2. 8. 3) is determined by 
the number of multiplication operations. For- the classification of a single 
resolution element the minimum number of multiplications and additions in 
terms of the number of channels n and the number of classes M is given in 
Table 4. 1.1(1). ' ' 

Mathematical Operations for Classification 

Reduced Form Matrix Form 

M(3n + n^) /2 M(n^*+ n + 2) 

(n /2 + n/2) (M+1) + Mn M(n^ + n +1) 

Regardless of how 'these operations are distributed throughout the pro- 
cessing system, for a given data rate the system complexity is proportional 


Table 4. 1.1(1) 
# Operations 

N.’ 


N, 


M 


138 



"to the processing time per pixel. Table 4.1.l(I) shows that for either 
formula the processing time per pixel is proportional to M(n + n) since 
multiplication time greatly exceeds addition time. A higher data rate can 
be achieved by increasing the number of processing elements over which the 
processing is distributed so that system complexity increases proportionally 

to the pixel rate. 

The sensitivity to the number of bits per pixel depends on the multipli- 
cation method used in the processor. The example architectures in Section 
2.2 use the Intel 8080 microprocessor unit but other manufacturers produce 
bitwise expandable units. For these units the complexity of each processing 
element is essentially proportional to the number of bits processed in 
parallel by the microprocessor array. If the multiplications are done by 
ROM table- looR- up the amount of ROM required varies drastically with the 
number of bits in the words to be multiplied, e.g. , to generate a (b+2)-bit 
product from two b-bit inputs requires 

ROM bits = (b+2) 2^^ (4. 1.1-1) 

For b > 6 the processor complexity is dominated by the ROM requirements, and 
system complexity increases exponentially with the number of bits. On the 
other hand, if a PLA or hardware multiplier is used, the circuit complexity 
is proportional to the number of bits (at least for b ^ 16); we conclude’ that 
the system complexity is proportional to b. 

The resulting system complexity function for the microprocessor-based 
maximum-likelihood classifier (]iPML) is 

?! = kj^ M(n^ + n) • r • b = kj^ M(n + 1) R (4. 1.1-2) 

where is a system normalizing constant to be determined- 

4-1.2 Hardware Maximum Likelihood Method (HML) 

The primary difference between the hardware and microprocessor methods 
is that the multiplications and additions that were handled by software in 
the microprocessor method are here replaced one-for^one by hardware units. 

The relationship of total system complexity to the parameters M, n^ b, and 
r remain, essentially the same- However, the optimum pipeline processing 
unit has a data rate which exceeds that required by many users and system 
complexity can be reduced by time-sharing some of the processor arrays. 

Time- sharing requires additional circuitry to multiplex the shared pro- 
cessors and this increases the complexity function as the inverse of the 
product of the pixel rate r and the number of channels n. The additional 
hardware complexity for time sharing is proportional to the word length b. 
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The system complexity fiinction for this method is, therefore, 

?2 = CIC 2 M(n+1)R'+ kg b^/R] >(4.1. 2-1) 

For all systems of interest in this study kgMCntl) R » kg b /R, i.e., 

for any realistic data rate, the reduction in complexity by time-sharing pro- 
cessing subsystems swamps the additional multiplexing complexity. Conse- 
quently, the complexity function for the hardware maximum likelihood pro-' 
cessors is given by 

= k 2 M(n+1) R (4. 1.2-2) 

4.1,3 Microprocessor Table Look-Up (yPTLU) 

Relative to the previous methods, the table look-up algorithm trades 
arithmetic hardware for memory storage. The system complexity depends 
almost entirely upon the memory required to store the feature space decision 
boundaries which depends largely upon the training sample statistics' and' the 
classification algorithm used. For most situations, the bulk of the memory 
is used to make comparisons' in the last dimension searched. Boundaries must 
be stored for each point in the projection of the decision volume on to 
the space spanned by the previously searched dimensions. Examination of 
equations (1), (2), (3), and (4) of Section 2. 2. 8. 2.1 indicates, at least for 
uncorrelated spectral measurements, that the memory requirement increases 
exponentially with the number of dimensions (channels). The empirical re- 
sults of Table 2.2.8(XIII) for a typical agricultural classification problem 
may 'be used in conjunction with Eq. (2. 2. 8-7) to tabulate the, storage re- 
quirements versus number of dimensions used, as illustrated in Table 4.1. 3(1). 
F(n) is a simple mathematical function that is a reasonable model of the ' 

number of memory bytes required as a function of n. 

' i 

Table 4. 1.3 (I) Table Look-Up 'Memory Requirements 
No, of Spectral Bands (n) 12345 

Memory Requirements (Bytes) 18 342 4824 61,000 1,050,000 

F(n) = 1.38 (15)“ ' 21 311 4666 70,000 1,050,000 

An increase in the number of channels n also causes a proportional in- 
crease in execution time since the additional dimensions must also be 
searched. To maintain the same pixel rate the system complexity is in- 
creased to handle the additional processing load. 

If ' the number of bits per channel is increased while keeping the 
system dynamic range constant, the number of pointers which must be stored 
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is increased. The new memory requirements may be expressed in terms of the 
old memory requirements by the relation 

= 4 H [1+2^^^ + „ + (0.5)2^^^ (4.1. 3-1) 

B new Q . 1 n-2 n-1 

•Equation (4. 1.3-1) is a modified version of Eq, (2. 2. 8-7), N^, ..., 
are for the original system, is for the bit-increased system, n is^,, , 

the number of channels, and Vb is the word-size change in bits. In any . . 
practical system most of the memory is used to store the last dimension 
boundaries, so that the last term in equation (4. 1.3-1) predominates. 

Equation (4.1. 3-1) reflects only the increase in memory words (bytes)' 
required to store the pointers and boundaries. The word length in fact must 
also change, perhaps even past the byte level. Pointer length must be 
changed to handle the added memory; increasing the word length by Vb will 
handle both the boundary and pointer requirements. 

Increases in the pixel rate are handled by multiplexing identical sys- 
tems so that the total system complexity is directly .proportional to the 
pixel rate r. Changing the number of classes M changes the total execution' 
time proportionately regardless of the distribution of the processing load. 
The number of parallel systems must therefore be increased proportionately 
to maintain the same pixel throughput rate. 

Taking all of the above factors into account the system complexity 
function for microprocessor table look-up (liPTL) is given by 

Pg = kg M r b n = kg M R X^”^ (4.1. 3-2) 

where kg and X^^ are system constants. 

4.1.4 Hardware Table Look-Up (HTLU) 

The system complexity for the microprocessor system is equally valid 
for the hardware approach except for the normalizing constant , i . e . , the 
complexity function .for hardware table look-up (HTLU) is given by 

= kj^ M R X^’^^ (4. 1.4-1) 

4.2 COMPLEXITY FUNCTION DEPENDENCE ON TIME 

It was shown in Section 3 that computer cycle times, memory access ‘ 
times, and logic propagation delays in newly available equipment decrease 
exponentially with time. To project the complexity functions into the 
future, therefore, each complexity function requires a multiplicative 
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factor of the forrn 

f(t) = A (4.2-1) 

o 

where is the normalized yearly decrease in processing delays and T is the 
number of years into the future. It is implicitly assumed that X^ is approach- 
independent, i.e., the rate of these decreases will occur equally for each of 
the four processor architectiires . This assumption is reasonable since each 
of the four processor architectures utilize the same semiconductor technology. 

We must also ta"ke into account the fact that the components increase in 
complexity and computational power with time. The system designer deals with 
these components on a macroscopic level and the design complexity which he 
must face is essentially independent of the microscopic complexity; it de-' 
pends only on the total number of components used in the design and not on 
whether they are SSI, MSI, or LSI chips. Figure 3.4(1) illustrates the 
evolution of the data processing systems for the IMP series of spacecraft 
over the 10-year period 1964-1974. This is a spacecraft data handling and ' 
telemetry system, hut the design techniques are comparable to our on-board 
processor designs and, indeed, the later systems utilized a computer-like 
processor. All designs were near the feasibility limits for a reliable 
spacecraft system at the time of design. Note that the number of components 
stayed nearly constant between 750 and 1000, while the gate count increased 
exponentially (two orders of magnitude in less than 10 years). Figure 3.2.2 
(1) gives a semiconductor manufacturer’s projection of the number of com- 
ponents on an LSI chip for 1975-1985. The projection indicates a tenfold 
increase over the next ten years. 

If the per IC power dissipation, cost, weight, design effort, etc. 
remain fairly constant over this period as they have over the last 10 years, 
then the complexity functions P should reflect the savings due to increases 
in the computational power of the components, i.e., the increase in the 
equivalent number of gates per IC. ■ This requires an additional exponential 
term in each of the system complexity functions. Alternatively, the value 
of X^ in Eq. (4.2-1) can be modified to accommodate this additional factor. 

The resulting complexity functions are summarized in Table 4.2(1). 

4.3 EVALUATION OF SYSTEM CONSTANTS 

In Sections 4.1 and 4.2 we determined the manner in which the complexity 
functions for each of the four processor architectures depend on the system 
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Table 4.2(1) On-Board Processor Complexity Functions 


Processor 

Microprocessor Maximum Likelihood (yPML) 
Hardware Maximum Likelihood (HML) 
Microprocessor Table Look-Up (yPTLU) 
Hardware Table Look- Up CHTLU) 


Complexity Function 
P, = k, M(n+1) R X 

1 O rp 

P^ = M(n+1) R X^~ 

^3 = ^3 “ h’t 

P, = k, HRX^" V 


parameters M, n, b, ‘r, R, and T. The remaining problem is to determine the 
system constants X^, X^, k^, k^s and k^^. 

4.3.1 Evaluation .of X 

o 

In Section 4 . 2' we described two factors that contribute to the expohen- 
-T 

tial term X^ ; (1) the exponential increase in processor speed due to de- 
creases in computer cycle times, memory access times, and propagation delays 
and (2) the exponential increase in the number of gates per IC. 

From the technology forecast results of Section 3 we conclude that com- 
puter cycle times memory access times, and propagation delays will decrease 
by one order of magnitude from 1974 to 1984. We further conclude that the 
equivalent number of gates per IC chip will increase between one and two 
orders of magnitude during the same period; we choose to use the more conser 
valiive nuni)er. Consequently 5 we choose to allow two order of magnitude 
improvement in system complexity (one order of magnitude for speed and one 
for the number of gates per i.e., we set 

X = .01 (4. 3. 1-1) 

o 

so that 

X = 1.58 (4. 3. 1-2) 

o 


4.3.2 Evaluation of.. Xi 


X^ can be accurately evaluated from the empirical results of Table 
2. 2. 8 (XIII) and Table 4.1. 3(1). That is, for table look-up crop classifica- 
tion for agricultural data with b - 6 bit data words we can equate the X^"^ 
term in Eqs. (4.3. 3-2) 'and (4.1. 4-1) with the F(n) function of Table 4.1.3- 
(I) to obtain 

X^”^ = (isf ' (4. 3. 2-1) 
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The scale factor 1.38 in F(n) can be absorbed into the system constants 
and Setting b = 6 in Eq. (.4, 3. 2-1) we find 

= 1.6 ( 4 . 3 . 2 - 2 ) 

Strictly speaking, this value is valid for b = 6 and is an average 
value for agricultural data. However, almost all users surveyed in the user 
requirement survey summarized in Section 1 indicated b = 6 bit data words 
and the few exceptional cases did not deviate significantly from b = 6. 

(The maximum mentioned was b - 7.) Empirical studies of the table look-up 
algorithm for other than agricultural data are not available. However, the 
decision boundaries for other data have the same general characteristics as 
for agricultural data, i.e. , they can -be modeled reasonably well by second 
order pol3momials as attested to by the fact that the Gaussian maximum like- 
lihood classifier works well for other than agricultural data. Consequently, 
the memory requirements for other than agricultural data are not expected to 
vary significantly from the average for agricultural data. We conclude that 
= 1.6 is a reasonable value for the purposes of this study. 

4.3.3. Evaluation of k^, k^, k^^ 

The complexity functions are now defined except for the normalizing 
constanrs k^. These normalizing constants can be obtained from the data 
points resulting from the detailed processor designs described in Section 
2.2.8. In that section we determined the number of IC packages, the power 
requirements, the weight, the volume-, and the cost of each of the four pro- 
cessor architectures for the base line system using 1975 technology. These 
results along with the definition of the baseline system are summarized in 
Table 4. 3. 3(1) 

Table 4. 3. 3 (I) Complexity Functions for the Baseline System 

• •. Described in Section 1.4; i.e., h=4 Spectral Bands 
M=12 Classes, b=6 Bits, R=31.2M Bits/Sec and 
T=0 (1975 Technology) 


Processor 

• ■ Complexity 
, Function 

#IC's 

(K) 

Power 

(Kw) 

Volume 

(m3) 

Weight 

(Kg) 

Cost 

(K$) 

UPML 


30.0 

3.5 

.50 

250 

1700 

'HML 

^2 

.6 

.3 

.03 

15 

15 

VlPTLU 


90.0 

45.0 

.75 • 

340 

2250 

HTLU 

"4 

1.2 

.6 

.10 

36 

32 
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By setting = l.SS, = 1.6,- n = 4, M = 12, b = 6j R = 31.2 x 10^ 
and T = 0 in each, of the eq^uations listed in Table 4. 2 Cl) and equating to 
the values of in Table 4. 3. 3 CD we determine the scale factors k^, kg 
and kj^' for each of the' performance measure listed in Table 4.3.3(,I). The re 
suiting scale factors for each performance measure are presented in Table 
4.3.3CII). 


Table 4.3.3(11) Scale Factors for the Complexity Functions 
of Table 4 .2(1) for Each Performance Measure. 


Scale Factor 



#IC’s 

Power (w) 

Volume (m^) 

Weight (kg) 

Cost 

1.6x10"^ 

l.Sxlo”^ 

2.67xl0"^® 

1.34xl0“’ 

-4 

9.08x10 

3.21x10“’ 

1.6x10"’ 

1.6x10"^^ 

8.01x10*"^ 

8.01xl0“® 

3.03x10“^ 

1.52x10"^ 

-14 

2.53x10 

1.15x10“^^ 

7.59x10"® 

4.05x10"^^ 

2.02x10“^^ 

3.37xl0“^^ 

-12 

1.21x10 

' 1.08x10"® 


4.4 SENSITIVITY ANALYSIS 

The sensitivity of the complexity functions to the various sytem para- 
meters is illustrated by the parametric curves in Figures 4.4(1), 4.4(2), 
and 4.4(3). These curves were obtained by setting all system parameters to 

their baseline values n = 4 spectral bands, M = 12 classes, b = 6 bits, R = 
0 

31.2 x 10 megabits, and T = 0 except for the parameters listed in the 
figures which were then varied from their baseline values. Consequently', 
these curves illustrate the deviations in the complexity functions as the 
system parameters vary from their baseline values. 

-T 

All of complexity functions involve the terms R A^ , i.e. , all four 
coropleKity functions increase linearly with the bit rate R and decrease 
exponentially with the time T. Figure 4.4(1) illustrates the behavior of’ 
the complexity functions as a function of T as R is varied in multiples of 
two from 7.5 megabifs/sec to 120 megabits/sec. ‘ The normalized complexity 
factor, valid for all four processor architectures, is plotted on a linear 
scale . 

The ML and TLU processors differ in their dependence on the number of 
spectral bands n and the data word length b. Recall that R = r • n * b so 
that the ML processor complexity functions increase as bn(n+l). The nor- 
malized ML complexity function is presented in a linear scale in Figure 
4.4(2)- as a function of n between 1 and 20 for b = 3, 4, ..., 8. Doubling 
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NORMALIZED COMPLEXITY FACTOR 


b doubles the complexity function whereas the dependence on n is essentially 

2 

as n . 

The normalized TLU complexity function is plotted on a logarithmic 
scale in Figure 4.4-(3) as a funetion of h and n, Except for very small 



T > ^ 


0 1* 2 3 4 5 6 7 8 9 

T YEARS 

Figure 4.4(1) Complexity Function Sensitivity to R and T (all Processors) 
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NORMALIZED COMPLEXITY FACTOR 


Figure 4.4(2) Complexity Function Sensitivity to b and ri (ML Processors) 


values of b and n the dependence is essentially exponential in both b and n. 
FEASIBILITY ANALYSIS 


The feasibility of on-board processing for the applications areas de- 
fined in Section 1 depends on the particular set of system parameters (user 


I I I I • I I * I I I , i 1 

0 2 4 6 8 10 12 14 16 18 20 

n CHANNELS 

/ 

Figure 4.4(3) Complexity Function Sensitivity to b and n (TLU Processor) 

requirements) n, M, b, R, the processor, the time T, and the definition of 
"feasible”.' 

4.4.1 What is Feasible? 

For a particular processor to be "feasible” at a particular point in 
time requires that it meet certain constraints on performance, complexity. 
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volume 5 weight, power, cost, reliability, environment, etc. Each- of the four^ 
system architectures meets the performance constraint since each was designed 
to accomplish the required task. All four processors use standard integrated 
circuit technology and meet the data throughput rates by adding more compon- 
ents (IC^s) in parallel • The volume, weight, and power dissipation of inte- 
grated circuits can be kept within limits by simply* keeping the number of 
integrated circuits within limits. The radiation, temperature, and other 
environmental constraints can be met by each processor as discussed in* 

Section 2, The limiting factors are cost and reliability vrhich can also be 
kept within bounds by imposing a constraint on the number of components. 
Consequently, we conclude that on-board processing using a particular pro- 
cessor is feasible provided we constrain the number of IC^s in the processor 
to a reasonable nuinber. 

4.4.2 Feasibility Curves 

In this section we plot the number of IC’s as a function of the time T 
for each of the four processors for each user application described in 
Section 1.1. Each figure contains two graphs. The first graph is for the 
maximum data-rate requirements listed in Table 1.1. 2 (I) and the second graph 
is for the minimum requirements listed in the same Table. Each curve is 
labeled \^th the user application symbol described in Section 1.1 and Table 
1.1.2(11). The time scale shows launch date and allows for a five year lag 
between conception and launch. Other time lags can be handled simply by 
adding or subtracting the appropriate number of years from the time axis. 

The parts cost of the on-board processor increases linearly as the 
number of IC^s. The costs associated with check out increase as the square 
of the number of IC’s. Limiting the number of IC’s in the on-board pro 7 
cessor to about 1000 appears to satisfy all constraints, i.e., cost is 
reasonable relative to the total system, cost (launch, sensors, telemetry, 
etc*), reliability is pushing the limits of present day technology as dis- 
cussed in Section 3, while volume, weight, and power dissipation do not 
appear to present serious difficulties. This is discussed in- greater detail 
in Section 3.4. A horizontal line is drawn on each graph at 500 IC^s, 1,000 
IC^-s (which corresponds to our suggestion for a reasonable (feasible) num- ' • 
ber of IC’s), and 2,000 IC’s. 

Some applications are feasible in^l980, some are not feasible by 1990, 
and some become feasible at some time between 1980' and 1990, depending on 
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the type of processor. For example, some of the applications requiring many 
spectral bands Ce.g., oceanography) are never feasible using the TLU pro- 
cessor since the complexity Cnumber of IC.'s) goes' up exponentially with the 
number of spectral bands . • The TLU processors can handle 4 or 5 spectral 
bands at the most." 

A few of the curves indicate’ a fraction of an IC. The reason is that 
an IC is considered equivalent to a certain number of gates that increases 
an order of magnitude over the ten year period. A fraction of an IC (log # 
IC*s 0) means that the on-board processing task can be accomplished with 
fewer gates than available in the "average” IC of that year. 

If we accept the 1000 IC definition of feasibility, then we can list ■ 
the year in which each user application first becomes feasible. Table 4.4.3 
(I) indicates the first year that each processor can accomplish each of the 
applications listed in Table 1.1. 2(1) for the worst case (maximum require- 
ments) situation. The S 3 rmbol N indicates not on or before 1990. All but 
one of the maximum requirements can be met on or before 1986, and about half 
can be met in 1980 by the HML processor. 

Table 4,4.3(11)- is similar to Table ‘^.4.3(1) except that the minimum 
requirements fpr each application were used. The HML processor can meet all 
but one requirement in 1980 the pPML processor can meet all but one require- 
ment by 1990. The HTLU processor can meet most requirements in 1980, but a 
few cannot be met by 1990 due to the number of spectral bands. 

Tables 4.4.3(111) and 4.4.3(IV) are similar to Tables 4.4. 3(1) and 
4.4.3(11) except that 500 IC’s was used as the definition of feasible. 
Similarly, Tables 4.4.3(V) and 4.4.3(VI) are for the case of 2000 IC*s. A 
factor of 1.58 in the number of IC’s corresponds to one year in the date the 
processor becomes feasible. Multiplying the number of IC's by 1.58 makes the 
processor feasible one year earlier. Or, stated the other way, waiting one 
year means the processor can be designed with 1/1.58 = .67 as many IC's. 
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Figure 4*4 ,2(1)3 Number of IC’s vs. Launch Date to Implement the Micropro-* 

cessor Maximum Likelihopd (yPML) Processor for the Minimum 
Data Rate Requirement in Table 1.1. 2 (I), Applications: 
Agricultiire, Forestry and Geology. 
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Figure 4. 4. 2 (2) A Nun)ber of IC’s vs. Launch Date to Implement the Micropro-- 
cessor Maximum Likelihood CyPML) Processor for the Maximuin 
Data Rate Requirement in Table. 1,1,2 (T); .Applications: 
Coastal-r Zone « Studies and Global Oceanography, 
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Figure 4.4.2(2)B Number of IC*s vs. Launch Date to Implement the Micropro- 
cessor Maximum Likelihood (yPML) Processor for the Minimum 
Data Rate Requirement in .Table l'.l.'2(I). Applications: 
Goastcil-rZone Studies and Global Oceanography". 
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Figure 4.4. 2 (3) A Number of IC‘s vs. Launch Date to Implement the Micropro- 
cessor Maximum Likelihood (]jPML) Processor for the Maximum 
Data Rate Requirement in Table- 1.1. 2(1-) : ■ Applications: 
Geography, Hydrology and Meteorology. 
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Figure 4.M-.2(3)B Number of IC’s vs. Launch Date to Implement the Hicropro-' 

cessor Maximum Likelihood (yPML) Processor for the Minimum 
Data Rate Requirement in Table 1.1.2(I). Applications: 
Geography j Hydrology and Meteorology. 
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Figure 4*4.2(4)A Number of IC’s vs. Launch Date to Implement the Hardware 

Maximum Likelihood (HML) Processor for the Maximum Data 
Rate Requirement in Table 1.142(1 ) . Applications : Agri^ 

culture. Forestry and Geology. 
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Figure 4. 4. 2 (5) A Number of IC!s vs. Launch Date to Implement the Hardware 
Maximum Likelihood CHML) Processor for the Maximum Data 
Rate’ Requirement • in Table 1.1. 2(1). Applications: Coastal 

Zone Studies and Global Oceanography. 
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Figiare 4.4.2(5)B Number of IC’s vs. Launch Date. to Implement the Hardware 

Maximum Likelihood CHML) Processor. for the Minimum Data 
Rate Requirement- in Table 1.1.2(-I)r Applications: Coas- 

:aJ.“Zone Studies and Global Oceanography. 
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Figixre 4,4.2C6)A Humber of IC^s ys; Launch Date to Implement the Hardware 
Maximum Likelihood , CHML ) Processor for the Maximum Data 
, Rate * Requirement in Table 1.1. 2(1): Applications: Geo- 

graphy, Hydrology, and' Meteorology. 
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Figure 4.4.2(6 )b; Number of IC's ys. Launch Date to Implement the Hardware 
Maximum Likelihood . CHML ) Processor for the MihiTfiura Data 
Rate Requirement- in Table 1.1.201). Applications: Geo- 

graphy j Hydrology and Meteorology. 
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Figure 4.4.2(7)A Number of IC’s vs. Launch Date to Implement the Micropro- 
, cessor Table Look-Up (|iPTLU) Processor for the Maximum 
Data Rate Requirement in Table 1.1. '2 Cl),- Applications: 
Agricultures' Forestry: and Geology. 
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rignre Number of IC’s ys. Launch Date. to Implement the Micropro- 

cessor Table Look-Up. (liPTLU) Processor . for the Maximum * 
Data' Rate . Requirement . in Table 1 • 1J2 Cl ) J Applications : 
Coastal-Zone . Studies and Global Oceanography. 
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figure 4.4.2(8)B Number of IC's vs. Launch Date to Implement the- Micropro 
cessor Table . Look-Up (pPTLU) Processor for the Minimum 
Data Rate ■- Requirement in Table 1.1. 2(1).'- Applications: 
Coastal-Zone Studies ' and Global Oceanography '. 




Figure 4. 4. 2 (9) A Number of IC’s ys. Launch Date to Implement the Micropro 

cessor Table Look-Up. (yPTLU) Processor for the Maximum 
Data Rate -Requirement in Table 1.1. 2 (.1).'" Applications: 
Geography,' Hydrology and Meteorology. - 
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Figure 4.4.2(9)B Number. of IC’s .vs; .Launch Date' to Implement ^ the Micropro 

cessor Table LooK-Up..(iiPTLU) .Proces'sor for the' Minimum 
Data Rate Requirement • in Table 1 . 1 . 2 ( I ) , Appllcajrions ; 
Geography, Hydrology and Meteorology. 
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Figure 4.4,2G10)A Number of IC^s ys^ Launch Date. to Implement ^ the' Hardware 

Table Loqk-UpCHTLU) Processor for the . Maxiitiuin Data Rate 
Requirement in Table 1.1. 2 Cl).' Applications: Agricul- 
ture, Forestry and Geology. 
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pTo;ux«e M-.4.2Cll)A Number of IC’s vs. Launch Date to Implement the Hardware 

Table Look-Up (HTLU) . Processor for the Maximum Data Rate 
Requirement . in . Table 1.1. 2 Cl). Applications; .Coastal- 
Zone N Studies . and ' Global ■ Oceanography. 
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Figure 4.4.2(11)3 Number .of JC’s vs. Launch Date. to Implement < the Hardware 

Table Look-Up ‘.‘(HTLU) Processor for .the. Minimum Data. Rate 
-Requirement- in. Table' 1.1. 2(1). Applications: Coastal- 

, Zone Studies and -Global Oceanography. 
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Figure 4. 4. 2(12) A 


Number’ of lC*s vs. Launch Date to Implement > the Hardware 
Table Look-Up :(HTLU) .Processor for .the Maximum Data -Rate 
Requirement 1 in • Table 1.1, 2(1). Applications: Geography, 
Hydrology and' Meteorology. 
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Table' 4.4. 2(V) For Each. Application Listed in the First Column j the' Succeed- 
ing Columns List ,th.e Year that the . processor Becomes Feasible 
C 2 OOO iC’s) for the’ Jiaximum Requirements Listed. in Table 
l.l.'2Ci)J N means Not Feasible' by, 1990.. ‘ 


Application 

Micro- 

processor 

Maxi'mum 

Likelihood 

(UPML) 

Hardware 
Maximum 
■Likelihood 
■ CHML) 

Micro- 
processor 
•Table 
Look-Up . 

■ '(pPTLU) ■ ' 

Hardware 
Table 
•Look-Up 
' • ' (HTLU) 

A1 

1988 

1980 

N 

N' 

A2 

1988 

1980 

N 

N 

A3 

1988 

1980 

N 

N 

A4 

N 

1981 

N 

N 

A5 

1988 

1980 

N 

N 

Cl 

N 

1983 

N 

N 

C2 

N 

1983 

N 

N 

C3 

N 

1983 

N 

N 

C4 

N' 

1981 

M 

N 

C5 

N 

1981 • 

N 

N 

C6 

N 


N 

N 

C7 

N 

1984 

N 

1984 

FI 

1985 

1980 

N 

N 

F2 

N 

1983 

N 

N 

F3 

N 

1984 

N 

N 

F4 

1988 

1980 

N 

N 

G1 

1986 

1980 

•1988 

1980 

G2 

1986 

1980 

1988 

1980 

G3 

1983 

1980 

1986 

1980 

G4 

1989 

1981 

N 

1982 

G5 

1989 

1981 

N 

1982 

LI 

1984 

1980 

N 

1983 

L2 

1984 

1980 

N 

1983 

L3 

1984 

1980 

N 

1983 

L4 

1984 

1980 

N 

1983 

L5 

1986 

1980 • 

N 

1984 

L6 

1981 

1980 

1990 

1980 

HI 

1983 

1980 

1980 

1980 

H2 

1985 

1980 

1981 

1980 

H3 

1986 

1980 

1983 

1980 

H4 

1985 

1980 

1981 

1980 

H5 

1982 

1980 . 

1980 

1980 

Ml 

1980 

1980 

1980 

1980 

M2 

1980 

1980 

1980 

1980 

01 

1980 

1980 

N 

N 

02 

N 

1984 

N 

N 

- 03 

1980 

1980 

N 

N 
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Table 4.4.2(VI) For Each, Application Listed in the First Column, the Suc- 
. ceeding Columns • List the. Year .that the Processor Becomes 
Feasiblh.C2000 IQ' s) -.for the ^linimum Requirements Listed in 
Table 1.1.' 2Cl) 5 H. -means Not Feasible by 1990,- •' 


Application 


Micro- 

processor 

Maximum 

Likelihood 

(pPML) 


Hardware 
Maximum 
Likelihood 
• (HME) 


Micro- 
processor 
Table 
Look-Up ' 
(pPTLU) 


Hardware 

Table 

. . . Look-Up 
• ' (HTLU) 


A1 

1984 

1980 

1987 

1980 

A2 

• 1984 

1980 

1937 

1980 

A3 

1984 

1980 

1987 

1980 

A4 

1983 

. 1980 

1985 

1980 

A5 

1984 

1980 

1987 

1980 

Cl 

1985 

1980 

N 

1990 

C2 

1985 

1980 

N 

1990 

C3 

1985 

1980 

. N 

1990 

C4 

1982 

1980 

H 

1986 

C5 

1982 

1980 

K 

1986 

C6 

1988 

1980 

N- 

N 

C7 

1980 

1980 

N 

1982 

FI 

1980 

1980 

1983 

1980 

F2 

1985 

1980 

1988 

1980 

F3 

1986 

1980 

1988 

1980 

F4 

1983 

1980 

1986 

1980 

G1 

1983 

1980 

1986 

. 1980 

G2 

1983' 

1980 

1986 

1980 

G3 

1980 

1980 

1983 

1980 

G4 

1985 

1980 

1987 

1980 

G5 

1985 

1980 

1987 

1980 

LI 

1980 

1980 

1980 

1980 

L2 

1980 

1980 

1980 

1980 

L3 

1980 

1980 

1980 

1980 

L4 

1980 

1980 

1980 

1980 

L5 

1980 

1980 

1980 

1980 

L6 

1980 

1980 

1980 

1980 

HI 

1980 

1980 

1980 

1980 

H2 

1980 

1980 

1980 

1980 

H3 

1980 

1980 

1980 

1980 

H4 

1980 

1980 

1980 

1980 

H5 

1980 

1980 

1980 

1980 

Ml 

1980 

1980 

1980 

1980 

M2 

1980 

1980 

1980 

1980 

01 

1980 

1980 

1980 

•1980 

02 

1980 

1980 

1983 

1980 

03 

1980 

1980 

1980 

1980 
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4^,5 POSSIBLE EFFECTS OF NASA STIMULUS TO INDUSTRY * 

Landsat-1 and 2 remote-sensing satellites have effectively demonstrated 
the value of monitoring of the earth *s processes and resources • Earth has 
only a limited ability to support life^ and we must continuously determine 
our effect on the environment. To be effective, this monitoring must be 
done on a long term basis and, while the rewards will be great, the costs 
will not be cheap • On-board processing would be very important to reduce 
the flow of data to the ground stations and consequently reduce the costs. 
Under practical cost and power considerations, the technology is not here 
today to build and fly a practical system, although one could be flown to 
satisfy a few users. 

It is reasonable to assume that the technological development of cer- 
tain critical components would hasten the application date of a cost-effec- 
iiive system. NASA funds could be set aside for the development of critical 
components of the on-board processor. It is the recommendation of this 
study for NASA not to stimulate Industry through increased funding in those 
areas which are receiving a natural stimulation. As an example, we feel 
that it would be unwise for NASA to undertake the development of faster 
memory, logic modules and microprocessor chips for on-board processors. 

There is already considerable economic pressure on industry to develop these 
items and improvements in their performance will occur without NASA’s aid 
to industry. This philosophy would allow NASA to spend its resources in 
areas which receive little stimulus from outside sources and yet are vital 
to the continuation of NASA programs. The special technologies discussed 
in Section 3.6 are example of the type of technology which would benefit the 
greatest through NASA stimulation and support ^ 
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5 SUMMARY, CONCLUSIONS AND KECOMMENDATIONS 

This section contains an overview of the entire study reported in 

detail in the preceding sections* Significant results and conclusions are 

/ 

. discussed', and necoimnendations- for future actions by NASA are presented; 

Section 5.1 defines the study objectives, gives a description of the 
tasks, and a description of the significant results of each of the tasks. 

Section 5.2 describes the over.all conclusions resulting from the total 

study. 

Section 5.3 contains a number of recommendations to NASA as a result 
of this study. 

5.1 SUMMARY 

In this section we give a brief description of the problem of on-board 
earth resources data processing. The problem can be subdivided into a number 
of tasks, some of which are interdependent. The pr>ocedure for carrying out 
each of the tasks and the significant results are discussed. 

5.1.1 Study Objective 

Most of the past effort in the field of earth resources data processing 
has been research oriented. Earth resources- imagery has been pro- 
vided by NASA to a number of researchers who have processed the data in various 
ways in order to determine what, if any, useful information could be extracted 
from it. These experiments have demonstrated that useful information can in- 
deed be extracted from aircraft and satellite multispectral scanner imagery 
of the earth's surface. Economic studies have indicated potential cost 
effective systems based on these techniques. Consequently, it is anticipated 
that during the 1980-1990. decade earth resources satellites will be designed 
and flown for specific pmpposes, i.e. , to monitor severe weather systems, to 
monitor water pollution, to stirvey and monitor world food production, etc. 

In these applications it may be more cost effective to process the data on- 
board the satellite and transmit the data products directly to the users 
rather than transmit the raw data to a ground processing station for genera- 
ting the data products and then distributing the data products to the users 
via another satellite system. 

The pin?pose of this study was to investigate the feasibility of an on- 
board earth resources data processor launched during the 1980-1990 time 
frame. Since about five years are required to design, build, check out, and 
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launch such a system^ a 1980 system would be based on 1975 technology ^ and* a 
1990 system would be based on 1985 technology^ 

In order to determine the feasibility of on-board processing we must 
first define the on-board processor. This requires that we define 
both the' technology available for use in the design and the computational 
requirements required of the processor. The computational requirements . 
depend on the* - algorithms that the processor must implement which in 
turn, depend on the data products that must be extracted from the data to 
satisfy the users. Consequently, in order to determine the feasibility of 
on-board data processors we must start with a study of projected user appli- 
cations to define the data format (data throughput rate, number of spectral 
bands, etc.) and* the information extraction algorithms the processor must 
implement. Based on these constraints and the constraints imposed by the 
available technology we can design some on-board processors and evaluate their 
feasibility. The study plan is summarized in Figure 5. 1.1(1). 



Figure 5. 1.1(1) Study Plan 
5.1.2 User Applications Survey 

A number of studies to define earth resources applications areas, their 
data requirements, priorities, etc., have recently been completed. The major 
source and reference documents include; 

OA — Earth Resources Programs Summary 

Results and Projected Applications-ERTS-1 Applications and Investi- 
gations 

TERSSE volumes 1 to 8 
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Advanced Scanners and Imaging Systems for Earth Resources 
EOS - Paylord Discussion Group Report 
EOS - Mission Review Group Report 
Interplan Cost Benefit Study 
Dynatrehd Cost Benefit Study 
EARTHSAT Cost Benefit Core Studies 
From these references and from a number of personal interviews with 
researchers at a number of government and university laboratories, we deter- 
mined eight potential classes of earth resources data users: Agriculture, 

Coastal-Zone Studies, Forestry, Geography, Geology, Hydrology, Heterology, 
and Global. Oceanography, and a number of- subclasses within each of these 
classes. For each potential user we determined the- minimum and the maximum 
resolution, the minimum and maximum field of coverage, the minimum and maxi- 
mum number of spectral bands, and finally, based on these, the minimum and 
maximum data rate out of the' multispectral scanner. These results are tabu- 
lated in Table 1.1. 2(1). The resolutions range from a minimum of 3 m to a 
maximum of 10 km. The fields of coverage range from 15-800 km. The number 
of spectral bands ranges from 1-20 and the resulting data rates range from a 
minimum of 312 bits /sec to 3470 megabits/sec. 

Since the data requirements for the various users cover such a wide range 
wC' selected a single candidate data format for the initial effort. This candi- 
date data' format has a swath width of 185 km, a resolution of 40 m, a satel- 
lite ground track velocity of 6500 m/sec, 7 spectral bands, and 6 bits per 
data word. The resulting data rate from the MSS is 32 megabits/sec. This 
data rate satisfies all but two of the minimum data rates and about half of 
the maximum data rates suggested by the users. 

5.1.3 Data Analysis Algorithm Survey 

Almost all of the data users surveyed indicated that their objectives 
could be satisfied using spectral signature analysis. Consequently, a detailed 
survey was made of algorithms for classifying n-dimensional vectors into one 
of M categories or classes, where n the number of spectral bands. As a 
result of this survey it was determined that four algorithms warranted de- 
tailed analysis. These are clustering, maximum likelihood, perceptron and 
table look-up. 

Clustering is an unsupervised data analysis technique used to determine th 
natural or inherent data classes in a set of -observations. Many such algorithms 
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have been siiudied. All essentially make a scatter plot of a siibset of the 
data to determine the different groupings within the data. Each group is 
assigned a label, and all of the data with this label are compared to ground 
truth to associate each label with one of the classes defined by the data 
user. After this training is completed, each data point is classified by 
measuring the distance between it and each of the cluster centers and class i- 
:^ing it according to the nearest cluster. 

The maximum likelihood algorithm is a statistical proceedure based on 
the probability density function of the data. For the case of Gaussian data, 
which is a valid approximation for multispectral imagery of the earth's 
surface, only first and second order statistics are required. Training is 
accomplished by measuring these statistics for some data samples from knoTO 
classes and then assuming that all data from the same class 'has these same 
statistics. Subsequently, data are classified by comparing their statistics 
to the statistics of each of the classes and deciding in favor of the class 
they most closely resemble. 

The perceptron algorithm is based on a set of decision functions which 
are adjusted by an iterative technique to fit training data of known class and 
then used to classify subsequent data. 

The table look-up algorithm essentially stores in a large table (compu- 
ter memory) all possible outcomes of the data and associates with each possible 
outcome one of the classes.. Training is required to associate one of the 
’ classes with each of the possible values of the input data. Subsequent data 
are -then classified by using the data point to address the memory to look up 
the classification. 

The clustering, maximum likelihood and perceptron algorithms require a 
significant amount of computation, mainly additions, multiplications and com- 
parisons. The table look-up algorithm requires a much smaller amount of 
computation, but significantly more memory. 

5.1. 4 Preprocessing Algorithms 

The extremely large volume of data generated by the MSS imposes a severe 
computational burden on the on-board processor. The possibility of using a 
preprocessor between the sensor and the processor to reduce the bulk of data 
by using data compression, feature selection, etc.,’ techniques was studied. 

Transform coding allows a data bulk reduction by a factor of 2 to 4 
•for most multispectral data without degrading the data quality. 
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The BLOB algorithm developed at Purdue Universiry achieves data bulk- 
reduction by a factor of 10 to 30, but requires more computation and more 
memory than transform coding. 

5.1.5 . Algorithm Computational Requirements 

Each of the -data analysis algorithms- (clustering, maximum likelihood, 
perceptron, and table look-up) and the preprocessing algorithms (transform 
coding, and BLOB coding) were analyzed in detail relative to their computa- 
tional requirements, i.e., the number of additions, multiplications, compari- 
sons, etc., required to implement these algorithms along with the requirements 
imposed by the sequence of operations (some operations can be done in parallel 
while others follow a sequence vfhere one operation must be completed before 
the next can begin). These algorithm computational requirements were tabulated 
for each of the data analysis and preprocessing algorithms. 

Using a preprocessor to reduce the load on the processor is not a lucra- 
tive alternative. Even though the preprocessor can reduce the data load by a 
factor from 2 to 30 and thus reduce the complexity of the data processor by 
this amount, the total system complexity is not reduced because the savings 
in processor complexity are more than overwhelmed by the increase in the pre- 
processor complexity. 

It was further determined that the perceptron and clustering algorithms 
require a more complex processor than the maximum likelihood and table 
look-up algorithms for all user requirements. Consequently, we concluded 
that only the maximum likelihood and table look-up algorithms are worthy of 
further consideration. 

5.1.6 Technology Forecast and Assessment 

A detailed survey of 1975 component technologies was completed. A 
number of 1975 microelectronics technology families are listed in Table 3. 2.1(1) 
The speeds, power, size, cost, reliability, etc., of each are tabulated. 

Component technology was also projected from 1975 to 1985 using esti- 
mates obtained from component manufacturers and other experts in the field. 

The major conclusions are that some parameters associated with microelectronic 
component technology are changing at rates between 1 or 2 orders of magnitude 
every 10 years, with the result that overall component performance is changing 
by sever^ orders of magnitude every 10 years. In particular, the number of 
components (gates, transistors, etc.) per chip increased by a factor of' 10 
be-fcween 1965 and 1975 and is expected to increase by another factor of 10 
•between 1975 and 1985. In addition, propagation delays decreased by one order 
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of magnitude between 1975 and 1985, and are expected to decrease by another 
order of magnitude between 1985 and 1995. With the equivalent number of gates 
in an IC chip increasing by- a factor of 10 and the processing speed increasing 
by a factor of 10, the total number. of computations per unit time Ccomputational 
power) increases by a factor of 100. . 

Projections for computer ' system technology resulted in similar estimates, 
i. e., microcomputer cycle times', add times, etc., are projected to decrease 
by one order of magnitude dviring the next 10 years as they have for the past 
10 years. The number of bits of memory contained in a given area on an IC 
chip are likewise projected to increase by an order of magnitude over the next 
10 years as they have over the past 10 years. Meanwhile, the size and power 
dissipation, per IC chip is expected to stay constant while the number of pins 
per package which increased by a factor of four between 1965 and 1975 is 
expected to increase by only a factor of two between 1975 and 1985. 

We also developed a computer model that uses input data from past years 
to predict future values of these parameters. These computer generated projec- 
tions are in close agreement to the predictions made by experts from the 
microelectronics industry. 

5.1.7 On-Board Processor Designs 

A nuniber of on-board processors capable of implementing the maximum 
likelihood and table look-up algorithms for the candidate input data format 
were designed. In order to operate in real time at the 32 megabit/sec data 
rate the designs are based on multiprocessor concepts using pipeline and 
parallel arrays of subprocessors. Sufficient subsystems were added in parallel 
to obtain the 32 megabit/sec throughput. 

Two different design approaches were investigated in detail. One is a 
hardware approach consisting .of logic circuits designed to effi- . 
ciently implement the mathematical operations required by the algorithms. 

One special purpose hardware design was done to implement the maximum likeli- 
hood algorithm and another special purpose hardware design was designed to 
/ 

implement the table look-up algorithm. 

The second design approach uses microprocessors which allows a number of 
different computations to be done with the same hardware under software con- 
trol. Programs for implementing all of the computations were written in order 
to determine the number of instruction cycles required to implement the 
algorithm. This established the throughput data rate and, consequently, the 
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-number of parallel subsystems required to handle the 32 megabit/sec rate.'' 

Applying both of these design approaches to both algorithms resulted in 
four system designs. Hardware maximum likelihood (HML), hardware table look- 
up (HTLU) , microprocessor maximum likelihood CuPML) and microprocessor table 
look-up (pPTLU). . For each of these designs the number of IC's, power > volume, 
weight, and cost were determined based on 1975 technology. 

Because microprocessors are significantly slower than TTL circuits the 
hardware approaches require fewer IC*s, less povrer and volume, and cost less 
than the microprocessor designs. 

5.1.8 Feasibility Trade-Off and Sensitivity Analysis 

Each of the processor designs handle the 32 megabit input rate by distri- 
buting the processing load between many similar subprooessors . Consequently, 
the number of IC's, power, weight, volume; and cost are all essentially pro- 
portional to the number of subprocessors. Therefore, we defined a system 
complexity function for each of the four processors and determined its de- 
pendence on the following parameters using 1975 technology: 

R data bit rate (bits/sec) 
n nun^er-of spectral bands (channels) 

b number of bits per resolution element per spectral band (bits) 

M number of classes 

r pixel rate (resolution elements/sec) 

From the results- of the component and system technology forecasts we. further 
took into account the complexity function dependence on time for 10 years into 
the future. The resulting complexity functions are listed in Table 5. 1.8(1) 


Table 5.1. 8(1) Processor Complexity Functions 


Processor 

microprocessor maximum likelihood (yPML) 
hardware maximum likelihood (HML) 
microprocessor table look-up (yPTLU) 
hardware table look-up (HTLU) 


Complexity Function 
Pf = k^ M(ntl) Rd.S)""^ 

P^ = k^ M(n+1) RCl.5)"'^ 

Pg = kg M R (1.6)"^ (l-S)""^ 
= kj^ M R (1.6)”^ (1.5)"'^ 


The scale factors k^^, k^, kg, and k^^ were determined for each performance 
measure (number of IC’s, power, volume, weight and cost) and are listed in 
Table 4,3.3(11). 
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These models for the four design approaches were then used to determine 

the sensitivity of the complexity to the various system parameters ♦ This was 

accomplished setting all system parameters to their baseline values n = 

0 

M=12,b=6,R=32xlO , and T = 0 (the candidate baseline format) except 
for one or more parameters which were then varied from their baseline values. 
Figures 4.4(1), 4.4(2) and 4.4(3) show the sensitivity of the designs to 
variations in the data throughput rate R, the time T, the number of bits per 
data word b and the number of spectral bands n. 

Any feasibility analysis depends on the definition of feasible. For a 
particular processor to be "feasible" at a particular point in time requires 
that it meet certain constraints on performance, complexity, volume, weight, 
power, cost, reliability, environment, etc. Each of the four system archi- 
tectures meets the performance constraint since each was designed to accom- 
plish the required task. All four processors use standard integrated circuit 
technology and meet the data throughput rates by adding more components (IC's) 
in parallel. The volume, weight, and power dissipation of Integrated circuits 
can be kept within limits by simply keeping the number of integrated circuits 
within limits. The radiation, temperature, and other environmental constraints 
can be met by each processor as discussed in Section 2. The limiting factors • 
are cost and reliability which can also be kept within bo\mds by imposing a 
constraint on the number of components. Consequently, we conclude that on- 
board processing using a particular processor is feasible provided we constrain 
the number of IC’s in the processor to a reasonable number. 

The parts cost of the on-board processor increases linearly as the number 
of IC's. The costs associated with check out increase as the square of the 
number of IC’s. Limiting the number of IC’s in the on-board processor to 
about 1000 appears to satisfy all constraints, i.e., cost is reasonable rela- 
tive to the total system cost (launch, sensors, telemetry, etc.), reliability 
is-- pushing the limits of present day technology as discussed in Section 3, 
while volume, weight, and power dissipation do not appear to present serious 

♦ Figures 4. 4. 3(1) through 4.4.3(12) show the number of IC packages 
required to implement each of the four hardware designs as a function of time 
from 1980-1990. 

If we accept the 1000 IC definition of feasibility, then we can list 
the year in which each user application first becomes feasible for each design 
approach. These results asre summarized in Table 4. 4. 3(1) and 4. 4. 3 (2). Sum- 
mary tables corresponding to other definitions of feasibility could easily be 
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generated from Figures 4.4. 3(1) through 4.4.3(12). These results are sum- 
marized in more compact form in Figures 5.1. 8(1) through 5. 1.8(3), which 
show the percentages of user applications that can be implemented by each of 
the four design approaches for both the minimum and the maximum user require- 
— -ments . 
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Figure 5.1. 8(1) Percent of User Applications Satisfied as a Function of 
Time (1,000 IC's). 
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Figure 5. 1.8(2) Percent of User Applications Satisfiea as a Function of 
Time (500 IC*s). 
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Figure 5.1. 8(3) Percent of User Applications Satisfied as a Function of 
' Time (2,000 IC's). 
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5.2- CONCLUSIONS 


Conclusions of this study are: 

(1) From results of the user applications survey we conclude 
that potential users will require a wide range of resolutions, 
a wide range of field of coverage, a wide range of number of 

channels, and these in turn result in a wide range of data 
throughput rates. 

(2) From the results of the survey of data analysis algorithms 
we conclude that the maximum likelihood and table look-up 
algorithms are superior to other known algorithms for all 
user requirements. The table look-up algorithm is superior to 
the maximum likelihood algorithm, except for situations re- 
quiring more than five spectral bands^, 

(3) From the results of the investigation of the possibility of 
using a preprocessor to reduce the data load on the processor, 
we conclude that the total on-board system complexity is mini- 
mized with no preprocessor j 

(4) From the results of the component and computer system tech- 
nology forecasts and assessment, we conclude that the on- 
board processor capability (the amount of throughput it will 
be able to handle) will increase by two orders of magnitude 
between 1975 and 1985. 

(5) From the on-board processor designs and evaluations we conclude 
that implementations utilizing specially designed hardware 
require less hardware, power, volume, weight, and cost less 
than microprocessor (software) based systems. 

(6) From the feasibility and sensitivity analysis, we conclude 
that most but not all user applications could be satisfied by 
an on-board processor sometime between 1980 and 1990. 

5.3 RECOMMENDATIONS 

5.3.1 Handling the On-Board Processor Output Data Products 

While this study was directed towards determining the feasibility of 
on-board processors for the 1980-1990 time frame, the question remains as to 
how the output of an on-board processor could be treated. Now that this 
study has established the feasibility of on-board processing, the problem of 
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compressing and distributing the on-board computer output needs to be 
addressed. It is recommended that a study be made to investigate the uses 
of on-board processor output with particular attention paid to data rates 
and formats. 

5.3.2 Dates for Cost-Effective Launches 

Our study concludes that some users could be satisfied vfith a processor 
designed today and flown in 1980. Other users cannot be satisfied until 
1990 and beyond. These conclusions are based on> technical feasibility and 
do not address the question of economics. It is recommended, therefore, 
that a study be made to establish a projected time frame for the launch. of 
cost-effective earth resources missions. 

'5,3.3 Stimulation of Industry 

Finally, it is recommended that no stimulus be given to industry to 
develop large-scale integration (LSI) technology for on-board earth re- 
sources processors. It is recommended instead that the resources be used to 
encourage the solution of problems peculiar to NASA which do not have a 
parallel in industry, such as improvements in multi-spectral scanners and 
special techniques such as parallel processing. While this philosophy may 
be contrary to NASA's "spin-off" philosophy, stilulii from other sources 
are already present in the LSI field and further stimulus by NASA would 
have little effect and would be wasteful of NASA resources. 
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